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SECTION  I 


INTRODUCTION 


1.0  PROBLEM  DEFINITION 

High  performance  ballistic  vehicles  are  subject,  during  re-entry,  to  an  intense 
fluctuating  pressure  field  which  can  affect  the  integrity  of  the  vehicle  structure  and 
impose  adverse  vibration  levels  on  internal  components.  It  has  been  observed  that 
maneuvering  re-entry  vehicles  possessing  the  added  complication  of  a control  device 
(such  as  a flap)  experience  vibration  responses  which  exceed  levels  measured  for 
ballistic  vehicles.  The  pressure  fluctuations  arise  from  instability  and  unsteady  mo- 
tions of  fluid  flow  within  the  transitional/turbulent  boundary  layer.  Several  hypotheses 
have  been  proposed  to  describe  the  random  motion  of  the  momentum-deficient  fluid 
which  focus  on  intermittent  eruptions  of  the  viscous  sublayer.  Research  in  this  field 
has  been  conducted  both  theoretically  and  experimentally  with  emphasis  on  incompres- 
sible flow.  Pressure  fluctuations  classically  have  been  examined  by  relating  the 
phenomenon  to  velocity  fluctuations  through  Poisson's  equation.  Moreover,  further 
simplification  has  been  invoked  by  considering  only  an  interaction  of  the  turbulent 
structure  and  mean  shear  stresses.  In  this  manner,  mean  square  pressure  fluctua- 
tions for  attached  flows  have  been  predicted.  However,  the  phenomenon  is  still  under- 
stood only  vaguely,  and  design  criteria  have  been  developed  primarily  on  the  basis  of 
experimental  data. 

Hie  design  resolution  to  this  problem  can  be  divided  into  two  parts:  first,  providing 
a definition  of  the  fluctuating  pressure  environment,  and  second,  predicting  the  vehicle 
structure  and  internal  component  response.  The  present  investigation  will  focus  ini- 
tially on  the  first  part.  Empirical  correlations  developed  then  will  be  applied  to  the 
second  part,  in  particular,  to  a ballistic  and  a maneuvering  re-entry  vehicle.  Empirical 
modeling  of  aeroacoustlc  loads  for  axlsymmetrtc  or  two  dimensional  bodies  currently 
exists;  however,  adjustments  are  required  to  obtain  agreement  with  experimental 
data.  These  data  were  acquired  from  wind  tunnel  tests,  aircraft  flight  tests,  rocket 
and  missile  flight  measurements  as  well  as  limited  re-entry  vehicle  flight  measure- 
ments. One  should  keep  in  mind  that  the  aeroacousttc  environmental  empirical  models 
are  subject  to  the  following  constraints;  namely:  (1)  they  have  been  developed 


w^Wffnin^w.n 


primarily  from  data  on  non-ablating  two  dimensional  shapes.  (2)  the  models  do  not 
consider  flow  over  control  surfaces . and  (3)  they  require  further  modifications  to 
achieve  agreement  with  existing  experimental  data. 

2.0  OBJECTIVE 

The  objective  of  the  present  investigation  is  to  formulate  methods  for  predicting 
the  fluctuating  pressures  associated  with  maneuvering  type  vehicles.  An  assessment 
has  been  made  concerning  applicability  of  both  existing  data  and  present  analytical 
capability  to  fluctuating  pressure  characteristics  associated  with  vehicle  control 
surface  regions.  The  current  program  tncludes  a review  of  scaling  laws  utilizing 
normalization  of  aeroacoustic  parameters,  definition  and  implementation  of  a wind 
tunnel  test  program  to  acquire  data  that  will  verify  or  redefine  existing  prediction 
schemes,  and  application  of  upgraded  aeroacoustic  load  prediction  methods  to  a 
ballistic  and  a maneuvering  re-entry  vehicle. 
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SECTION  n 


CONCLUSIONS  AND  RECOMMENDATIONS 

In  general,  the  maneuvering  re-entry  vehicle  model  that  was  studied  in  this  investi- 
gation displayed  characteristics  that  were  approximately  the  same  as  a ballistic  type 
re-entry  vehicle  except  for  conditions  local  to  the  control  surface  region.  The  following 
algorithms  are  recommended  to  describe  the  aeroacoustic  environment  relative  to  rms 
pressure,  power  spectral  density  and  cross-spectral  coefficients  for  various  flow  and 
model  geometric  conditions. 

1.  0 ATTACHED  TURBULENT  BOUNDARY  LAYER  FLOW 
1. 1 Cone/Slice  Region  (<y  = 0® ) 

Root-Mean-Square  Pressure 

(°/qe)0  “ eT  <9°) 

where 


(a/qe)  - 0.006;  incompressible  value  (69) 

f - (T*/T  jCam-a+njVP+n) 

T*/Te  - [1/2  (1+  TwAe)  ♦ . 22r  ivIe  ] (56) 

for 


m ■ 3/2  + 


In  [(Te  + 198, 6)/(Tw  + 198.6)] 
In  fr^/T^ 


r ■ turbulent  recovery  factor  3 . 896 

n * velocity  power  law  exponent  * 7 for  fully  developed 
turbulent  flow 

n ■ 9 for  most  wind  tunnel  data  and  values  of  TBL  just 
after  transition 


(68) 


3. 


Power  Spectral  Density 


2 2 

cp(ou)ue  (2/tt)  (c/qeJj  eT 
q|  o*  }c“  1 + (l/e|  Xuu  6*/ue)2 


Cross-Spectra  Coefficients 
Cone 


(S.w)  ■ e 


■.1ou?/uc(uj) 


[.58  *.42  e--04865''e']  \ 


A CT.»J  - e-l.1«V^«»)[e-0.»ll>5V»‘].  „ s ,7e 


(94) 


Seotion  V 
2.3.1 


Slice 


At(S,u>) 


-e'-2u,5/u^i[.»5  + 


.SSe-*073^ 


(98) 


(99) 


1.2  Flap  Region  (0  * -6F  * 20) 

Root-Mean -Square  Pressure 

- (1  + I (o/<lo)  eT 

Power  Spectral  Penalty 

<p(u>)u0  fVn)^/^)2  (1+  I 6pt/90)2 

0 1 + (*  fiVue)2 


(91) 


(95) 


4. 


In  the  above,  it  is  noted  that  fimct.ional  variation  for  extended  flap  angles  is  not  con- 
sistent from  rms  pressure  measurements  to  those  for  power  spectral  density  i.  e. , 
the  exponent  for  PSD  values  was  anticipated  to  be  unity  (square  of  the  rms  pressure 
value).  However,  the  above  is  recommended  on  the  basis  of  limited  data  analyzed 
in  these  experiments. 

Cross-Spectral  Coefficients 

(S  ,uj)  = Equation  (98) 

(T) , u>)  = Equation  (99) 

2.0  TRANSITIONAL  BOUNDARY  LAYER  FLOW 

As  a consequence  of  the  voluminous  potential  of  acoustic  data  that  could  be  generated 
from  the  wind  tunnel  test  program,  a data  sub-set  was  selected  for  analysis.  The 
data  sub-set  was  established  on  the  basis  of  angle  of  attack  and  flap  deflection  combi- 
nations that  are  compatible  with  re-entry  vehicle  applications.  As  a result,  the  data 
sub-set  did  not  provide  sufficient  transition  information  that  could  be  used  to  interpret 
control  surface  characteristics.  Accordingly,  only  conical  frustum  transition  measure- 
ments will  be  presented  in  this  section. 

2. 1 Conical  Frustum 

Root-Mean-Square  Pressure 

l*£Jl±n) 

(°/<le)c  ■ («/<le>/[a/2)(l  + Tw/Te)  + .llr  (y-IJM*  ] (3+n)  (92),  (93) 

n » 6 peak  transition  values 
n » 4 for  start  to  middle  of  transition  zone 

Power  Spectral  Density 

An  attempt  was  made  to  incorporate  the  compressibility  factor  (c^.)  into  the  power 
spectral  density  algorithm  for  transitional  flow  as  developed  for  rms  pressure  with 
the  appropriate  value  of  n.  However,  no  direct  fluid  dynamic  logic  appeared  feasible. 

On  the  other  hand,  the  zero  intercept  values  of  PSD  appeared  to  be  predictable  using 


the  inverse  of  the  compressibility  factor,  i.  e. , 


tp(o)Ue 

q|  6* 


P/nKo/qJ2  /e£ 


and  should  be  used  accordingly  with  n = 4. 

Cross-Spectral  Coefficients 

A§(§,a))-e“0,06lu?/uc^)[.6  + .4  e',024?/6*],  all  values  of  a 


(97) 


(100) 


AT,Cn,aJ)-e-68uJf?/uc^)[ 


- nodn/f,* 


(101) 


windward 

Or  > 0, 7 3q 


(102) 


3.0  ANGLE  OF  ATTACK  EFFECTS 


Although  limited  angle  of  attack  data  were  analyzed  in  this  program.  It  was  deter- 
mined that  angle  of  attack  effects  did  exist  in  the  measured  data.  In  particular, 
normalized  acoustic  in  the  cone  and  flap  regions  exhibited  significant  a variations 
compared  to  the  slice  region  In  which  angle  of  attack  effects  appeared  to  be  compen- 
sated for  by  changes  in  the  compressibility  parameter.  Consequently,  It  Is  recommended 
that  cone  and  flap  rms  pressure  equations  for  windward  ray  environments  be  divided 

bv  the  function  (1  + a/9  ) l.e. 

c 


Cone 


(o/Vc 


Equation  190) 

a ♦ »/9J 

0 


(o/Vo 


Equation  191) 
(1  ♦ a/9c) 


Flan 


4.0  GENERAL  COMMENTS 


An  examination  of  the  fluctuating  pressure  data  indicated  that  overall  sound  pressure 
levels  (OASPL)  measured  on  the  slice  were  less  than  those  obtained  over  the  conical 
frustum.  This  is  believed  to  be  a result  of  the  expansion  wave  at  the  cone/slice  inter- 
section. The  condition  was  consistent  for  all  acoustic  functions.  In  turbulent  flow, 
the  OASPL  for  the  cone  was  approximately  5 to  10  dB  higher  than  slice  values.  With 
respect  to  the  control  surface  region  characteristics,  there  did  not  appear  to  be  any 
upstream  effects  on  the  slice  or  the  cone  resulting  from  flap  deflections.  For  de- 
flections up  to  20° , all  flap  characteristics  were  local  except  in  the  case  of  zero  de- 
flection for  which  measured  acoustic  functions  were  virtually  the  same  as  3lic&  values. 

At  Maoh  4 with  20*  flap  angle,  the  OASPL  experienced  by  the  flap  was  approximately 
20  dB  greater  than  the  slice  and  15  to  20  dB  greater  than  conical  frustum  values.  On 
the  other  hand,  at  Mach  8 conditions,  the  20*  flap  angle  showed  OASPL  levels  of  20 
to  25  dB  greater  than  slice  values  and  15  to  20  dB  greater  than  cone  levels.  When 
transitional  flow  was  experienced  along  the  control  surface  region,  for  Mach  4 tests, 
the  20*  flap  deflection  exhibited  levels  30  to  35  dB  higher  than  turbulent  levels. 

When  root- mean- square  fluctuating  pressure  and  power  spectral  density  are  non- 

_ . .2 

dimenalonalized  using  the  edge  dynamic  pressure  (i.e. , p/q  & cp(f)u  /q  6*),  measure- 
ments  from  the  cur  lent  program  exhibit  a significant  Mach  number  effect  which  is 
consistent  with  previous  wind  tunnel  data.  It  is  evident  that  values  of  rms  pressure 
and  PSD  so-normalized  generally  decrease  with  increasing  edge  Mach  number.  Lower 
values  of  these  functions  are  noted  both  for  the  higher  Mach  number  measurements 
of  Tunnel  B,  and  for  the  slice  data  which  exhibit  effects  of  increased  velocity  associated 
with  the  flow  expansion.  Mach  number  effects  are  explicitly  introduced  into  the  pre- 
diction formulae  for  rms  pressure  and  PSD  through  the  compressibility  factor  c^, 

whose  dependence  on  M is  described  in  Equations  (69)  and  (56). 
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5.  0 RE-ENTRY  RANDOM  VIBRATION- BALLISTIC  AND  MANEUVERING 
RE-ENTRY  VEHICLES 

Prediction  of  vibration  response  characteristics  for  a ballistic  re-entry  vehicle 
at  component  mounting  locations  (under  a specific  aeroacoustic  environment)  demon- 
strates overall  validity  of  the  present  analytical  technique  by  exhibiting  good  corre- 

2 

lation  with  maximum  response  levels  (.  022  g /Hz)  derived  from  applicable  flight  data 
in  fully  turbulent  boundary  layer  flow. 

Two  sets  of  maneuvering  vehicle  response  predictions  were  made  in  order  to 
demonstrate  how  much  effect  the  fluctuating  pressure  environment  associated  with  a 
.deployed  control  surface  has  on  structural  response  for  a specific  maneuvering  con- 
figuration. When  a hypothetical,  symmetric  acoustic  excitation  is  applied,  maximum 
internal  levels  of  . 044  g2/Hz  at  2013  Hz  and  . 028  g2/Hz  at  1657  Hz  result.  When  the 
environmental  description  is  modified  in  the  control  surface  region  to  account  for 
deployed  flap  acoustic  characteristics,  a significant  shift  in  contributing  modes  occurs. 
Maximum  response  levels  are  now  directly  associated  with  those  modes  in  whioh  there 
is  coupling  between  the  flap  and  other  vehicle  structure.  Maximum  internal  (compo- 
nent) levels  under  the  appropriate  maneuvering  R/V  loading  are  now  Increased  to 
. 066  g2/Hz  at  154  Hz  and  .063  g2/Hz  at  454  Hz,  with  significant  motion  evident  at 
flap,  R/V  nose,  and  shell  locations.  Maximum  responses  therefore  increased  only 
by  a factor  of  approximately  2. 4 when  modes  above  2000  Hz  are  not  included.  However, 
at  frequencies  for  which  there  is  significant  flap-vehicle  coupling  R/V  internal  accel- 
erations went  up  by  an  order  of  magnitude  or  more.  It  is  to  be  emphasized  therefore 
that  detailed  characteristics  of  "abide  dynamic  response  (i.e. , explicit  coupling 
behavior)  must  be  known  before  specific  conclusions  can  be  reached  regarding  the  in- 
fluence of  control  surface  exoltatlons  on  a maneuvering  R/V. 

Shell  acceleration  levels  (accounting  for  0*“  and  1st  harmonics)  under  the  cone/ 

slice/deployed  flap  environment  reached  maximum  levels  of  . 1 g2/Hz.  If  effeots  of 

st 

harmonics  higher  than  the  1 had  been  Included,  maneuvering  shell  response  would 

have  been  significantly  higher,  resulting  in  levels  comparable  to  actual  flight  data 
2 

(.  8 g /Hz  maximum).  However,  Inclusion  of  higher  harmonics  would  not  appreciably 
alter  predicted  internal  levels,  since  component  packages  ore  generally  mounted  on 
support  structures  such  that  their  primary  behavior  is  simple  beam-column  motion. 

The  support  structure  does  not  transmit  vibration  levels  associated  with  high  harmonic 


shell  response.  It  Is  therefore  to  be  concluded  that  measured  shell  accelerations  on 

2 

the  order  of  . 5 to  1. 0 g /Hz  are  not  inconsistent  with  internal  levels  (for  components 
on  shelves  or  bulkheads)  of  less  than  . 1 g~/Hz,  as  computed  for  the  specific  maneu- 
vering configuration. 

On  the  basis  of  flight  data  as  well  as  vibration  levels  predicted  herein,  it  is  possible 
to  specify  general  environmental  levels  for  ballistic  and  maneuvering  R/V  components. 
The  levels  discussed  below  are  strictly  applicable  only  to  the  ballistic  and  maneuvering 
R/V  configurations  considered  in  this  report;  other  vehicles  should  be  similarly  eval- 
uated before  applying  recommended  environments  to  additional  R/V  programs.  For 

2 

components  in  the  ballistic  re-entry  vehicle,  a maximum  level  of  . 05  g /Hz  would 
effectively  envelope  both  measured  and  computed  internal  responses.  Re-entry  vehicle 
shell  levels  of  .036  g /Hz  were  computed  in  the  present  analysis;  however,  highc  r 

harmonics  are  not  accounted  for,  and  shell  responses  are  consequently  underpredicted. 

2 , 

Flight  data  indicate  that  0. 1 g /Hz  is  an  appropriate  upper  bound  for  axial  and  radial 

shell  responses  applicable  to  the  ballistic  vehicle. 

For  maneuvering  vehicle  components  mounted  on  internal  structures,  the  present 

2 

analysis  and  data  included  in  Volume  n indicate  that  . 1 g /Hz  represents  a reasonable 

upper  bound  on  re-entry  vibration  level.  On  the  other  hand,  actual  flight  response 

measurements  on  maneuvering  vehicle  shell  structure  reveal  that  a range  of  approxi- 
2 

mately  0.2  - 1. 0 g /Hz  would  envelope  anticipated  vibration  environments  for  shell- 
mounted  items  in  maneuvering  configurations. 


SECTION  HI 


EVALUATION  OF  EXISTING  DATA  AND  PREDICTION  METHODS 
1.0  PREVIOUS  WORK 

An  inspection  oi  the  literature  has  revealed  that  considerable  work  has  been  de- 
voted to  incompressible  How  fields  and  more  recently  to  compressible  flow  in  the  super- 
sonic range  (M<5).  With  the  exception  of  the  limited  work  by  Heller  et  al1^’  12 ’ 13 

very  sparse  information  was  developed  relative  to  aeroacoustics  in  hypersonic  flow 
prior  to  the  AFFDL  sponsored  detailed  work  by  Chaump  et  al1  in  1972.  The  investi- 
gation of  Reference  1 is  considered  a comprehensive  review  of  surface  pressure  fluc- 
tuation characteristics  of  ground  and  flight  test  data  relative  to  overall  magnitude, 
spectra,  cross-correlation  functions  and  convection  velocity.  The  data  were  categorized 
into  various  boundary  layer  type  flow  conditions  including  attached  turbulent,  attached 
transitional,  separated  turbulent  and  base  flow.  In  addition  to  this  review  of  ground 
and  flight  data,  Reference  1 also  reports  on  a wind  tunnel  program  in  which  pressure 
fluctuations  were  measured  on  a 7. 2 degree  half  angle  cone.  The  test  program  con- 
sidered transitional,  turbulent,  separated  (induced  by  crossflow  on  a cone),  and  base 
flow  conditions  at  various  angles  of  attack,  Mach  number  and  nose  radii.  Measure- 
ments were  correlated  with  data  obtained  from  the  literature  to  develop  improved 
prediction  methods.  Consequently,  with  the  exception  of  attached  or  separated  flow 
conditions  over  a control  surface  (such  as  a flap).  Reference  1 is  an  ideal  starting  point 
for  reviewing  existing  data. 

A list  of  references  that  are  directly  or  indirectly  related  to  the  present  study  has 
been  compiled.  While  the  list  emphasizes  work  completed  in  the  present  decade,  many 
well  known  documents,  both  experimental  and  analytical,  have  not  been  Included  inas- 
much as  accessibility  and  general  contents  are  noted  in  the  papers  presented  in  the 
Reference  section.  The  reference  list  can  be  sub-divided  into  several  categories  which 
include  the  following: 

1)  General  Electrlc/AFFDL  sharp  cone  pressure  fluctuation  data  (References 
1 through  8).  It  should  be  noted  that  Reference  7 re-examined  the  spatial 
and  temporal  properties  of  the  data  of  Reference  1 while  Reference  8 was 
a direct  result  of  the  literature  review  of  the  present  shady. 

2)  NASA  Ames  flat  plate  work  of  Raman  (Reference  9). 
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3)  Bolt,  Beranek  and  Newman  work  by  Heller  et  al  on  a sharp  eone  and 
straightwbig  orbiter  configuration  (References  10  through  13). 

4)  The  NASA  Ames  work  of  Coe,  Chyu,  Hanly,  and  Dods  on  an  ogive- 
cylinder  and  wind  tunnel  walls  (References  14  through  20). 

5)  Flight  data  from  several  sources  (References  21  through  26). 

6)  Several  auxiliary  studies  relating  to  specific  aeroacoustic  data 
(References  27  through  35). 

7)  Flow  separation  and  pressure  gradient  effects  on  pressure  fluc- 
tuations (References  36-39,  also  References  14-19  and  28-29), 

8)  Boundary  layer  survey  (hot  wire  probes)  work  for  validation  of 
aeroacoustic  phenomena  (References  40-45)  and  freestream  noise 
effects  (References  32-34  and  95). 

9)  Analytical,  semi-analytical,  empirical  correlation  work  describing 
aeroacoustic  characteristics  (References  46  through  80). 

The  present  investigation  will  be  primarily  concerned  with  supersonic/hypersonic 

flow  conditions  and,  as  such,  will  not  be  concerned  with  a detailed  literature  review 

of  subsonic  flow.  Concerning  the  latter,  interested  readers  are  referred  to  the  ex- 
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cellent  survey  paper  by  Willmarth  as  well  as  Reference  9. 

1 . 1 Evaluation  of  Existing  Data 

When  assessing  data  reported  in  the  majority  of  references  cited  above,  it  was 
difficult  to  extract  raw  data  from  the  graphical  representations.  This  was  generally 
due  to  the  lock  of  available  aerodynamic  parameters  that  were  used  in  the  references 
for  normalization.  Moreover,  model  size,  geometric  construction  and  model  mount- 
ing techniques  rendered  very  questionable  data  results.  For  example,  data  employing 
freestream  Instead  of  local  boundary  layer  conditions  with  the  state  of  the  boundary 
layer  being  transitional  at  best,  must  be  carefully  weighed  before  being  used  in  fully 
developed,  hypersonic  turbulent  boundary  layer  predictions.  As  a result,  the  present 
study  has  concentrated  on  attached  flow  conditions  where  emphasis  has  been  placed  on 
the  work  associated  with  References  1,  9,  and  18.  Finally,  it  should  be  noted  that  the 
prediction  methods  developed  to  date  have  been  founded  on  turbulent  attached  flow 
criteria  where  modifications  have  been  made  to  consider  transitional  and  separated 
flow  effects.  Also,  the  data  sub-set  analyzed  from  the  experimental  task  of  this  study 


considered  attached  flow  behavior  only.  This  was  a consequence  of  using  realistic  de- 
sign criteria  (i.  e. , flap  angle  versus  angle  of  attack  variations)  for  maneuvering  re- 
entry systems. 

1.1.1  Power  Magnitude 

The  overall  magnitude  of  pressure  fluctuations  in  attached  turbulent  flow  from 
a number  of  experiments  is  shown  in  Figure  1.  Here  the  root-mean-square  (rms) 
pressure  normalized  b>  the  local  dynamic  pressure  is  shown  as  a function  of  local  Mach 
number.  A significant  variation  in  the  data  is  noted  when  plotted  in  these  coordinates. 

It  is  interesting  to  note  an  increase  in  normalized  rms  pressure  with  Mach  number 
(M  > 5)  for  various  experiments  except  for  the  data  point  represented  by  Reference  11. 
Moreover,  data  were  obtained  on  several  geometric  configurations  including  flat  plates, 
sharp  and  blunt  cones,  wind  tunnel  walls  and  an  ogive  cylinder.  The  use  of  local  bound- 
ary layer  conditions  is  considered  more  reasonable  for  extension  to  flight  application 
than  freestream  values.  In  experiments  involving  wind  tunnel  walls  and  flat  plates, 
local  and  freestream  conditions  are  identical.  However,  for  those  data  points  reflect- 
ing geometries  other  than  the  above,  no  attempt  was  made  in  the  current  study  to  deter- 
mine the  local  boundary  layer  properties,  and  reported  values  corresponding  to  free- 
stream conditions  were  used  accordingly. 

1.1.2  Power  Spectral  Density 

Figures  2 through  4 display  various  normalized  spectra  formats  for  wind  tunnel 
data  for  subsonic  and  hypersonic  conditions.  These  data  show  the  normalized  power 
spectral  density  as  a function  of  Strouhal  number  for  attached  turbulent  boundary  layer 
flow  conditions.  Several  authors  have  noted  that  the  spectra  can  be  graphically  dis- 
played as  a function  of  frequency  in  terms  of  the  Strouhal  number,  fl/v,  utilizing  charac- 
teristic length  and  velocity  parameters.  However,  choice  of  the  proper  typical  length 
and  velocity  has  been  a topic  of  discussion.  In  general,  the  freestream  or  local  bound- 
ary layer  edge  velocity  has  been  chosen  as  the  characteristic  velocity  parameter,  while 
the  boundary  layer  thickness  or  displacement  thickness  are  commonly  used  charac- 
teristic length  parameters. 

Figure  2 shows  the  dimensionless  form  of  power  spectral  density  (PSD)  where 
the  boundary  layer  edge  velocity  and  displacement  thickness  parameters  have  been 


chosen  for  normalization.  While  a general  collapse  of  the  data  is  evident,  no  dis- 
cernable  Reynolds  number  or  Mach  number  effects  are  apparent.  Figure  3 shows  the 
PSD  as  a function  of  Strouhal  number  where  freestream  velocity  and  boundary  layer 
thickness  are  used  for  normalization  parameters.  Data  in  Figure  3A  represent  the 
PSD  on  wind  tunnel  walls  while  Figure  3B  compares  tunnel  wall  data  with  those  ob- 
tained on  an  ogive-cylinder  at  supersonic  conditions.  It  should  be  noted  that  the  tunnel 
wall  measurements  agree  with  the  model  data  using  freestream  parameters.  This  is 
a consequence  of  the  flow  expansion  along  the  ogive  section  where  the  local  boundary 
layer  edge  values  approach  freestream  values  as  one  approaches  the  cylindrical  section 
(the  recording  station  X/D  = 2.9  is  approximately  at  the  intersection  of  the  ogive  and 

cylinder).  One  should  also  note  the  dramatic  increase  in  the  normalized  PSD  in  the  low 
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frequency  range  (f  6/u  < 10  ),  For  data  that  have  been  reported  in  the  literature,  a 
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significant  variation  exists  in  this  range.  Data  have  also  been  reported  that  tend  to 
decrease  at  Strouhal  number  < 10  *.  As  a result  of  the  disparities  in  this  region,  many 
experimentalists  have  not  reported  low  frequency  data. 

Figure  4 shows  a plot  of  the  PSD  distribution  with  Strouhal  number.  Here,  the 
dimensionless  spectra  are  normalized  by  the  ratio  of  the  square  of  rms  pressure  and 
dynamic  pressure.  The  authors  of  Reference  1 chose  this  format  because  of  the  ap- 
parent coalescence  of  the  data.  However,  when  data  of  recent  experiments  are  added, 
a significant  variation  is  observed.  Assessing  the  spectral  measurements  as  a function 
of  Strouhal  number,  it  is  apparent  that  the  dimensionless  form  of  PSD  decreases  with 
Strouhal  number,  in  particular  for  values  of  the  Strouhal  number  greater  than  10  1 
when  represented  in  the  coordinates  of  Figures  2 through  4.  It  will  be  shown,  in  effect, 
that  normalizing  one  dimensionless  quantity  with  another  is  actually  an  attempt  to  trans- 
form compressible  data  into  an  incompressible  plane. 

Data  describing  spectral  distribution  of  flight  measurements  are  very  limited. 
Inasmuch  as  ground  test  data  have  been  primarily  used  to  develop  prediction  capability, 
the  flight  test  data  will  be  used  only  when  deemed  appropriate.  Consequently,  these 
limited  data  will  be  discussed  when  prediction  techniques  are  applied  to  re-entry  vehicles 
(Vol.  H). 


1.1.3  Cross  Correlation  Functions 


References  describing  broad- and  narrow-band  cross  correlation  function  mea- 
surements for  attached  turbulent  boundary  layer  flow  conditions  have  been  compiled 
by  Howe  and  Chaump  et  al.  Measurements  in  the  supersonic/hypersonic  range  are 

very  limited.  It  should  be  noted  that  data  from  the  experiments  of  Heller  and  Holmes 
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and  Raman  were  not  included  in  the  investigation  of  Reference  7.  This  was  a conse- 
quence of  the  graphical  format  which  required  assumptions  for  explicit  use  (Reference 
10).  On  the  other  hand,  Reference  9 displayed  only  one  trace  (for  a given  Mach  number 
condition)  of  a common  space-time  and  auto-correlation  function.  This  investigation 
will  not  repeat  any  data  that  has  been  compiled  in  References  1 and  7 except  for  com- 
parison purposes  in  the  data  evaluation  section.  However,  considerable  attention  will 
be  given  to  the  data  interpretation,  in  particular  to  the  effects  of  fluid  flow  compres- 
sibility which  must  explicitly  be  accounted  for  at  hypersonic  conditions. 

1.1.4  Convection  V eloc  ity 

Figure  5 shows  broad-band  and  narrow-band  convection  velocities  for  both  in- 
compressible and  compressible  flow  data  obtained  from  a variety  of  geometries.  It  is 
quite  interesting  to  note  that  the  broad-band  values  appear  to  coalesce  for  both  the  in- 
compressible and  compressible  states.  Moreover,  the  incompressible  prediction  tech- 
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nique  of  Lowson,  when  modified  to  include  displacement  thickness,  adequately  de- 
scribes the  variation  of  the  broad-band  behavior  with  separation  distance.  On  the  other 
hand,  the  narrow-band  character  istics  reflect  a variation  between  Incompressible  and 
compressible  data  when  expressed  as  a function  of  frequency.  It  should  be  observed, 
however,  that  the  incompressible  data  indicate  a specific  trend  of  attenuation  in  the 
normalized  convection  velocity  with  decreasing  separation  length. 

The  characteristic  distance  for  spatial  normalization  has  been  chosen  as  the 
boundary  layer  displacement  thickness.  Keeping  in  mind  that  incompressible  flow  allows 
for  the  interchangeability  of  6 and  6*  (i.  e. , 6 - 86*),  one  would  expect  the  same  results 
if  6 were  used.  Some  further  comments  concerning  characteristic  lengths  will  be  given 
in  the  following  sections.  Finally,  it  is  interesting  to  note  that  Reference  7 presents 
compressible  data  for  angle  of  attack  conditions  (using  local  parameters  for  normali- 
zation) which  are  consistent  with  results  achieved  under  symmetric  flow  conditions 


(1.  e. , a - 0).  This  phenomenon  appears  to  be  valid  for  angles  of  attack  less  than  or 
equal  to  approximately  70%  of  the  cone  half  angle.  This  result  demonstrates  the  impor- 
tance of  using  local  boundary  layer  properties  for  normalization  as  opposed  to  free- 
stream  values. 

1.2  Evaluation  of  Prediction  Methods 

It  was  previously  noted  that  prediction  techniques  for  describing  aeroacoustic  loads 
have  been  developed  from  sound  theoretical  principles  appropriate  for  attached  turbu- 
lent boundary  layer  flow  behavior.  Modifications  were  made  to  the  attached  flow  methods 
to  account  for  flow  separation,  transition  and  base  flow  effects  primarily  on  the  basis 
of  experimental  data.  Such  an  approach  was  adopted  in  Reference  1,  which  considered 
a sphere-cone  type  configuration  (ballistic  R/V).  Consequently,  these  results  will 
form  the  reference  point  for  the  maneuvering  type  configuration  of  this  study.  This 
section  will  discuss  briefly  existing  prediction  techniques  for  describing  aeroacoustic 
environments  subject  to  attached  flow  conditions  for  ballistic  type  R/V's.  While  it  is 
not  the  intent  to  develop  a historical  review  of  the  pressure  fluctuations  in  turbulent 
flow  relative  to  prediction  capability,  an  examination  of  current  techniques  will  be  made 
regarding  the  fluid  dynamic  phenomenology  and  subsequent  assumptions  invoked  in 
developing  the  prediction  methods. 

1.2.1  Power  Magnitude 

Reference  1 outlines  the  correlation  for  predicting  the  magnitude  of  pressure 
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fluctuations  as  developed  by  Houbolt  that  was  modified  further  to  include  wall  tern- 
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perature  effects  in  unpublished  work.  A careful  examination  of  the  work  of  Lowson 
indicated  an  approach  similar  to  that  of  the  unpublished  Houbolt  concept,  bias  much  as 

the  Lowson  result  is  available,  details  concerning  its  development  will  not  be  given. 

0 

It  should  be  noted  that  a synopsis  of  the  Lowson  work  was  given  by  Laganelli  et  al. 

Houbolt  considered  a fluid  where  the  eddy  velocities  were  proportional  to  the  freestream 
velocities,  and  assumed  the  local  density  (region  of  maximum  noise  potential)  as  the 
significant  variable  governing  noise  production.  The  region  of  maximum  noise  pro- 
duction is  that  which  corresponds  to  maximum  shear  flow  or,  by  Reynolds  analogy, 
where  maximum  temperature  occurs.  If  one  conaiders-the  momentum  deficit  of  the 
fluid,  this  phenomena  occurs  within  the  classic  law  of  the  wall  region. 


15. 


Houbolt  considered  the  divergence  of  Euler's  equation,  where  the  properties  con- 
sist of  a sum  of  the  mean  and  fluctuating  components,  such  that 
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where  D/Dt  is  the  substantial  derivative.  The  above  expression  is  recognized  as  a 
Poisson  type  equation  with  the  right  hand  side  representing  a source  or  sink,  but  with 
steady  pressure.  Equation  (1)  has  the  solution 
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where  rQ  is  a distance  to  a general  field  point.  The  above  is  recast  into  dimensionless 
form  by  introducing  the  boundary  layer  thickness,  and  a local  mean  density  (p^)  at  a 
position  in  the  boundary  layer  where  the  maximum  noise  production  occurs  (6^).  Equa- 
tion (2)  is  then  written 
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where  6,  6 and  u are  scalar  quantities. 
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The  mean  and  fluctuating  velocities  were  considered  to  be  proportional  to  the  free- 
stream  velocities  such  that  the  density  profile  remains  fixed;  hence,  the  terms  in  the 
bracket  remain  invariant  to  tbs  flow  velocity.  The  pressure  field  can  then  be  described 
as 
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which  Incorporates  the  far  field  velocity  (u  ) as  well  as  the  governing  density  (p  ) and 
location  (6^)  at  the  site  of  intense  eddy  for -notion.  While  Houbolt  recognized  that  6/6^ 
can  be  expressed  as  a function  of  Mach  number,  its  dependence  is  weak.  Consequently, 
consideration  was  given  only  to  the  influence  of  p, . The  rms  value  of  pressure  can 
then  be  expressed  as 
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where  q is  the  dynamic  pressure  and  c a constant.  Equation  (5)  is  the  starting  point 

® 55  A 

of  the  analysis  given  by  Lowson. 

For  an  adiabatic  flow,  one  can  define  the  recovery  temperature  as 


T«A.  - 1 + r <ar->  “2e 


(6) 


and,  in  an  analogous  fashion,  it  is  assumed  that  the  temperature  at  the  site  of  maximum 
noise  generation  can  be  defined  by  a recovery  type  factor 
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where  the  coefficient  a /a  is  analogous  to  the  recovery  factor,  r,  of  Equation  (6). 
If  we  consider  the  equation  of  state  together  with  the  boundary  layer  assumption 
dP/dy  = 0,  Equation  (7)  becomes 
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and  the  rms  pressure  is  expressed  as 
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In  order  to  evaluate  the  new  recovery  type  factor,  a^/ a^ , the  Crooco  linear  tem- 
perature-velocity relationship  is  used;  namely, 


T/Te  - Tw/Te  + (l-Tw/Te)u/ue  + (^)M2e  (l-u/ue)u/ue 
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For  adiabatic  wall  conditions  (T  =T  ) about  the  position  6, , the  above  becomes 

w aw  1 

T/re  “ 1 + ^ J M^e  W (11) 

A comparison  of  Equations  (7)  and  (11)  allows  for 
a^  * unity 

a2/al  " (1 " Wfr  + U/V  (12ft) 


For  a non- adiabatic  condition,  the  above  becomes 
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Lowson  considered  the  velocity  ratio  uVu  to  be  1/2  and  the  constant  c to  be 
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0. 006  which  was  experimentally  determined  from  subsonic  flow  data.  The  resulting 
rms  pressure  was  then  expressed,  using  Equations  (5)  and  (11),  as  follows: 


(prms/<»oomp-0-00e''<1  + 0-14Me) 
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On  the  other  hand,  Houbolt  considered  the  constant  c^  to  be  0. 007  and  found  the  new 
recovery  type  factor  (a  /a  ) to  be  0.  06  as  a conservative  choice  when  considering  the 
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general  Crocco  equation  over  a range  of  wall  temperature  ratio  T /T  (4  to  7)  and 
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velocity  ratio  u/u  (0. 4 to  0. 7).  The  rms  pressure  was  then  expressed  as 
6 
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which  when  generalized  to  include  wall  temperature  ef  fects,  gas  density,  and  composi- 
tion becomes 
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where  * a^a^  and  is  given  by  Equation  (12b). 
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If  one  considers  the  definition  of  the  dynamic  pressure  (q  = y /2  p M ) the  rms 

G 6 6 

pressure,  when  normalized  by  the  static  value  for  the  Lowson  and  Houbolt  results, 
becomes  respectively: 

prm/pe  >L  ~ 0042  **Va  + »• « ) (16) 

and 
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At  this  point,  several  factors  concerning  the  above  formulations  should  be  noted. 

These  are:  (1)  the  constant  c^  (.  006  or  . 007)  was  obtained  at  through  experimental 

observations  of  subsonic  flow  conditions,  (2)  the  compressible  state  had  been  developed 

by  use  of  a density  formulation  only,  (3)  the  use  of  the  Crocco  linear  temperature- 

velocity  relationship  is  questionable  in  turbulent  flow  (for  example,  data  have  indicated 

a quadratic  relationship  exists);  (4)  the  velocity  ratio  (u  /u  ) of  1/2,  while  considered 

as  a representative  average  of,  say,  the  viscous  sublayer  velocity  to  freestream  value, 

can  vary  in  the  law  of  the  wall  region  where  turbulent  intensity  prevails;  and  (5)  use  of 

adiabatic  wall  (hot  wall)  conditions  is  not  justified  for  ail  cases,  in  that  many  flights/ 

experiments  were  conducted  under  cold  wall  conditions. 

While  Lowson  considered  compressibility  through  a density  change,  he  also  exam- 
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ined  the  works  of  Bies  for  predicting  t through  heating  near  the  wall  whose  key 

w 

effect  results  from  viscosity  rather  than  density.  Good  agreement  was  noted  between 
the  two  methods;  however,  one  must  keep  in  mind  that  data  are  limited  to  the  low  super- 
sonic range  where  the  comparisons  were  made.  Moreover,  Lowson  noted  that  T /T 

W fllW 

would  probably  not  be  less  than  0. 8 in  practice;  hence,  effects  would  be  insignificant 
for  his  adiabatic  wall  prediction.  However,  it  should  be  noted  that  the  condition 
T /T  < 0. 8 can  exist  on  re-entry  vehicles  and,  for  that  matter,  in  wind  tunnel  tests 

W fiW 

where  T /T  can  be  less  than  0. 5 especially  at  hypersonic  conditions. 

W BW 

Figure  1 shows  the  prediction  methods  of  Lowson  and  Houbolt  compared  to  data. 

It  is  noted  that  the  adiabatic  wall  results  of  both  authors  appear  to  give  the  best  agree- 
ment to  the  measured  experimental  data. 


9 

Raman  considered  several  correlations  of  the  rms  pressure  using  the  dynamic 
pressure,  wal'  shear  stress  and  static  pressure  as  normalization  parameters.  Poly- 
nomial or  exponential  type  fits  to  his  data  indicated  that 

p ■ f (q  § M );  f (M  , Re.J;  f_  (Re.) 

Inasmuch  as  these  correlations  were  not  developed  on  the  basis  of  fluid  dynamic  prin- 
ciples and  no  comparisons  to  other  data  were  made,  the  functional  representations  of 
Raman  will  be  used  when  deemed  applicable  to  the  present  study.  It  should  be  noted 
however,  that  the  functional  form  representing  Prms/(100  can  be  derived  through  fluid  dy- 
namic principles. 

If  one  considers  the  concept  subsequently  developed  in  Section  3.0,  Equation  (59), 
it  is  easily  shown  that 


P /q 
rms  e 
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From  boundary  layer  theory  for  a flat  plate 
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Re4/5  - 2.275  e~l/4  Re^4 


(19) 


where  e_  is  defined  by  Equation  (57).  The  ratio  of  displacement  thickness  to  momentum 

g3 

thickness  is  given  as 


6*/e  - -1  + C(^)TwAaw  + 1 3 [1  + 0.208  M2  ] 


(20) 
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If  one  considers  A to  be  equal  to  2,  an  average  of  the  range  predicted  by  Lilley, 
together  with  a 1/7  velocity  power  law  (i.e. , n = 7),  Equation  (18)  subject  to  Equations 
(19)  and  (20)  will  have  two  ranges.  These  are  the  cold  wall  case  and  the 

adiabatic  wall  case  (Tw*T&w).  Hence,  Equation  (18)  becomes 
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and 


prms/qe  >aw  " °- 052  * . 48  M2  ] * . ,3  M2  ]4/5 


(22) 


In  Equations  (21)  and  (22),  the  exponent  of  the  displacement  thickness  Reynolds 
number  (-1/4)  is  a consequence  of  a 1/7  fully  developed  turbulent  velocity  profile 
power  law.  A comparison  to  the  data  of  Reference  9 indicates  that  the  data  are  over- 
predicted for  Mod=  5.2  and  7.4  and  underpredicted  for  the  = 10.4  case.  It  should 
be  noted  that  the  cold  wall  prediction  was  used  inasmuch  as  the  adiabatic  case  was 
significantly  underpredicted.  The  slope  of  the  Raman  data  indicated  a 1/5  distribution 
which  is  synonomous  with  a 1/9  velocity  profile,  a condition  indicating  a non-fully 
developed  turbulent  boundary  layer. 

1. 2. 2 Power  Spectral  Density 

The  mean  square  fluctuating  pressure  in  terms  of  the  power  spectral  density 
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function  is  defined  as 
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where  the  mean  square  value  is  equal  to  the  total  area  under  a power  spectrum  corres- 
ponding to  a sample  time  history  record.  For  attached  turbulent  boundary  layer  flows, 
power  spectra  have  been  found  to  scale  with  Strouhal  number.  The  Strouhal  number 
represents  a normalized  frequency  using  characteristic  length  and  velocity  of  the  flow 
field.  The  choice  of  the  characteristic  length  and  velocity  parameters  has  been  a sub- 
ject of  considerable  debate;  in  particular,  for  high  speed  applications.  This  topic  will 

be  further  considered  in  a following  section. 
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Lowson  and  Houbolt  ’ developed  empirical  representations  of  the  power 
spectrum  on  the  basis  of  subsonic  and  available  supersonic  flow  data.  Koubolt,  apply- 
ing engineering  logic,  considered  a system  of  rolling  eddies  to  characterize  the  con- 
vective velocity  which  are  located  a distance  above  the  wall  at  a site  of  intense  noise 
distribution  to  formulate  the  Strouhal  number.  It  was  determined  that  spectra  appeared 
to  scale  in  terms  of  rms  pressure  as 
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Houbolt  further  noted  from  the  definition  of  the  rms  pressure  and  PSD  the  following 
relationship 
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where  dZ  = dui/ai.  Hence,  a plot  of  the  density  distribution  of  the  power  magnitude 
together  with  the  product  oucp(o>),  namely 


n.  uc  uo cp  (a.) 

2 6*  p3 

rms 


(26) 


would  maximize  the  power  magnitude  when  the  Strouhal  number  was  unity.  It  should 
be  noted  that  the  above  concept  is  quite  similar  to  the  first  moment  of  the  power  mag- 
nitude when  employing  the  definition  of  the  latter;  i.  e. , 
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This  approach,  which  is  similar  to  that  of  Black,  " is  essentially  a weighted  residual 

technique  and  wili  be  further  discussed  in  the  following  section. 

55 

Lowson  developed  an  expression  similar  to  Equation  (24)  that  considered  sub- 
sonic data  and  the  supersonic  data  of  Reference  37.  It  was  determined  that  the  data 
appeared  to  coalesce  when  the  Strouhal  number  was  based  on  boundary  layer  thickness 
rather  than  displacement  thickness.  The  expression  developed  by  Lowson  is  given  by 

2 o 3/2 

Vn(lJ,"trms/Vri  + 1 (28) 
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where  * Su^/6.  Robertson  re-examined  the  correlation  developed  by  Lowson 
when  comparing  the  above  formulatioi  to  supersonic  data  obtained  in  References  18 
and  19.  It  was  found  that  the  Lowson  prediction  underestimated  spectral  levels  at  low 
Strouhal  numbers.  As  a consequence,  Robertson  modified  the  Lowson  result  to  include 
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boundary  layer  displacement  thickness  and  freestream  velocity  and  also  changed  the 
exponents.  Robertson's  result  is  given  as 


„ Prms/q» 


(29) 


where  uu  =l/2u  /6*. 
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Figures  2 and  6 show  spectral  data  of  several  experimenters  compared  to  pre- 
diction methods  of  Houbolt  and  Lowson,  respectively.  There  does  not  appear  to  be  any 
advantage  in  using  either  concept  (i.e. , displacement  thickness  or  boundary  layer 
thickness).  It  should  be  noted  that  the  convection  velocity  in  Equation  (24)  was  replaced 
by  the  local  boundary  layer  edge  value.  Moreover,  the  Mach  number  effect  represented 
in  these  figures  was  a result  of  the  use  of  Equation  (13)  for  both  methods.  It  is  also 
noted  that  neither  method  agrees  well  with  the  data  over  the  spectrum  range.  One 
feature  of  the  data  variation  is  the  apparent  Reynolds  number  effect. 
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Raman  found  that  the  PSD  distribution  for  the  Mach  numbers  of  his  experiment 
followed  the  expression 


9 (QUqo 
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(30) 


where  Aj  is  a constant.  The  above  expression  is  recognized  as  the  formulation  by 
Lowson55  with  slight  modification.  Raman  integrated  the  above  and  noted  that  the  left 
hand  side  contained  the  definition  of  the  rms  pressure,  such  that  over  the  limits  of 
Integration  the  following  results: 
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Equation  (30)  is  then  written  as 
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As  previously  indicated,  the  advantage  to  using  the  above  format  for  displaying  PSD 
distribution  is  to  incorporate  compressibility  effects  with  an  incompressible  algorithm. 
Moreover,  Equation  (30)  pre-supposes  a solution  of  the  PSD  as  a function  of  Strouhal 
number  with  the  arbitrary  constant  , a methodology  suggested  by  Houbolt.  It  should 
be  noted  that  the  left  hand  side  of  Equation  (30)  considers  the  dynamic  pressure  as 
opposed  to  the  rms  pressure  for  normalization.  Manipulation  of  the  constant  A^  togeth- 
er with  the  definition  of  the  rms  pressure  yielded  the  result  obtained  by  Lowson  and 
Houbolt. 

1„2.3  Cross-Correlation  Functions 

Because  of  the  importance  of  the  narrow-band  spatial  correlation  function  to 
describe  the  impinging  effect  of  a fluctuating  pressure  field  on  a structure  and  its 
response,  a detailed  examination  of  the  field  must  be  made  at  many  points.  Moreover, 
the  complexities  of  describing  cross -correlation  and  cross -spectrum  functions  for  an 
aeroacoustic  environment  (equivalent  representations  related  through  their  Fourier 

transforms),  has  required  a semi-empirical  approach  coupled  with  physical  reasoning. 
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One  such  method  was  developed  by  Lowson  and  later  modified  by  Robertson  * 

and  Chaump  et  al. 1 The  analytical  development  of  Lowson  was  primarily  based  on  the 

incompressible  experiments  of  Bull. 80  Coe  et  al14  considered  a different  approach 

whereby  an  attenuation  coefficient  function  was  introduced  that  considered  a decaying 

exponential  of  the  moduli  of  the  cross-spectra.  A complete  discussion  concerning  the 

analytical  development  and  subsequent  modifications  of  the  cross  spectra  coefficients 

7 

is  given  by  Howe.  The  following  is  concerned  with  the  application  of  analytical  tech- 
niques developed  from  incompressible  flow  behavior  to  compressible  flow  situations. 
Consequently,  some  repetition  of  the  work  from  Reference  7 will  be  required  for  con- 
tinuity. 

Inasmuch  as  the  criteria  developed  by  Lowson  was  based  on  sound  mathematical 
and  physical  principles,  its  foundation  together  with  the  data  of  Bull  will  be  emphasized 
here.  Lowson  found  it  mathematically  convenient  to  assume  a product  solution  of  the 
correlation  coefficients  where  the  spatial  dependence  was  separated.  In  this  manner, 
the  real  contribution  of  the  complex  cross-power  spectral  density  function  (f^)  re- 
lating points  1 and  2 assumed  the  form 
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1/2 
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where  cp  and  cp  are  the  power  spectral  densities  at  locations  1 and  2,  respectively. 
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Lowson  allowed  for  a correction  to  the  separable  solution  by  a factor  of  tt/2  to  com- 
pensate for  possible  underestimation  of  the  correlation  area.  This  was  a consequence 
of  more  probable  elliptical  contours  of  the  data  about  tie  origin  as  opposed  to  the 
diamond  patterns  suggested  by  a separable  form  of  the  equations.  The  separable 
cross -spectral  density  factors  and  A^  (i.  e. , correlation  coefficients)  were  further 

separated  into  exponentials  of  spatial  and  mixed  spatial/frequency  terms;  namely 

A (s.ou)-e'cl,s,/6YC2,s,uj/uc  (33) 

s 

where  s is  a generalized  separation  distance  (T|  or  § ) and  and  c ^ are  empirical 
constants. 

(a)  Longitudinal  Correlation  (Incompressible) 

Lowson  considered  the  work  of  Bull  to  obtain  the  spatial/frequency  asymptotic 
high  frequency  exponential  function  of  Equation  (33),  shown  in  Figure  7.  Moreover, 
one  notes  a divergence  in  the  narrow-band  correlation  function  at  the  low  frequency 
end  of  the  spectrum.  To  account  for  this  divergence,  the  data  were  extrapolated  to 
the  ordinate  (<u  - 0)  to  obtain  the  asymptotic  values  of  the  coefficients  at  constant  § . 
These  asymptotic  values  are  shown  as  a function  of  %,/ 6*  in  Figure  8.  Also  shown  is  a 
curve  fit  to  the  data.  The  empirical  formula  representing  the  incompressible  data  of 
Bull  was  expressed  by  Lowson  as 

I A§|  " eXp  t "{ (°.  hi)  ?/uo)2  + (0. 27  S/6)2}  1/2 ] (34) 

It  should  be  noted  that  Lowson  considered  the  boundary  layer  thickness  rather  than 
displacement  thickness  for  normalization.  For  a fully  developed  turbulent,  adiabatic 
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flow  with  a 1/7  velocity  power  law,  the  ratio  of  boundary  layer  and  displacement  thick- 
ness 6/6*  is  8,  hence  the  change  in  constants  shown  in  Figure  8 (i.e. , 0.034  5/6*  = 

0.27  |/6). 

28 

Robertson  considered  a modified  version  of  Equation  (34)  in  the  form 

(§^)»e'0,1“^uOe'°,275/6  (35) 

which  was  also  adopted  by  Chaump  et  al. 1 If  one  were  to  plot  both  Equations  (34)  and 
(35),  a negligible  difference  results.  Equation  (35)  is  considered  a more  readily  usable 
formulation. 

Before  proceeding,  some  further  discussion  concerning  the  Bull  data  is  in  order. 

Relative  to  the  longitudinal  correlation  function,  Bull  notes  that  the  data  of  Figure  7 

exhibit  a coalescence  of  all  values  of  5/6*  at  high  frequency  which  suggests  identifying 

frequencies  with  convected  wave-numbers  according  to  the  relation  u>  = k_u  (a>).  Hence, 

5 c 

the  longitudinal  high  wave -number  components  will  lose  coherence  in  times  which  are 
proportional  to  the  times  required  for  them  to  be  convected  distances  equal  to  their 
wavelengths.  On  the  other  hand,  the  amplitude  of  the  correlation  coefficients  tends  to 
be  independent  of  frequency  for  low  values  of  Strouhal  number  at  a given  5/ 6*.  Accord- 
ingly, Bull  hypothesizes  a rapid  loss  of  coherence  in  a period  of  the  development  for 

the  initial  wave-number  component  where  similarity  of  uj5/u  is  required.  Once  these 
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high  wave-number  components  have  lost  coherence,  the  remaining  large  scale,  low 
wave-number  components  of  the  field  are  still  correlated  until  complete  development 
of  the  field  entails  their  gradual  loss  of  coherence.  Here,  the  general  shape  of  the 
spectrum-producing  components  remain  fairly  constant  in  the  process.  Moreover, 
this  narrow-band  asymptotic  trend  was  found  to  be  consistent  with  broad-band  cross- 
correlation  results. 

In  order  to  justify  the  above -postulated  asymptotic  state,  Bull  noted  the  low  fre- 
quency divergence  occurs  at  higher  values  of  u>5/uc  for  corresponding  higher  values 
of  5/6*.  This  implies  that  the  divergence  occurs  at  approximately  the  same  dimension- 
less frequency  w 6*/u  for  all  value3  of  5/‘*»  thereby  yielding  a characteristic  Strouhal 
c 

number  and  scale  consistent  with  the  final  phase  of  evolution  of  components  in  the  pres- 
sure field.  To  find  this  scale  in  the  pressure  field,  associated  with  low  frequency 
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turbulence,  one  considers  values  of  id?  /u^  at  which  the  asymptotic  high  frequency  curve 

(e  ‘ *u>  ?/uc)  has  a value  equal  to  R (5, 0, 0)  for  a given  value  of  ?/6*  to  obtain  the 

PP 

corresponding  value  of  u>  6*/u  . This  Strouhal  number  based  on  6 * is  considered  a 

c 

measure  of  the  largest  value  associated  with  the  large-scale  pressure  eddies.  With 
reference  to  Figure  7,  for  the  case  in  which  Mo=  0.  5 and  ?/6*  = 19.7b,  the  high  fre- 
quency component  of  the  low  frequency  asymptotic  value  of  R > 0.  5 is  <d?/u  5-7.2. 

PP  c 

The  corresponding  Strouhal  number,  based  on  6*  becomes 

5*/u  (cu)  ss  0.36 
c 

It  was  found  that  the  Strouhal  numbers  for  the  range  of  displacement  thickness  of 

Bull's  investigation  were  nearly  constant  with  an  average  value  of  approximately  0.36. 

Hence,  the  large  scale  pressure  eddies,  with  wavelength  \ in  the  stream  direction, 

are  those  greater  than  \^/6*  = 2rr/0. 36  = 17.4.  If  we  again  consider  the  incompressible 

ratio  6*/ 6 = 1/8,  the  wavelengths  would  be  greater  than  about  twice  the  boundary  layer 

thickness.  On  the  other  hand,  the  amplitude  of  the  high  frequency  response  appears  to 

reach  a value  of  0. 05  when  id  ?/u  is  24.  5.  Here,  the  component  has  been  convected 

c 

a distance  24. 5/2n  c-  4.  Consequently,  high  frequency  components  in  incompres- 

sible flow  appear  to  lose  their  identity  when  convected  a distance  approximately  4 times 
their  wavelength. 

(b)  Lateral  Correlations  (Incompressible) 

The  lateral  cross-correlation  coefficients  were  found  to  have  a more  simple 
decay  characteristic.  Figure  7 shows  the  asymptotic  high  frequency  envelope  which 
can  be  expressed  by  the  curve  e * ^ ^ ^Uc.  As  in  the  longitudinal  case,  the  data  re- 
flect a divergence  at  the  low  frequency  end  of  the  spectrum.  Lowson  considered  the 
extrapolated  amplitudes,  shown  plotted  on  Figure  8 from  Bull  results,  and  suggested 
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a curve  fit  to  the  data  in  the  form  e . Again,  using  the  incompressible  relation- 
ship 6*/ 6 equal  to  1/8  (such  that  e 2T^6  = e ,25T^'  6 ),  this  exponential  form  of  Lowson 
is  shown  compared  to  the  Bull  data.  It  is  quite  apparent  that  the  asymptotic  form  of 
the  data  for  u>  - 0 does  not  decay  as  suggested  by  the  data  fit  of  Lowson.  It  should  be 
noted  that  Robertson  considered  the  same  form1  as  Lawson,  which  was  believed  to  be 
valid  for  both  subsonic  and  supersonic  speeds.  However,  this  condition  will  be  shown 
to  be  invalid  for  compressible  data.  Also  shown  in  Figure  8 for  the  lateral  narrow- 


band  coefficient  is  the  expression  proposed  in  Reference  1 which  provides  a much 
better  fit  to  the  data.  Accordingly,  the  expression  representing  the  lateral  cross 
spectrum  coefficient  for  incompressible  flow  becomes 
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As  in  the  longitudinal  case,  the  scaling  for  the  high  and  low  frequency  contributions 

to  the  correlation  coefficients  are  desired.  The  amplitude  for  the  lateral  correlation 

coefficient  at  a value  of  0.  05  is  u)  n /u  =4.7  which  indicates  a convected  distance 

c 

0.75.  Hence,  a high  frequency  component  with  a given  longitudinal  wavelength 
is  laterally  coherent  over  a distance  of  approximately  1 1/2  times  its  wavelength.  In 
Figure  8,  one  notes  that  the  asymptotic  values  of  the  lateral  cross-correlation  function 
(as  a)  — 0)  did  not  decay  in  a pure  exponential  characteristic  as  might  be  expected.  By 
considering  the  low  frequency  correlation  amplitudes,  one  can  find  the  corresponding 
points  on  the  high  frequency  curve  (e  ' 72uj  ^ ^UC)  where  values  equal  to  Rpp(°»  "H » °) 
represent  the  smallest  of  the  low  frequency  pressure  eddies  in  which  the  correlation 
amplitude  is  independent  of  frequency  (see  dotted  line,  for  example,  in  Figure  7).  The 
average  value  of  the  Strouhal  number  was  determined  to  be  0. 36  as  in  the  longitudinal 
case.  Here,  it  is  important  to  note  that  the  low  frequency  longitudinal  and  lateral 
correlation  amplitude  characteristics  are  due  to  the  same  source  of  eddies  within  the 
pressure  field,  i.  e. , the  pressure  source  with  longitudinal  wavelengths  greater  than 
twice  the  boundary  layer  thickness. 

The  above  methodology  can  be  extended  to  surfaces  that  are  not  flat,  in  particular 
to  structures  that  are  axlsymmetric  or  feature  irregularities  such  as  control  surfaces. 
However,  additional  complexities  must  be  accounted  for,  such  as  those  effects  arising 
from  a non-symmetric  transition  front  as  developed  on  an  R/V  at  angle  of  attack.  In 
general,  any  surface  obstruction,  mode  of  ablation  or  angle  of  attack  situation  that 
will  produce  a non-uniform  flow  pattern  over  a surface  should  reflect  different  charac- 
teristics in  the  low  frequency  range. 


(c)  Compressible  Effects  (Longitudinal  Correlation) 
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As  previously  indicated,  Howe  re-examined  the  spatial  and  temporal  properties 
of  fluctuating  pressure  data  obtained  in  Reference  1.  Figures  9 and  11  show  the  nor- 
malized longitudinal  Co  and  Quad  functions  of  the  cross  spectral  data  of  Reference  7. 
The  data  are  compared  to  the  exponential  high  frequency  asymptotic  envelope 
e O’  1 (2n)f5/uc  wjj£C|j  characterized  the  incompressible  data  of  Bull.  ^ Howe  noted 

that  at  finite  values  of  ? the  asymptotic  amplitude  for  f§/u  -0  should  be  < 1.  Con- 

c 

sequently,  measured  cross-PSD's  should  be  compared  to  an  exponential  envelope  which 

has  a factor  to  account  for  spatial  attenuation.  Accordingly,  a factor  of  0. 8 or  0. 9 

was  incorporated  in  the  figures  to  account  for  measured  attenuation  corresponding  to 

§ » 0.  5 inches.  He  attributed  the  apparent  increasing  correlation  for  |f/u  > 1. 5 to 
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possible  tunnel  tare  noise,  as  seen  in  Figure  9.  Moreover,  at  the  higher  Mach  number 
conditions,  the  cross-spectral  properties  appear  to  exhibit  a greater  attenuation  com- 
pared to  the  incompressible  asymptotic  high  frequency  response. 

Figure  12  shows  the  normalized  Co  asymptotic  values  [$_(S,0,f-0)]  as  a function 
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of  §/ 6*  at  Mach  4 and  8.  As  in  the  case  of  Lowson,  these  values  were  obtained  by 

extrapolating  the  amplitudes  to  §f/u  = 0.  The  data  are  shown  compared  to  the  asymp- 
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totic,  incompressible  equation  developed  by  Lowson.  It  is  quite  apparent  that  the  data 
do  not  exhibit  a pure  exponential  trend  as  in  the  incompressible  case.  Howe  postulated, 
with  reservation,  that  the  high  frequency  components  could  be  characterized  by  the 
incompressible  exponential  envelope  but  modified  the  asymptotic  low  frequency  contri- 
butions (as  shown  in  Figure  12).  When  combining  frequency  and  spatial  effects,  he 
proposed  that  the  longitudinal  cross-spectral  coefficient  in  fully  turbulent  flow  takes 
the  form 


A?(S.f) 


#-O.Xu»5/uo*s)[0#5  + 0>5  e-0-07  5/6*^ 


(37) 


The  fact  that  the  compressible  data  show  a different  characteristic  is  not  at  all 
surprising.  The  effects  of  compressibility  on  the  fluctuating  pressure  field  should 
increase  the  wavelength  and  change  the  corresponding  small  and  high  scale  coherence 
characteristics. 
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As  an  indication  of  this  phenomena,  consider  the  methodology  employed  by  Bull 

where  a characteristic  non-dimensional  frequency  and  scale  can  be  obtained.  Here, 

we  consider  the  pressure  field  associated  with  the  low  frequency  asymptotic  state  for 

which  a Strouhal  number,  cug/u  , has  a value  at  which  the  high  frequency  asymptotic 

c 

curve  is  equal  to  $(§,0, f-*0)  for  given  values  of  5/**.  The  corresponding  wavelength 

can  then  be  obtained.  From  Figure  9,  consider  the  data  point  with  spatial  distance 

§ = 1.1  inches.  The  high  frequency  contribution  gives  ?f/u  ^ 0. 4 about  $ (§,  0,  f,  0)K-  0. 7. 

c 

The  Strouhal  number  associated  with  the  large  scale  eddies  becomes  u>6*/uc  = 0.25 
which  is  less  than  the  incompressible  value  obtained  by  Bull.  The  large  scale  pressure 
eddies  having  wavelength  will  be  greater  than  \^/ 6*  = 2rr/.  25  ^ 25.  For  the  Mach 
4 condition,  the  ratio  of  boundary  layer  thickness  to  displacement  thickness  is  */ 6*^  2; 
hence  the  wavelengths  would  be  greater  than  approximately  13  6,  a significant  departure 
from  the  incompressible  case. 

If  we  consider  the  asymptotic  high  frequency  amplitude  for  the  compressible 

state,  it  appears  to  be  characterized  as  in  the  incompressible  case.  From  Figure  9, 

9 (§,  0,  f)  - 0. 05  for  §f/u  - 4 which  corresponds  to  a value  of  §id/u  %25,  implying 
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that  the  pressure  component  has  convected  a distance  = 25/2n  ^4.  Again,  the 
high  frequency  component  appears  to  lose  identity  when  convected  a distance  approxi- 
mately 4 times  its  wavelength;  the  same  condition  experienced  in  the  incompressible 
case. 

From  the  above  it  appears  that  the  effects  of  compressibility  tend  to  decrease 
the  wave-number  whereby  the  low  frequency  (large  scale)  components  convect  more 
rapidly  in  the  stream  direction  than  the  corresponding  incompressible  result,  such  that 
they  lose  their  identity  even  more  slowly.  On  the  other  hand,  the  high  frequency  (low 
scale)  components  are  characterizeo  by  large  wave-numbers,  travel  slowly  in  the  stream 

direction,  and  subsequently  lose  their  identity  rapidly.  These  findings  are  consistent 
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with  the  hypothesis  of  Bull  and  the  compressible  experiments  of  Raman.  The  effect 

of  compressibility  is  quite  apparent  when  considering  the  high  frequency  decay  envelopes 

of  Figures  10  and  11.  A postulated  envelope  has  been  constructed  to  fit  the  data  at  Mach 

numbers  8 and  10.  It  is  clear  that  when  a characteristic  Strouhal  number  u>6*/u  is 

c 

obtained  corresponding  to  points  where  the  high  frequency  asymptotic  curve  has  values 
of  9(Z,0,f-  0),  (which  represents  the  smallest  of  the  low  frequency  pressure  eddies 
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for  that  portion  of  the  field  where  the  correlation  amplitude  is  independent  of  frequency) 

smaller  values  of  u>fi*/u  or  k 6*  (wave-number)  will  be  obtained.  This  in  turn  will 
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yield  high  values  of  wavelength.  It  is  interesting  to  note  that  such  a postulated  enve- 
lope yields  values  of  the  characteristic  Strouhal  number  (associated  with  the  large 
scale  pressure  eddies)  that  are  consistently  less  than  the  incompressible  values,  and 
subsequently  yield  wavelengths  an  order  of  magnitude  greater  than  6 for  the  compres- 
sible flow  case  compared  to  the  incompressible  flow  values. 

(d)  Compressible  Effects  (Lateral  Correlation) 

Figure  13  shows  a comparison  of  typical  circumferential  Co  cross-power  spectral 
density  for  incompressible  data  to  the  compressible  data  of  Reference  1 which  were 
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re-examined  by  Howe.  It  is  noted  that  the  incompressible  data  acquired  at  a given 
frequency  for  various  separation  distances  were  implicitly  attenuated  as  a result  of 
the  pure  spatial  dependence  suggested  by  Equation  (36).  Hence,  comparing  these  data 
to  a function  of  mixed  dependence  such  as  e c can  be  misleading.  Accordingly, 

Howe  compiled  the  Mach  4 and  8 data  of  Reference  1 at  several  frequencies  for  fixed 
separation  distances.  These  data  are  shown  in  Figure  14  where  one  notes  that  the  de- 
cay function  is  not  asymptotic  to  unity.  A comparison  to  Figure  13  indicates  that 
neither  exponential  form  shown  fits  the  compressible  data  of  Howe.  Consequently,  in 
Reference  7 a curve  fit  was  applied  to  the  compressible  data  and  is  shown  in  Figure  14. 

As  in  the  longitudinal  case,  variation  of  circumferential  cross-power  spectral 
data  can  be  determined  by  plotting  asymptotic  values  (f  - 0)  as  a function  of  separation 
distance  Tl/6*.  Figure  15  shows  the  results  for  Mach  4 and  8 conditions  from  Reference 
7.  Also  shown  is  the  low  frequency,  asymptotic,  incompressible  equation  suggested 
in  Reference  1.  It  is  apparent  that  the  compressible  data  exhibit  a higher  correlation 
than  that  implied  by  the  incompressible  data  fit  expression.  A curve  fit  to  the  com- 
pressible data,  developed  by  Howe,  is  also  shown.  Consequently,  if  one  considers  the 
high  frequency  asymptotic  amplitudes,  the  recommended  equation  for  circumferential 
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cross-power  spectral  density  for  a fully  developed  turbulent  boundary  layer  is 


A^cn.f)-  e 


■1.4wT|/uc(<JU) 


[0.3 


, n _ -0. 05571/6* 

+ °*7  6 3,  a £ 0.7  9 


(38) 
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The  lack  of  attenuation  in  the  circumferential  direction  of  the  cross-power  spec- 
tral density  indicates  the  strong  effect  of  compressibility.  As  previously  noted,  the 
low  frequency  longitudinal  and  lateral  (here,  circumferential)  correlation  amplitude 
characteristics  are  a result  of  the  same  source  of  eddies  within  the  pressure  field. 
Inasmuch  as  the  low  frequency  longitudinal  components  (large  scale)  lose  their  identity 
slowly,  the  same  phenomena  can  be  expected  in  the  lateral  direction  for  compressible 
flow. 

1.2.4  Convection  Velocity 

It  has  been  determined  that  the  convection  velocity  is  a function  of  frequency 
and  spatial  separation.  The  variation  with  frequency  is  a consequence  of  momentum 
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(different  size  eddies)  variation  of  the  pressure  field  in  the  boundary  layer.  Lowson 
noted  that  the  variation  with  spatial  separation  could  be  a result  of  accelerated  eddy 
trajectories  within  the  boundary  layer  as  well  as  varying  coherence  lengths.  This  is 
consistent  with  turbulent  boundary  layer  phenomenology  whereby  transition  and  turbu- 
lence are  a result  of  bursting  of  fluid  particles  from  the  viscous  sublayer  with  subse- 
quent exchange  of  momentum  with  the  surrounding  fluid.  Lowson,  recognizing  the 
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difficulty  of  defining  convection  velocity,  accordingly  accepted  the  definition  of  Bull. 
Moreover,  the  incompressible  data  of  Bull  were  used  to  develop  an  empirical  curve 
for  the  narrow-band  and  broad-band  convection  velocities. 

Bull  expressed  the  convection  velocity  as  a function  of  Strouhal  number  based 
on  displacement  thickness  wSVu^.  Lowson  obtained  the  broad-band  values  by  extra- 
polating broad-band  results  to  their  asymptotic  values  and  plotting  results  in  terms  of 
the  spatial  coordinate  5/6.  A cross-plot  of  the  Bull  data  at  a given  value  of  uu  6*/u 

CD 

along  lined  of  constant  5/6*  provided  r.he  convection  velocity  as  a function  of  5/6*. 
Empirical  curves  through  the  data  yielded  the  following  results 

«>)/«.-  0.675  + 0.3  e‘°*iiuj6/u«  _ 0.25  e~1,25/6  (39) 

where  the  broad-band  convection  velocity  curves  correspond  to  u>5/u  = 8.  such  that 

00 
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If  we  consider  an  adiabatic,  incompressible,  1/7  velocity  power  law  flow,  6/6*  = 8 
and  Equations  (39)  and  (40)  become 


u 
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(id)/ u ■ 0, 67  5 + 0. 3 e 
00 


•0,  88oi  6*/u0 


-0.25 


-0.15  5/6* 
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(41) 
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0. 25  e 


■0.15  5/6* 


(42) 


Equation  (42)  is  shown  to  adequately  describe  the  data  shown  in  Figure  5 for  both  the 
incompressible  and  compressible  flow  states. 

When  the  incompressible  formulation  of  Lowson,  based  on  6,  was  compared  to 
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the  compressible  data  of  Howe  a departure  was  observed  that  waB  not  reflected  when 
6*  was  used.  Figure  16  shows  the  normalized  broad-band  convection  velocity  as  a 
function  of  spatial  distance  5/6.  As  a consequence  of  this  disparity,  Howe  suggests 
a modified  version  of  the  Lowson  result  to  account  for  the  apparent  compressible 
effects.  However,  we  must  keep  in  mind  the  usage  of  boundary  layer  characteristic 
lengths.  In  Figure  16,  6 was  used  for  normalization  whereas  Figure  5 employed  6*. 
Again  if  we  consider  the  incompressible  relation  between  6 and  6*,  the  value  of  5/ 6 
where  the  incompressible  and  compressible  curves  join  in  Figure  16  is  about  8.  This 
value  corresponds  to  5/6*  equal  to  unity  in  Figure  5,  a region  where  the  data  show  a 
slight  deviation.  In  effect,  Figure  16  has  enlarged  the  scale  by  the  ratio  of  the  boundary 
layer  thicknesses  (i.  e. , by  a factor  of  8)  for  the  region  5/6*  s 1 of  Figure  5.  Further 
discussions  concerning  choices  of  normalization  length  and  velocity  parameters  will 
be  given  in  the  next  section.  It  will  be  shown  why  6*,  when  used  as  in  Equations  (41) 
and  (42),  is  a more  effective  parameter  than  6 when  considering  compressible  flow 
effects. 

Bull  noted  that  it  was  possible  to  assign  a unique  value  of  u = u (u  ) for  a given 
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value  of  the  Strouhal  number  uj  6*/u  , which  we  previously  indicated  had  near  constant 

ao 

values.  This  implies  that  the  fluctuating  pressure  components  have  characteristic 

velocities  at  some  distance  from  the  wall  with  a mean  velocity  of  u i (#»).  An  estimate 

c 

of  the  mean  velocity  was  taken  to  be  about  0. 6 u at  a position  from  the  wall  of  approxl- 
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mately  2%  of  the  boundary  layer  thickness.  This  corresponds  to  law  of  the  wall  coordi- 
nates y+-  yu  /u^  ^100  (see  sketch  in  Section  2.0).  Inasmuch  as  u^ (uj ) does  not  appear 

to  exceed  0. 9 u , it  was  postulated  that  the  eddy  system  whose  location  from  the  wall 

00 

exceeded  1/2  6 would  make  no  significant  contribution  to  the  wall  pressure  fluctuations. 

Consequences  of  the  above  observation  are  two-fold.  First,  it  establishes  a rela- 
tive position  in  the  boundary  layer  where  phenomenological  laws  can  be  used  to  develop 
prediction  methods  that  characterize  the  fluctuating  pressure  field,  as  for  example, 
in  References  8,  55,  56,  63,  and  67.  Second,  the  law  of  the  wall  coordinates  are  appli- 
cable to  both  incompressible  and  compressible  flow  conditions.  This  characteristic 
of  boundary  layer  behavior  could  be  used  to  interpret  the  coalescence  of  incompressible 
and  compressible  broad-band  velocity  distribution  with  spatial  distance  as  shown  in 
Figure  5. 

2,  0 EVALUATION  OF  AERODYNAMIC  NORMALIZATION  PARAMETERS 

The  empirical  basis  of  the  theory  of  turbulent  boundary  layers  has  hampered  de- 
velopment of  turbulent  aeroacoustic  technology  because  the  independent  variables  con- 
trolling production  of  turbulent  noise  cannot  be  derived  on  a strictly  theoretical  basis, 
but  must  be  inferred  from  data  correlations  or  from  a phenomenological  view-point. 

A number  of  empirically  developed  prediction  techniques  have  been  proposed  for  the 
calculation  of  turbulent  boundary  layer  induced  acoustic  levels  and  spectra.  As  previously 
noted,  the  choice  of  characteristic  length  and  velocity  associated  with  the  boundary 
layer  is  not  a standard  selection  process.  This  problem  is  further  compounded  for 
supersonic  flow  conditions  where  the  typical  characteristic  lengths  vary  significantly 
with  Mach  number  and  wall  temperature  ratio. 

Early  investigators  considered  the  boundary  layer  thickness  (6)  and  displacement 
thickness  (6*)  as  obvious  choices  of  the  characteristic  lengths.  In  particular,  the 
choice  of  the  latter  together  with  the  boundary  layer  edge  velocity  appeared  to  scale 
the  power  spectra.  The  use  of  6*  appears  to  be  satisfactory  for  both  subsonic  and 
supersonic  flow  conditions  as  determined  from  wind  tunnel  experiments  over  the  last 
decade.  However,  its  use  in  hypersonic  flow  conditions  (for  flight  applications)  must 
be  done  with  caution.  This  is  a consequence  of  the  sensitivity  of  6*  for  highly  cooled 
wall  and  pressure  gradient  regions  such  as  the  nosetip  where  the  displacement  thickness 
oan  become  negative. 
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In  an  early  paper,  Willmarth  recognized  that  the  power  spectral  density  qp  (<d) 
could  be  expressed  in  non-dimensional  form  as 


Cp  (ULI  )Ucd 

q2  6* 

CO 


F<V-*'  M-'  Re) 
06 


Several  investigators  have  used  the  rms  pressure  Prms»  or  wall  shear  stress 
rather  than  the  dynamic  pressure  q^ , and  5 or  6 (momentum  thickness)  rather  than 
6*.  The  above  can  be  re-written  as 

(Sf*> 


q2  6" 

*■  >^03 


rms  « 


wherep  /q  =f  (M,  Re,T  /T  , y ) and  is  consistent  with  the  findings  of  Willmarth. 
rms  ' w aw 

One  notes  that  the  functional  form  of  p _/q  , whioh  will  be  developed  in  Section  3. 0, 

rms  <*> 

depends  on  the  wall  temperature  ratio  as  well  as  properties  of  the  flow  environment. 
Moreover  the  above  representation,  as  previously  mentioned,  allows  for  a comparison 
of  compressible  data  with  incompressible  algorithms. 

Other  normalizing  characteristic  velocities  of  the  boundary  layer  consist  ol  the 
shear  velocity  u ( = /t  "/p ) or  convection  velocity  u rather  than  u . In  a recent 
paper,  Willmarth  noted  that  low  frequency  (large  scale)  wall  pressure  fluctuations 
scale  with  u and  6*  while  high  frequency  (low  scale)  fluctuations  scale  with  the  wall 
parameters  u and  v _ (kinematic  viscosity).  This  concept,  originally  proposed  by 
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Corcos  and  later  advanced  by  Black,  has  as  its  basis  the  fact  that  the  turbulent 

boundary  layer  is  not  a self-similar  flow,  in  the  sense  that  average  quantities  are  not 

functions  of  a single  non-dimensional  length  or  a single  non-dimensional  velocity. 

The  outer  part  of  the  flow  has  a characteristic  velocity  V (or  u for  flight  appli- 
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cat  ion/non -flat  plate  wind  tunnel  models).  The  characteristic  length  in  the  outer  wake 
region  is  the  boundary  layer  thickness.  Since  incompressible  flow  was  being  considered, 
6*  was  essentially  interchangeable  with  5.  In  the  law  of  the  wall  region  (that  region 
associated  with  the  site  of  intense  eddy  formation)  the  characteristic  length  suggested 
by  Corcos  is  u^/u^.  Corcos  noted  that  in  the  law  of  the  wall  region,  some  of  the  length 
scales  of  the  eddies  are  impressed  upon  it  from  the  outer  region.  On  the  other  hand. 
Black  showed  that  using  the  length  scales  of  the  outer  flow  correlates  the  low  frequency 
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power  spectra  while  failing  to  do  so  in  the  high  frequency  range.  The  opposite  was 
shown  to  be  true  in  the  use  of  the  inner  region  scales. 

To  further  support  the  concept  of  two  distinct  regions  characterizing  low  and  high 
frequency  scaling,  classic  boundary  layer  phenomenology  will  be  used  to  indicate 
choice  of  length  and  velocity  parameters.  If  one  considers  the  momentum  boundary 
layer  equation  for  a turbulent  flow,  there  results 


...  3u  du  dp  d 
PU  + PV  — “ - T*  + — (T  ) 

dy  r 


dx 


dy 


dx 


where  t is  the  combined  molecular  (laminar)  and  turbulent  shear  stress.  If  one  in- 
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troduces  the  Prandtl  mixing  length  concept,  the  shear  stress  can  be  expressed  as 


du 

ay 


t -puV  + Pll  (f:) 


2 .du  2 
m 'dy 


where  the  unprimed  terms  represent  the  mean  motion  (viscous  stress)  and  the  prime 


terms  the  turbulent  fluctuations  (Reynolds  stress  terms).  The  term  1 is  the  classical 

m 


mixing  length. 


Prandtl  considered  the  mixing  length  to  be  proportional  to  the  coordinate  from  the 
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wall,  such  that  1 = ky  for  k equal  to  0.41,  a universal  constant.  VanDrlest  ” argued 
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that  the  Reynolds  stress  is  of  considerably  greater  magnitude  away  from  the  wall  than 
the  viscous  stress,  such  that,  as  the  wall  is  approached  the  effect  of  the  viscous  stress 
should  start  to  become  more  significant  until,  at  the  wall,  viscosity  dominates.  Con- 


sequently, a damping  function  was  suggested  that  modified  the  mixing  length  to  give 
y/D, 


1 = ky  (1-e 

m 


).  The  term  D represents  the  properties  of  the  fluid  and  the  frequency 
of  oscillation  of  its  movement. 

If  one  introduces  the  following  dimensionless  groups 


u+  ■ u/u 


y ■ put/u 


then  the  mixing  length  becomes 


■;  -ky‘(i-.-yVD+) 


where  1 + ■ pu  1 /u  and  D+  - ou  D/u.  For  an  incompressible  flow,  the  constant  D+ 
m t m T 

has  been  determined  empirically  to  be  approximately  26.  For  compressible  flow. 
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Cebeci  has  modified  the  above  formulation  to  give 

1+  » k y+  (1  - e~yw/D  v 
m Jw  * • 


where 

+ + ] /o 

yw  - y %/uw and  D " 26  (fijzy 

The  universal  law  of  the  wall  plot  is  schematically  shown  below  together  with  a typical 
turbulent  velocity  profile. 


Cebeci  considered  two  regions  where  scales  of  the  eddy  terms  have  different 
characteristics.  The  two-layer  closure  model  has  been  compared  to  a number  of  ex- 
perimental data  indicating  feasibility  of  the  concept  for  predicting  turbulent  flow  charac- 
teristics. The  inner  region  is  characterized  by  an  eddy  viscosity  model  based  on 
Prandtl  mixing  length,  namely 
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and  represents  the  region  0 s y+  a;  700  on  the  above  sketch.  In  the  outer  region,  re- 
presented by  y+  > 700,  a constant  eddy  viscosity  is  used  that  considers  the  Klebanoff 
intermittency  factor  (y  ).  The  eddy  viscosity  for  this  region  given  by 

K 

emo“0-0168  VK  ®*ue 

where 

\K  " + 5-5  (y/6)6]"1 

Since  the  two-iayer  closure  model  of  Cebeci  is  synonomous  with  the  high  and  low 
frequency  regions  observed  from  acoustic  data,  an  analogy  between  the  turbulent 
closure  model  and  the  spectrum  of  the  fluctuating  pressure  field  can  be  made.  The 
mixing  length  will  be  the  characteristic  length  of  the  high  frequency  (low  scale)  region 
represented  by  the  law  of  the  wall,  including  the  viscous  sub-layer.  On  the  other  hand, 
6*  represents  the  characteristic  length  in  the  low  frequency  (high  scale)  region  of  the 
law  of  the  wake.  Moreover,  characteristic  velocities  are  the  shear  velocity  (u  ) for 
the  wall  region  while  boundary  layer  edge  velocity  appears  reasonable  for  the  wake 
region. 

The  characteristic  Strouhal  numbers  are  accordingly 


and 


)1/UT-  (D  U-w/ UT 


w 


low  scale 


u)  6*/ue  high  scale 

This  result  is  consistent  with  those  reported  previously  with  the  exception  of  Black89 
who  considered  ut  and  6 as  representative  velocity  and  length  for  the  low  scale  region. 

Non-dimensionalizing  the  pressure  fluctuations  at  the  wall  has  been  accomplished 
throughout  the  literature  by  using  local  dynamic  pressure,  local  static  pressure  and 
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wall  shear  stress.  Corcos  and  Laufer  note  that  normalizing  with  minimizes 

variation  of  the  results  with  Reynolds  number.  Moreover,  Laufer  reported  that  the 

radiated  energy  spectrum  contains  considerably  less  high  wave-number  components 

61 

than  the  pressure  spectrum  measured  at  the  wall.  Bull  proposed  the  same  phe- 
nomenological reasoning  which  was  implicitly  supported  by  Corcos'  and  Black's  view- 

46 

point  and  more  recently  by  Willmarth.  The  use  of  as  a normalizing  parameter 

for  the  fluctuating  pressure  has  theoretical  basis  as  proposed  by  the  early  works  of 
70 

Lilley.  This  approach  will  be  used  in  the  current  investigation  to  be  presented  in 
the  following  section.  Inasmuch  as  the  shear  stress  is  a difficult  parameter  to  deter- 
mine, a method  has  been  prescribed  that  converts  the  shear  stress  into  local  boundary 
layer  parameters.  In  the  development,  it  was  found  that  Reynolds  number  was  a weak 
function  compared  to  Mach  number  for  the  fluctuating  pressure  field.  Moreover, 
viscous  effects  are  included  in  the  analysis  together  with  compressible  effects,  whereas 
heretofore  investigators  considered  only  the  latter. 

When  assessing  various  fluid  parameters  used  for  non-dimensionalizing  acoustic 
parameters  in  attached  boundary  layer  flows,  it  has  been  found  that:  (1)  the  boundary 
layer  thickness  and  displacement  thickness  have  been  the  primary  characteristic  lengths, 
(2)  local  boundary  layer  edge  velocity  or  a convection  velocity  based  on  some  fraction 
of  the  edge  value  are  characteristic  velocities  used  almost  exclusively,  and  (3)  the  dy- 
namic pressure  has  been  primarily  used  for  normalizing  fluctuating  pressure,  with 
shear  stress  and  local  static  pressure  seldom  employed.  Until  recently,  the  wall  shear 
stress,  the  friction  velocity,  and  the  inner  layer  thickness  u /u^  have  been  neglected 
parameters.  However,  inasmuch  as  the  low  frequency  energy  is  the  dynamic  driving 
function  in  most  aerostructural  applications,  and  the  source  of  the  low  frequency  energy 
is  from  the  outer  portion  of  the  boundary  layers,  those  parameters  associated  with 
this  region  appear  to  be  preferable  for  design  considerations. 

2.1  Further  Comments  Concerning  6*  and  6 

It  has  been  previously  indicated  that  for  incompressible  flow  conditk  s,  charac- 
teristic boundary  layer  thicknesses  6 and  6*  can  be  readily  interchanged.  However, 
for  supersonic  flow  conditions  the  effects  of  compressibility  do  not  allow  for  such  an 
interchange  process;  hence  interpretation  of  prediction  methods  based  on  incompressible 
flow  behavior  must  be  carefully  weighed.  Because  of  the  effect  of  compressibility,  a 
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position  on  a body  where  the  incompressible  boundary  laysr  thickness,  6 , would  equal 

that  for  a compressible  flow  case,  6 would  necessarily  be  a different  position.  This 

c 

is  shown  in  the  following  sketch. 

x > x,  for  4*6, 
c i c i 


For  an  adiabatic,  incompressible,  1/7  velocity  power  law  turbulent  flow,  the  ratio 

6*/6  is  1/8.  On  the  other  hand,  for  a compressible,  adiabatic,  1/7  velocity  power 

law  flow,  the  ratio  6*/6  > 1/8.  According,  the  ratios  of  6 */6  * and  6,/6  should  afford 

i c i c 

some  insight  as  to  the  proper  choice  of  length  when  considering  compressible  effects. 

Consider  then  the  boundary  layer  thickness  developed  for  an  arbitrary  velocity  power 
83 

law  by  Laganelli  et  al 

{ - (n+1)  + [ (”  > ■£*•  + 1]  Cl  + 0.303  M2  ] (43) 

“ naw  e 

together  with  Equation  (20)  for  the  boundary  layer  displacement  thickness. 

For  an  adiabatic  wall  and  1/7  power  law  (n  = 7),  Equations  (20)  and  (43)  become 

2 

6/6  • 10.286  + 0,693  M (44) 

© 


6*/8-  1.  286  + 0.475  M2 


(45) 
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From  boundary  layer  theory,  the  boundary  layer  thickness  can  be  expressed  as 


(>,  *0.371  x Re 
i x 


■1/5 


(46) 
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and  the  compressible  momentum  thickness  as 
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C - 0.0371  x Re  " 1/5  (47) 

T x 

where  eT>  the  compressibility  parameter  defined  by  Equation  (57),  is  unity  for  incom- 
pressible flow  conditions.  Substitution  of  Equation  (47)  into  Equation  (45)  for  an  in- 
compressible flow,  gives 


6 * - (9/7)(.  0371)  x Re  "l/S 
1 x 


(48) 


On  the  other  hand,  substitution  of  Equations  (47)  and  (58)  in'.D  Equation  (45)  gives 


5*  ■ [9/7  + . 0475  M0 


] [ . 0371  x (1  + . 13  M2  )*  64Re 

G X 


(49) 


The  ratio  of  the  above  two  expressions  yields 

9 fi4 

(9/7)  (1  + ,13*rr 
6*/6*  * o 

9/7  + .475  M 

e 


(50) 


One  notes  that  for  M - 0,  the  ratio  8 */b  * - unity.  However,  the  ratio  decreases 
e i c 

rapidly  with  increasing  values  of  M . It  was  assumed  in  the  development  of  Equation 

1/5  e 

(50)  that  the  value  x/Re^  is  invariant  between  the  incompressible  and  compressible 
state,  i.  e. , a fixed  position  on  a model  subject  to  tbe  two  different  flow  conditions. 

In  a similar  manner,  the  ratio  of  boundary  layer  thickness  becomes 


6t/6 


o ~ 


10  (1  + .13  M2  )’64 
0 

10.286  + .693  M2 
e 


(51) 


Equations  (50)  and  (51)  are  plotted  below  to  illustrate  the  functions'  behavior  in  com- 
pressible flow. 
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Here,  one  clearly  sees  the  effect  of  compressibility  on  displacement  thickness  com- 
pared to  boundary  layer  thickness.  The  boundary  layer  thickness  does  not  change  sig- 
nificantly until  hypersonic  conditions  are  approached  (i.e. , M >5).  Consequently 

6 

the  use  of  6 in  prediction  schemes  founded  on  incompressible  flow  principles  does  not 
appear  to  be  a valid  parameter  as  compressible  conditions  occur  over  a surface.  On 
the  other  hand,  the  displacement  thickness  reflects  a significant  change  with  compres- 
sible effects  and  should  be  used  accordingly. 

Based  on  the  preceedlng  developments,  the  characteristic  lengths  and  velocities 
that  are  recommended  for  compressible  flow  are  the  same  as  those  used  for  incom- 
pressible flow.  These  consist  of: 

low  scale  (high  frequency) 
v /u  length 

W T 

u velocity 

T 

high  scale  (low  frequency) 

6*  length 

ug  velocity 

3. 0 MODIFICATION  OF  PREDICTION  METHODS 

This  section  is  concerned  with  development  or  modification  of  prediction  techniques 
that  express  the  power  magnitude,  power  spectral  density  and  cross-spectral  coefficients 
as  functions  of  local  boundary  layer  characteristics.  The  objective  here  is  to  develop 
methods  that  are  logically  constructed  from  fluid  dynamic  and  acoustic  principles  while 
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maintaining  an  easy  to  use  format  for  use  in  design  analyses.  The  rationale  is  pri- 
marily based  on  conoepts  and  data  existing  prior  to  the  experiments  of  this  investi- 
gation. This  is  a consequence  of  the  requirement  to  establish  a baseline  case  of 
attached  turbulent  boundary  layer  flow  behavior. 


3.1  Power  Magnitude 


The  following  methodology  will  consider  the  power  magnitude  for  an  attached  tur- 
bulent boundary  layer  flow  using  the  theoretically  based  concept  of  Lilley  thatp  /t 

rjQ  rms  w 

is  bounded  from  subsonic  to  supersonic  flow  conditions.  Inasmuch  as  the  shear 
stress  is  not  a readily  derivable  term  for  engineering  type  predictions,  other  r ormal- 
ization  terms  will  be  used  as  suggested  by  experimental  data  trends  and  the  phenom- 
enological laws  of  boundary  layer  flow.  The  procedure  wili  be  generalized  for  an 

arbitrary  power  law  turbulent  boundary  layer  and  takes  into  account  wall  temperature 

56 

and  viscous  effects.  The  results  will  be  compared  to  the  works  of  Houbolt  and  Low- 
55 

son,  who  considered  density  variations  only,  and  will  demonstrate  how  to  obtain  the 

constants  which  were  previously  derived  only  from  experimental  data.  It  should  be 

noted  that  a similar  approach  was  developed  in  Reference  67. 

70 

Lilley  found  that  the  normalized  power  magnitude  for  subsonic  conditions  waB 
bounded  in  the  range 


1.7  < p /t 
rms  w 


@ M < 1 

Oft 


whereas,  for  a compressible  flow,  the  ratio  p / t ranged  from  2.2  at  zero  frustum 

rms  w 

Mach  number  to  5. 6 at  a Mach  number  of  10.  The  above  limits  were  experimentally 
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observed  by  Kistler  and  Chen,  Raman  and  Martellucci  et  al.  These  investigators 
noted  that  the  ratio  appeared  to  be  a weak  function  of  the  Reynolds  number  and  Mach 
number.  Based  upon  the  above  premise,  the  ratio  can  be  expressed  as 


PrmgAw  = A,  where  A is  a parameter.  (52) 

Using  the  definition  of  skin  friction  coefficient  together  with  dynamic  pressure, 
the  wall  shear  stress  can  be  expressed  as 


Tw  " <C/2)(PU\  " 2VC/2) 


(53) 


Considering  the  Blasius  form  of  skin  friction  generalizing  to  include  a variable 
power  law  (see  Reference  83),  there  results 


Cj/2  ■ W K (n)  eT  (Res)'2/(3+n) 


(54) 


where  MF  is  the  Mangier  factor  (unity  for  flat  plates),  K(n)  is  a parameter  (.  0296  for 
a 1/7  power  law),  and  the  compressibility  factor  is  defined  as 


- (p*/p J(1+n)/(3+n)(n*/u  )2/(3+n) 


(55) 


In  the  above,  starred  properties  are  based  on  the  classic  Eckert  Reference  Tempera- 
ture method,  namely 


T*/Te  " 1/2  + Tw/Te)  + 0.22  r (^)  M 


(56) 


which,  in  effect,  represents  an  average  through  the  boundary  layer.  However,  it  is 

important  to  note  that  the  compressibility  parameter  includes  both  density  and  viscosity. 

If  one  considers  a constant  mean  pressure  boundary  layer  (dp/dy  = 0)  together  with  the 

4/5 

equation  of  state  and  the  Sutherland  power  viscosity  law  (p,  « T ),  Equation  (55)  be- 
comes 


eT  " CT*/Ter 16/25  - [1/2  (1  + Tw/Te)  + .22r  M2  ]‘‘64 


(57) 


where  a 1/7  (n  - 7)  velocity  power  law  was  assumed.  Figure  17  shows  the  compres- 
sibility factor  as  a function  of  Mach  number  with  wall  temperature  ratio  as  a parameter. 
It  is  quite  interesting  to  note  that  incompressible  data  (M  < 1)  could  reflect  wall  tem- 
perature effects  for  T /T  <1.0. 

w aw 


For  an  adiabatic  flow  (T  =T  ),  Equation  (57)  becomes 

W ftw 


«T  - '1  + 0. 13  M2  )"• 64 
x aw  e ' 


(58) 
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Considering  a flat  plate  geometry  (NF  = unity)  and  a 1/7  power  law.  Equation  (52)  can 
then  be  expressed  as 


P /q 

Frms  e 


2A 


2 A (.  0296) 


f Re  ' [1/2  (1  + T /T  ) + . 22  r ^ 

S W 6 " 


M2  r-64 

e 


(59) 


and  for  adiabatic  flow 


p /q 

rms  e 


2A  (.  0296) 


aw  Re  1/5  [1  + .13  M2  ]*64 
s e 


(60) 


Keeping  in  mind  that  the  objective  here  is  to  lend  some  rationale  to  constants  ob- 
tained from  experimental  data,  we  will  allow  for  the  following  values:  A = 2 which 

6 7 

represents  an  average  of  the  incompressible  results  and  10  < Re  < 10  . Specifically 

6 ^ l/5 

for  an  average  Reynolds  number  condition,  i.e. , Re  ^3.2x10  , Re  ^20.  Hence, 

s s 

Equation  (60)  becomes 

Prms/qe  > " °* 006/(1  + • 13  >* 64  (61) 

aw  e 

A comparison  of  the  above  to  the  Lowson  result  (Equation  (13))  indicates  that  the  co- 
efficient 0. 006  is  a reasonable  choice.  Moreover,  the  essential  difference  in  the  two 
techniques  appears  in  the  exponent  of  the  denominator.  Here,  the  true  effect  of  vis- 
cosity is  noted  inasmuch  as  the  viscous  effect  tends  to  lower  the  value  of  the  exponent 
from  the  Lowson  result.  The  value  chosen  for  the  parameter  A and  Reynolds  number 
(realistic  for  turbulent  flow  requirements),  although  selected  to  acquire  the  value  of 
0. 006,  is  considered  valid  for  comparative  purposes.  This  is  a consequence  of  several 
combinations  of  the  two  parameters  which  remain  in  the  observed  limits  of  P__.0/T,„  with 

rms  w 

corresponding  wind  tunnel  local  Reynolds  number  values. 

Figure  18  shows  Equatlou  (61)  compared  to  predictions  of  Lowson  (adiabatic)  and 
Houbolt  (cold  wall).  It  should  be  noted  that  Houbolt  has  modlfiod  his  cold  wall  prediction 


to  an  adlabatio  result,  which  Is  essentially  the  same  as  Lowson.  Also  shown 

are  data  from  several  experimental  studies  for  both  incompressible  and  compressible 

flows.  While  a large  scatter  is  evident  in  the  experimental  data,  quite  clearly  the 

effects  of  viscosity  and  wall  temperature  appear  to  be  significant  with  increasing  Mach 

number,  a result  that  is  not  surprising.  This  is  believed  to  be  a consequenoe  of  the 

increase  in  p resulting  from  the  viscous  layer  adjacent  to  the  wall, 
rms 

As  an  indication  of  the  importance  of  wall  conditions  Equation  {58}  is  rewritten 
as 


p /q  - 0.006/C1/2CT  /T  • T /T  + 1)  + . 22  r M2  ]'64  (62) 

rms  e ■ ' w aw  w e 1 2 e 

where  for  recovery  factor  r - C.  9 and  y = 1.4  together  with  the  definition  of  the  re- 
covery temperature  (Equation  (6)),  Equation  (62)  becomes 

prms/qe  * °* 006/[1/2  + (VTRW)(i/2  + * 09  + • 04  Me  3* **  <63> 

Sever;il  values  of  T^/T^  were  chosen  and  the  results  are  plotted  in  Figure  19.  A 

comparison  to  the  Lowson  result  indicates  the  importance  oi  viscosity  particularly 

for  the  region  M > 1. 
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3.1.1  Normalizing  with  Static  Pro ssure 

The  previous  arguments,  together  with  experimental  verification,  have  shown  that 
root-mean-square  pressure  fluctuation  when  normalized  by  dynamic  pressure  tends 
to  be  a function  of  Mach  number  and  wall  temperature  under  the  influence  of  viscosity. 
In  order  to  determine  the  magnitude  of  p relative  to  the  local  statio  pressure  for  re 

entry  conditions.  Equation  (63)  was  cast  into  a different  format  using  the  definition  of 

/ 2 

dynamic  pressure,  i.e. , q =y/2P  M , such  that 

6 6 6 


p /p 
rms  e 


t/nkmi 

[1/2  Me"3*  125(1  + Tw/Taw)  + Me_1*  125{.  04 + .09  T^T^  T 64 


(64) 


The  two  limits  that  can  be  developed  from  the  above  equation  arc  as  follows: 

Cold  Wall  (T  /T  « 1) 

w aw ' 


^rms 


cw 


= 0.  0042  [1/2  Me"3* 125  + 0. 04  Me~^ 125]“ °* 64 


for  M > > 1 
e 


/Pc). 


P ' Vr> 
rms  c cw 


18/25 

- 0. 033  Me 
» 1 


Hot  Wall  (T  = T ) - Adiabatic 
w SiW  


(65a) 


(65b) 


P /Pe)  « 0. 0042  125  + 0.13  1L_1* 125 1“°’ 64  (66a) 

rms  ° aw  ° ^ 

for  M >>  1 
e 

P /Pe)flW  - 0.0155  Me18/25  (66b) 

rms  aw 

Me»i 

Here  one  sees  that  the  power  magnitude  will  increase  with  increasing  Mach  number, 
which  now  includes  viscous  effects  through  the  exponent.  Figure  20  shows  rms  pres  - 
sure normalized  by  static  pressure  as  a function  of  Mach  number.  Data  from  several 
experiments  are  shown  compared  to  Equations  (65a)  and  (66a).  Again,  one  must  keep 
in  mind  that  no  attempt  was  made  to  fine-tune  the  analysis  relative  to  a generalized 
velocity  profile,  1.  e. , hilly  developed  attached  turbulent  flow  was  assumed.  Good 
agreement  is  noted,  and  most  importantly,  the  theory  predicts  the  correct  trend.  Finally, 
one  notes  that  the  present  concept  shows  an  order  of  magnitude  difference  from  that 
predicted  by  Lowson  for  M > > 1. 

Also  shown  in  Figure  20  are  data  of  freestream  noise  measurements  as  well  as 
data  obtained  on  the  surface  of  a cone  in  laminar  boundary  layer.  It  is  quite  interesting 
to  note  the  difference  in  level  of  the  rms  pressure  between  the  freestream  and  cone 
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surface  in  the  works  of  Stainback  et  al.  The  authors  noted  that  at  low  shock  strengths, 
the  model  shock  had  little  effect  on  the  ratio  of  rms  sound  pressure  and  local  pressure 
and  that  the  disturbances  behind  the  shock  are  still  predominantly  sound.  Also,  fluc- 
tuating pressure  levels  measured  underneath  the  laminar  portion  of  the  boundary  layer 

differed  significantly  for  the  various  facilities  where  data  were  obtained.  An  exami- 
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nation  of  the  unpublished  measurements  of  Donaldson  as  well  as  Laderman  indi- 
cates a high  level  of  tunnel  noise  when  compared  to  the  cone  data  of  Reference  1 
(M  = 6.  7)  which  were  obtained  in  the  same  facility. 

As  noted  previously,  use  of  the  Crocco  linear  temperature-velocity  relationship 
may  be  questionable  in  turbulent  boundary  layer  type  flows  (Reference  83).  If  the 
Crocco  relation  is  modified  to  a quadratic  distribution,  as  suggested  by  turbulent  pro- 
file data,  the  corresponding  temperature  distribution  at  the  site  of  most  intensive  eddy 
formation  can  be  expressed  as 

(T“  - T )/(T“  - T ) « (u,/u  )2 
1 w'  * e w'  ' 1 e' 

For  adiabatic  wall  conditions,  the  static  temperature  distribution  becomes 

W - w ^ M2  ri  - (u/u  fi 
aw  co 

and  for  u./u  = 1/2,  the  quadratic  relation  gives 
1 6 

T./T  ) - 1 + .135  M2 

l e e 

aw 

which  is  essentially  the  same  as  the  linear  result. 

In  order  to  assess  the  difference  between  the  Lowson  result  (density  change)  and 

the  present  method  (which  combines  both  density  and  viscous  changes),  the  GE  Equili- 

84 

brium  Non-Similar  Boundary  Layer  Program  (ENSBL)  was  used.  This  program 
provides  detailed  boundary  layer  profile  characteristics  using  a finite-difference  solu- 
tion of  the  boundary  layer  equations  and  an  eddy  viscosity  model  for  closure  in  turbulent 
flows.  An  examination  of  profiles  subjected  to  wind  tunnel  conditions  of  Reference  1 
indicate  that  the  choice  of  u./u  is  reasonable  for  the  site  of  eddy  intensity;  however, 

1 G 

the  temperature  distribution  as  provided  by  Lowson  at  this  site  appears  to  be  adequate 
for  the  low  supersonic  range  only  (i.e. , M < 4).  On  the  other  hand,  while 
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under-predicting  the  ENSBL  result,  the  present  method  gave  reasonable  temperature 
predictions  throughout  the  super sonio  and  hypersonic  range. 

3.1.2  Generalization  to  Arbitrary  Velocity  Power  Law 

A generalization  of  the  present  prediction  technique  to  include  both  an  arbitrary 
velocity  power  law  exponent  and  a viscosity  power  law  exponent  will  be  made  in  this 
section.  Equation  (54)  will  retain  its  form  subject  to  Equations  (55)  and  (56).  More- 
over, assumption  of  a constant  pressure  boundary  layer  and  use  of  the  equation  of  state, 
p*/p  = (T  /T*)  is  also  maintained.  The  Sutherland  viscosity  law  is  now  generalized 
to  give85 


uVu  - (T*/T  )m 
e e 


(67) 


where 


T + 198.6 

r_e 3 

In  lT  + 198.6  J 


m - 3/2  + /T  j 

' w e ' 


(68) 


The  compressibility  factor  expressed  by  Equation  (55)  for  an  arbitrary  velocity 
power  law  exponent  when  transformed  into  the  Eckert  reference  temperature  to  include 
the  variable  viscosity  exponent,  becomes 


fT*/T  j£2»»  - (l+n)]/(3+n) 


(69) 


where  TVT  is  given  by  Equation  (56).  The  Mangier  factor,  subject  to  the  arbitrary 
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velocity  exponent,  becomes 


*Fm[2  (2+n)/(l+n)]2/(3+n) 


(70) 


and  Equation  (59)  is  generalized  to  read 


P /q 
rms  e 


2A  K (n)  Mr'  (n) 


Re82/(3+n)  Cl/2  (1  + + .22r  Xy-  M*  ] 


(71) 
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The  following  table  gives  values  of  K(r)  and  hF  (n)  as  functions  of  n. 


VELOCITY  POWER  LAW  PARAMETERS 


n 

7 

8 

9 

10 

K(n) 

.0297 

. 02276 

. 01852 

.0153 

MF(n) 

1.176 

1.156 

1.1404 

1.1275 

K(n)  hF  (n) 

.0332 

. 0263 

.0211 

.0173 

It  should  be  noted  that  hf  is  applicable  to  axisymmetric  shapes  and  has  the  value  of 
unity  for  flat  plates  and  wedges. 

3. 1.3  Sensitivity  Considerations  Concerning  Ve locity  Power  Law  Exponent  -n 

As  previously  noted,  the  power  law  velocity  exponent  has  a value  of  7 for  a fully 
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developed  turbulent  boundary  layer  flow.  It  has  been  demonstrated  ’ that  this 
value  can  be  as  high  as  16  for  flows  corresponding  to  the  end  of  transition  as  defined 
by  surface  beat  transfer.  An  examination  of  Equation  (71)  indicates  that  the  state  of 
the  flow  development  will  not  change  the  rms  pressure  inasmuch  as  the  ratio 

^/(3^a^  invarient  f°r  all  values  of  n. 

Re 

s 

Hence,  variations  in  rms  pressure  will  be  a consequence  of  changes  in  the  parameter 

A,  Mach  number,  or  wall  temperature  ratio. 

In  the  previous  development,  it  was  demonstrated  that  inclusion  of  viscous  effects 

is  the  primary  variation  from  the  work  of  References  55  and  56.  In  particular,  the 

exponent  on  the  compressibility  factor  of  0. 64,  due  to  the  viscosity  contribution,  differed 

from  the  value  of  unity  when  considering  density  changes  only.  The  value  of  m was 
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taken  to  be  0. 8 for  the  method  presented  within.  An  examination  of  flight  data  has 
indicated  that  m has  an  average  value  of  approximately  0. 64  while  wind  tunnel  data 
from  References  1 and  8,  as  well  as  from  the  present  experiments,  has  indicated  that 
the  average  value  of  m is  0. 90.  If  one  considers  a fixed  value  of  n for  the  compres- 
sibility parameters  as  well  as  the  range  of  wall  temperature  and  Mach  number  reported 
in  References  1 and  9,  the  parameter  m was  allowed  to  vary  between  0. 6 and  1. 0.  It 


50. 


was  found  that  the  low  Mach  number  cases  showed  a small  variation  in  the  compres- 
sibility factor  (approximately  10%)  due  to  changes  of  m,  a result  that  had  been  noted 
55  56  67  85 

by  other  analysts.  ’ ' * However,  for  hypersonic  conditions,  variations  in  the 

exponent  m resulted  in  significant  changes  in  the  compressibility  factor  (as  high  as 
50%). 

If  on  the  other  hand,  one  fixes  the  value  of  m and  Mach  number  (for  an  adiabatic 
wall  situation),  a variation  in  the  velocity  power  law  exponent  indicated  only  a 6% 
change  from  M = 4 to  10.  Consequently,  with  the  exception  of  the  choice  of  the  para- 

CO 

meter  A,  the  most  significant  effect  in  use  of  the  compressibility  parameter  appears 
to  be  in  the  coupling  of  density  and  viscosity  as  opposed  to  density  changes  only  (i.e. , 
exponent  0. 64  instead  of  unity).  This  result  appears  to  be  amplified  for  high  Mach 
number  values. 

70 

The  premise  of  the  method  presented  within  is  based  on  the  work  of  Lilley  that 

p /t  is  bounded  over  a small  range  from  subsonic  to  supersonic  conditions.  Others 
rrns  w 

have  noted  that  the  parameter  A,  of  Equation  (52),  was  probably  a weak  function  of 

Reynolds  number  but  should  be  a function  of  Mach  number.  No  one  has  reported  the 

g 

potential  wall  temperature  ratio  effect  on  A.  Raman  correlated  the  p /t  data  as  a 

rms  w 

function  of  momentum  thickness  Reynolds  number  and  noted  that 


p /t 

rms 


w 


a Re, 


(72) 
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for  a and  b constants.  It  is  easily  demonstrated  from  boundary  layer  theory  that 
Re  = Re  (Re  ),  such  that 

U U S 


Reg  ■ constant  eT 


(73) 


where  one  notes  the  compressibility  factor,  e_ , is  a function  of  T /T  and  M . 

1 W RW  0 

When  data  from  Reference  1 were  plotted  according  to  the  above  formulation,  a 

significant  variation  in  slope  was  noted  (the  constant  b in  Raman's  result).  Moreover, 

the  dimensionless  ratio  p /t  was  plotted  as  a function  of  Mach  number  in  Reference 

rms  w 

9 and  showed  a decrease  with  -'creasing  Mach  number,  a trend  that  was  opposite  to 
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the  theory  of  Lilky  as  well  ee  the  tunnel  wall  data  of  Kistler  and  Chen.  However, 


18  37 

data  of  Chyu  and  Hanly  (ogive  cylinder).  Speaker  and  Allman  (wind  tunnel  wall), 
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Willmarth  (tunnel  wall),  and  Chaump  et  al  (sharp  cone),  also  indicated  a down- 
ward trend  with  increasing  Mach  number.  Evaluation  of  the  oarameter  A will  require 
more  details  of  thermodynamic  and  aerodynamic  characteristics  from  reported  experi- 
ments before  attempting  to  assess  a definition  for  more  exacting  prediction  capability. 

When  assessing  the  sensitivity  of  the  velocity  power  law  exponent  (n)  and  the  vis- 
cous power  law  exponent  (m),  it  was  found  that  the  state  of  development  of  the  turbulent 
boundary  layer  had  no  effect  on  terms  that  constitute  the  empirical  constant  0.  006. 

The  most  dramatic  effect  of  the  power  law  exponents  was  founu  to  be  on  the  compres- 
sibility parameter;  in  particular  for  hypersonic  flow  conditions.  The  main  contribution 
was  due  to  viscous  changes  at  these  higher  speeds.  It  is  apparent  that  wall  tempera- 
ture effects  and  subsequent  viscous  changes  in  the  wall  region  of  the  boundary  layer 
must  be  considered  in  prediction  techniques  for  hypersonic  flow  conditions. 

It  is  recommended  that  Equation  (63),  which  was  subsequently  generalized  for  ar- 
bitrary velocity  and  viscosity  power  laws  in  Equation  (71),  be  used  to  predict  the  rms 
fluctuating  pressure.  The  normalizing  parameter  considered  for  the  prediction  tech- 
nique consists  of  the  local  dynamic  pressure.  In  order  to  determine  the  relative 
magnitude  of  the  pressure  fluctuating  field  compared  to  the  mean  static  value,  Equation 
(€4)  can  be  used. 

3.2  Power  Spectral  Density 
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This  section  discusses  the  concepts  of  Black  and  Houbolt  who  considered  a 
first  moment  type  representation  of  the  power  spectral  density  (PSD).  Black  had  noted 
that  the  strength  of  a pressure  signature  beneath  a vortex  system  will  generally  be 
determined  by  local  wall  shear  stress,  t . Moreover,  for  a particular  class  of  vortex 
systems  of  mean  spacing  1 , and  wave  velocity  u , the  mean  frequency  of  passage  for 
this  system  over  the  wall  will  be  u^/l^  As  a result,  the  contribution  of  the  vortex 
system  to  the  fluctuating  wall  spectrum  should  scale  with  and  Uj/1^,  such  that 


d prms  /Tw 
«i 


.wlx 


- f (=3*  ) 

Ui 


(74) 
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For  small  vortex  systems,  where  heights  are  of  the  order  of  the  viscous  sub- 
layer thickness,  length  and  velocity  scales  are  governed  by  near  wall  conditions  and 
Equation  (74)  is  represented  as 


2 , 2 
dp  /t 

„rms_w  _ .omav. 

d(^W)  l(^ 

i»2 


(75) 
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On  the  other  hand,  large  voitex  systems  will  scale  to  boundary  layer  edge  values  and 
Equation  (74)  can  be  written  as 


. 2 / 2 
dp  /t 
*rms  w 

id  6 . 

d(“ ) 
Utw 


- F, 


(»i., 

* ii 


'w 


(76) 


It  should  be  noted  that  most  investigators  have  chosen  6*  and  u as  the  representative 

G 

scaling  parameters  for  the  region  associated  with  Equation  (76). 

Black  found  it  more  convenient  to  examine  the  first  moment  of  the  PSD  rather  than 
the  power  spectrum  per  se.  In  terms  of  first  moment  representation,  the  functional 
forms  of  Equations  (75)  and  (76)  become 


and 
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Black  further  noted  that  dimensionless  frequenciei  represented  above  are  related  by 
the  parameter  Su-^/kw,  such  that  the  overall  rms  value  of  the  pressure  field  should 
have  tho  functional  form 


<p2  )1/2  UTul 

- F.  ( Tw°  ) 


'w 


'Vv 


(79) 
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Equation  (79)  again  suggests  a possible  means  of  evaluating  the  parameter  A as  deter- 
mined by  Equation  (52).  However.  If  the  boundary  scaling  parameters  6*  and  u aro 
used  for  the  far  field  region  (y  > 6^),  the  right  hand  side  of  Equation  (79)  becomes 

6*UTW/Uwue' 

Houbolt  considered  the  empirical  representation  of  ike  PSD  as  a function  of 

Strouhal  number,  Equation  (24),  such  that  the  density  distribution  of  a plot  of  p2  as 

rms 

given  by  tvcp(u>)  would  maximize  when  the  Strouhal  number  is  unity.  TUlr  was  accom- 
plished by  a change  in  variable,  Equation  (25),  to  the  definition  of  the  rms  pressure. 
The  rationale  used  by  Houbolt  is  as  follows.  Suppose  one  considers  Equation  (24) 

<P  M/Prms  uc(iu)/6*  - 1/[1  + MVu,)2] 

One  can  take  any  function,  say  F,  of  an  arbitrary  variable  x to  give 

F - 1/(1  + x2) 

analogous  to  the  above  empirical  formulation  of  the  PSD.  If  one  plots  the  function  F 
versus  x on  a linear  scale,  the  ordinate  has  a value  of  unity  (x  = 0)  and  becomes 
asymptotic  to  the  abscissa.  If  one  multiplies  both  sides  of  the  equation  by  x and  plots 
xF  versus  x on  a linear  scale,  the  ordinate  has  a maximum  value  of  1/2  about  x ~ unity 
and  becomes  slowly  asymptotic  to  the  abscissa  once  again.  Houbolt  then  noted  that  if 
the  product  xF  (linear  scale)  were  plotted  as  a function  of  the  Inr  (leg  scale),  that  the 
xF  product  would  maximize  at  1/2  when  x - unity,  and  rapidly  decay  to  asymptotic 
values  of  the  abscissa  for  In  x < 3 and  In  x > 1.  This  Is  shown  helow  schematically 


Houbolt  considered  the  change  of  variable  represented  by  Equation  (25)  together 
with  the  form  of  the  PSD  correlation,  Equation  (24)  and  suggested  the  functional 
relationship  represented  by  Equation  (26)  with  the  moment  product  oucp(uj).  Although 
Houbolt  did  not  compare  any  measurements  v/ith  the  above  format,  Black  compared 
his  method  to  data  of  several  experiments,  both  subsonic  and  supersonic.  Inspection 
indicated  that  tht  first  moment  type  representation  of  data  and  Black’s  method  display 
the  characteristics  suggested  by  Houbolt  and  by  the  sketch  provided  above. 

The  following  is  the  analogy  adapted  in  the  present  investigation  regarding  first 
moment  type  representation  of  the  PSD.  Consider  a typical  plot  of  the  PRD  spectrum, 
shown  below  schematically 


By  definition,  the  rms  pressure  (0^  moment)  is  expressed  its 

p2  ■ f"  cp(tu)d(D"  tp(uj)du>  + f cp(u>)dw 
Frms  Jo  o w 

Much  of  the  data  considered  throughout  the  literature  is  analyzed  out  to  approximately 
20  KHz.  Hence,  If  one  considers  the  second  integral  on  the  right  hand  side  above 
negligible,  the  first  moment  can  be  approximated 

(81) 


prm«L  ~ 

If.  m. 


and  differentiation  gives 


dp  /dm  ■ oi cc 

rms 


(82) 
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One  immediately  reeognizes  the  similarity  of  the  above  to  both  Black's  and 
Houbolt's  interpretation.  From  Equation  (82),  the  functional  form 

2 1 

uu  <p  (ul  )/ p * function  (—)  (83) 

rms  v 

should  be  representative  of  the  PSD  distribution.  Here  the  characteristic  length  1, 

and  velocity  v take  the  required  near  field  (wall)  values  v and  uT  and  far  field 

w w 

(y  > 6,)  parameters  6*  and  u , respectively.  Moreover ,p2  is  the  rms  value  obtained 
1 e rms 

to  the  frequency  cut-off  point.  One  notes  that  the  above  is  actually  a weighted  resi- 
dual type  fit  to  data  where  the  spectrum  is  weighted  in  the  frequency  range  where  data 
were  obtained. 

As  an  indication  of  the  first  moment  type  representation  for  displaying  data,  Fig- 
ures 21  and  22  were  prepared.  Data  were  taken  from  References  1,  9,  14,  and  18 
for  attached  flow  conditions.  Moreover,  the  data  represent  fluctuating  pressure  en- 
vironments over  cones,  plates,  ogive-cylinders  and  wind  tunnel  walls.  Figure  2 
shows  the  PSD  as  a function  of  Strcuhal  number  for  high  scale  coordinates.  Spectra 
have  been  normalized  using  conventional  techniques  in  the  literature.  It  is  quite  ap- 
parent that  there  is  a significant  scatter  in  data  displayed  in  this  fashion;  in  particular 

the  ogive-cylinder  in  the  low  frequency  portion  and  the  M =4  sharp  cone  measure- 

00 

ments  mid-range  portion  of  the  spectrum.  The  latter  is  believed  to  be  a consequence 
of  transitional  flow  behavior. 

The  same  data  were  replotted  in  first  moment  type  format  in  Figures  21  and  22 
for  the  far  field  and  near  field  (wall)  forms  of  the  Strouhal  number,  respectively.  In 
Figure  21  one  notes  that  the  spectra  tend  to  coalesce  rather  significantly  within  each 
data  group.  On  an  overall  basis,  the  data  appear  to  be  more  uniformly  represented 
in  the  coordinates  shown.  From  these  representations,  no  discernable  Mach  number 
or  Reynolds  number  trend  is  noted. 

Figure  22  shows  the  first  moment  type  representation  with  the  high  frequency 
scaling  of  the  Strouhal  number.  Again,  a coalescense  of  data  is  noted  as  well  as  a 
definite  pattern.  Here,  a Mach  number  effect  is  apparent  for  the  data  group  represen- 
ted by  Reference  1.  These  data  were  taken  because  of  the  difficulty  in  obtaining  the 
wall  conditions,  and  uT# , for  the  experiments  of  References  9,  14,  and  18. 
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It  is  also  quite  interesting  to  note  that  data  appearing  in  first  moment  type  format 
displayed  the  characteristics  suggested  by  Houbolt;  for  both  the  near  and  far  field 
Strouhal  numbers.  For  example  the  ordinate  appears  to  peak  at  approximately  one- 
half  for  a Strouhal  number  approaching  unity.  It  is  noted  that  if  the  circular  frequency 
were  used,  the  abscissa  would  shift  the  data  to  the  right  by  a factor  of  2rr. 

3.2.1  Further  Comments  on  Power  Spectral  Density 
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Houboltd ' examined  the  data  of  several  experiments  * ' and  noted  that  the 

spectra  could  be  represented  by  the  empirical  formula 
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in  the  range  of  0. 2 < u>  6*/u  < 20  which  is  schematically  represented  below 
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Inasmuch  as  the  data  indicate  a peak  in  the  low  frequency  range  (point  A),  the  empirical 
representation  was  considered  to  be  flat  in  this  range.  This  situation  is  considered  a 
conservative  approach  in  structural  applications  where  low  frequency  responses  (f  < 
1000  Hz)  tend  to  predominate.  From  the  definition  of  the  mean  square  pressure  and 
utilizing  Equation  (84),  there  results 

°2  " Prms  " A 2 x 10-5  f q* 
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Here,  we  note  that  the  constant  tt/2  is  a consequence  of  the  functional  form  of  the 
empirical  formula  representing  the  spectrum.  Normalized  rms  pressure  is  deter- 
mined to  be  o/q  = 0. 0056  whioh  is  approximately  the  value  of  measured  incom 
6 

pressible  flow  data. 

On  the  basis  of  the  above  development,  it  appears  that  a reasonable  approach  to 

evaluate  empirical  formulae  for  power  spectral  density,  i.e. , Equations  (24),  (28), 

and  (29),  is  to  assess  the  zero  frequency  intercept  value  and  the  normalized  rms 

ores  sure  values  resulting  from  integration  of  the  various  equations.  Moreover,  the 
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incompressible  data  of  Bull  and  Blake0  will  be  used  as  a baseline  for  testing  the 
techniques.  It  should  be  noted  that  these  data  have  been  well  documented  and  reviewed 

in  the  scientific  community.  Moreover,  the  recent  assessment  of  these  data  by  Will- 
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marth  and  counterargued  by  Bull  and  Thomas  indicate  the  credibility  of  the  Bull 

data  for  a baseline  case. 

To  integrate  the  various  empirical  power  spectral  density  equations,  consider 
the  following  definite  integral 


a,  r-1 

f ~.dx 

°(p+qxs)n+1 
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for  the  condition  0 < r/s  < n + 1.  In  the  above  F (z)  is  the  gamma  function  = r°e-ttz_1dt. 

Jo 

The  three  equations  to  be  evaluated  are  Houbolt,  Equation  (24),  Lowson,  Equation  (28), 
and  Robertson,  Equation  (29)  and  will  be  designated  as  H,  L and  R,  respectively.  Con- 
sidering incompressible  flow  and  the  format  of  Equation  (84),  the  three  equations  become 
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where  $ (0)  is  the  normalized  zero  intercept  (low  frequency)  value  of  the  power  spectra. 
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Figure  23  shows  the  data  of  Bull,  Blake  and  Schloemer.  It  appears  that  the 
intercept  value  of  the  Bull  data  is  approximately  2.2  x 10  while  that  of  Blake  approxi- 
mate ly  1.3  x ID'5 

If  we  consider  the  definition  of  the  rms  pressure  and  integrate  each  of  the  above 
expressions,  there  results 

(c/qJ2H  ■ fv(oy,  (c/qj*  - <P  (o);  (c/qm)^  - J$(o) 


which  respectively  yields 
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Hence,  it  appears  that  the  Houbolt  functional  form  best  represents  the  measured  incom- 
pressible normalized  rms  pressure  data.  Relative  to  the  zero  intercept  value,  if  we 
allow  o/q  to  take  on  the  value  0. 006  and  consider  the  following 
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Thus,  considering  a limiting  representation  (M  - 0)  of  each  empirical  format,  it  appears 
again  that  the  Houbolt  concept  best  represents  the  data  of  Bull.  One  should  note  that 
o/q  - f [M,Re,T^]  and  approaches  a value  of  approximately  0.006  for  incompressible 
flow  conditions. 
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If  the  Houbolt  algorithm  Is  reconsidered  over  the  entire  spectrum 
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where  K = k 1/V,  for  1 and  V some  arbitrary  length  and  velocity  and  k a parameter.  A- 
gain  using  the  definition  of  the  power  magnitude  and  Equation  (86)  there  results 
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and  Equation  (86)  becomes 
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where  k = k/2rr.  For  zero  intercept  values,  the  above  is  written  as 
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Now  Bull  reported  values  of  o/q  = 0. 005  while  Blake  reported  values  of  o/q  * 0. 00876; 
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the  latter  result  is  considered  quite  high.  If  we  allow  for  k = unity,  the  normalized 
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zero  intercept  values  of  PSD  become  1.  59  x 10  for  Bull  and  4. 89  x 10  for  Blake. 

One  sees  from  Figure  23  that  these  values  are  not  in  concert  with  the  data. 

In  order  to  match  the  data  as  u>  - 0,  E would  have  values  of  1. 3825  and  0. 266  for 
Bull  and  Blake,  respectively.  If  these  values  were  used  in  Equation  (87)  with  the 

measured  values  of  a/q  , the  equation  would  underprediet  the  Bull  data  while  signifl- 
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cantly  overpredicting  the  Blake  data.  If  we  allowed  o/q^  = 0. 006  and  k =*  unity  (Houbolt 
result),  the  zero  intercept  value  by  Equation  (88),  would  be  2.29  x 10  a result  thar 
matches  the  Bull  data  quite  nicely.  Consequently,  Equation  (87)  with  1c  « unity  is  con- 
sidered the  most  appropriate  algorithm  for  predicting  incompressible  power  spectral 
density. 

Keeping  in  mind  that  Houbolt's  functional  form  does  not  preclude  compressible 
conditions,  the  following  is  offered  as  an  appropriate  generalization.  Using  Bull's  data 


aa  a baseline  as  well  as  the  Houbolt  functional  form.  Equation  (87)  can  be  written  as 


cpCwJue 
q%  6* 


, 2 A~ 

(a/qaJj  n k 
1+  (k)2(w6*/ue)2 


k ^ unity 


Raman  determined  a similar  expression  for  compressible  flow  that  considered  the 

functional  form  of  Lowson,  namely  Equation  (32).  It  was  previously  stated  that  the 

. 2 

incorporation  of  the  (a/q  ) hi  the  above  acts  like  a transformation  function  for  com- 
pressible  effects  as  given  by  Equation  (53).  Thus  for  <u  - 0,  we  can  write 


[cp(o)ue/qe  5*  Jc 

4 


(89) 


We  recognize  that  the  right  hand  side  of  the  above  is  constant,  as  well  as  ( a/q  ) which 

0 1 

has  a value  of  0. 006.  Inasmuch  as  the  compressibility  factor,  e,p,  is  an  inverse  func- 
tion of  Mach  number  (Figure  37),  then  the  zero  intercept  value  for  compressible  flow 

should  decrease  with  increasing  M . Moreover,  use  of  the  incompressible  formulation 

e 2 2 

for  PSD  via  the  Houbolt  format  that  considers  q instead  of  a for  normalization  should 

0 

allow  for  compressible  data  to  be  transformed  into  the  incompressible  plane  using  the 
compressibility  factor.  Figure  24  shows  the  results  from  several  compressible  ex- 
periments as  well  as  the  incompressible  data  of  Bull  and  Blake.  The  compressible 
data  Indicate  that  the  zero  intercept  valuec  ere  indeed  less  than  the  incompressible 

_5 

value  (2.29  x 10  ).  Hence,  Equation  (89)  will  serve  as  the  basis  for  which  compres- 
sible data  will  be  evaluated  in  the  current  investigation. 

Some  comments  are  in  oraer  relative  to  gage  size  and  measuring  errors  in  the 
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high  frequency  range.  Wlllmarth  noted  that  earlier  investigations  (prior  to  Blake) 

used  large  transducers  that  could  not  resolve  small  scale  fluctuations  even  when  using 

gage  size  corrections.  Figure  23  shows  such  a disparity  between  the  data  of  Bull  and 
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Blake  for  Strouhal  numbers  > 1.  Bull  and  Thomas  performed  an  experiment  to  deter- 
mine the  difference  between  pinhole  microphone  measurements  (such  as  Blake)  and 
piezoelectric  transducers  as  used  previously.  It  was  found  that  the  pinhole  caused 
spurious  contributions  to  cp(<u),  up  to  factors  of  4.  Consequently,  assuming  some  error 
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in  the  Bull  data  such  that  the  integrated  area  of  the  spectra  would  yield  slightly  higher 

values  of  a/q  than  the  measured  values  (.  005).  we  will  assume  that  the  normalized 

power  magnitude  has  a value  of  0. 006  for  incompressible  flow.  This  appears  to  be  a 

reasonable  choice  due  to  the  possible  high  frequency  errors  from  both  experiments. 

Hence,  use  of  the  above  value  for  a/q  and  k = unity  yields  the  desired  result  for  the 

00 

zero  intercept  value.  When  these  value&are  used  in  Equation  (87),  the  Houbolt  format 

fits  the  Bull  data  over  his  predicted  range  0. 2 < m 6*/u  < 20. 

c 

It  is  observed  that  a significant  dispersion  of  the  data  occurs  for  Strouhal  numbers 
-2 

< 5 x 10  . Due  to  the  high  velocities  attained  at  hypersonic  conditions  the  compressible 

Strouhal  number  can  be  less  than  the  incompressible  values  by  over  an  order  of  mag- 

15 

nitude.  This  is  reflected  in  the  compressible  data.  Moreover,  Dods  and  Hanly 
noted  that  at  low  frequencies  spurious  contributions  are  accorded  to  the  pressure 
fluctuating  field  from  tunnel  noise.  We  have  already  noted,  Figure  20,  that  the  AEDC 

14 

von  Karman  A and  B Facilities  are  considered  high  noise  level  tunnels.  Also  Coe  et  al 
recognized  that  statistical  accuracy  is  reduced  in  the  low  frequency  range.  Wlllmarth 

in  his  review  commented  that  many  experimenters  terminate  the  spectra  for  w6*/u  < 

-1  00 
10  because  measurements  are  obscured  due  to  tunnel  noise.  Other  factors  leading 

to  spurious  results  are  attributed  to  exact  definitions  of  the  transition  zones,  vibration 

modes,  surface/gage  roughness,  pressure  gradient,  gage  size  (high  frequency  range) 

and  the  interaction  of  turbulence  with  turbulence. 


62. 


SECTION  IV 


EXPERIMENTAL  PROGRAM 


1.0  INTRODUCTION 

This  section  presents  a brief  review  of  the  ground  test  program,  described  in  de- 
tail in  Volume  U.  Wind  tunnel  experiments  were  conducted  at  Arnold  Engineering 
Development  Center  (AEDC)  in  the  von  Karman  Gas  Dynamics  Facilities  Tunnels  A 
and  B at  Mach  numbers  4 and  8,  respectively.  The  objective  of  the  test  program  was 
to  obtain  acoustic  data  for  the  purpose  of  verifying  or  redefining  aeroacoustic  pre- 
diction techniques  currently  used  in  the  aerospace  community.  In  particular,  tests 
were  made  on  a maneuvering  configuration  that  featured  a control  surface.  Tunnel 
operating  conditions  and  acoustic  sensor  (Gulton  and  Kul lie  gages)  locations  were 
chosen  such  that  measurements  were  made  in  regions  of  laminar,  transitional  and  tur- 
bulent boundary  layer  flow.  Variations  resulting  from  the  effects  of  Mach  number, 
Reynolds  number,  angle  of  attack,  yaw,  nose  bluntness  and  control  surface  flap  angles 
were  examined.  Several  acoustic  gages  were  located  in  regions  where  laminar  bound- 
ary layer  conditions  prevailed  for  the  purpose  of  measuring  tunnel  tare  noise  levels. 

Figure  26  shows  a schematic  of  the  model  with  surface  sensor  locations.  Stations 
3.6  and  8. 35  represent  the  flap  and  slice  regions,  respectively.  Further  details  con- 
cerning the  instrumentation  are  given  in  Volume  Q.  Table  I shows  a test  summary  of 
tunnel  conditions,  model  attitude  and  model  geometric  variations  for  which  data  were 
obtained.  Because  of  the  potential  volume  of  data  that  could  be  generated  from  all  the 
tests  run,  a data  sub-set  was  established  for  analysis  hi  this  investigation.  Table  n 
shows  the  conditions  for  this  sub-set.  Also  Included  in  this  table  are  the  transition 

zones  (onset  and  end)  as  determined  from  heat  transfer  data.  Tables  III  and  IV  present 
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nominal  local  flow  properties  generated  from  the  GE-RE3D  3DFF  and  3DV  pro- 
grams for  the  data  sub-set  tunnel  and  model  geometric  conditions.  It  is  noted  that 
properties  shown  are  those  Identified  as  important  normalization  parameters  for  aero- 
acoustic phenomenology.  Tables  V and  VI  present  the  calculated  displacement  thickness 
values  of  the  data  sub-set  for  main  ray  sensor  locations.  It  should  be  noted  that  Tables 
HI  and  IV  consider  nominal  values  of  the  properties  along  the  Important  sections  of  the 
model  (l.e. , conical,  slice  and  flap  regions)  where  acoustic  sensor  arrays  exist.  While 
it  is  recognized  that  some  variation  does  exist,  reported  values  can  be  used  for  quantitative 
results. 


63 


2.  0 SUMMARY  OF  TEST  DATA 


2. 1 Thermal  and  Acoustic  Distributions  - Mach  4 

Figure  26  shows  a typical  result  where  the  extent  of  transition  is  determined  from 
a plot  of  the  Stanton  number  as  a function  of  social  distance.  Transition  onset  is  defined 
as  the  point  where  the  local  heat  transfer  is  a minimum  and  the  end  of  transition  is  de- 
fined as  the  point  where  the  local  heat  transfer  is  & maximum.  Also  noted  in  this  figure 
is  the  decrease  in  heat  flux  over  the  slice  areas  which  results  from  flow  expansion 
from  the  conical  to  the  flat  (parallel  to  freestream)  surface. 

Figures  27  through  32  summarize  the  overall  measured  sound  pressure  level  (dB, 
20  Hz -20  KHz)  as  well  as  the  corresponding  heat  transfer  distribution  along  the  model 
surface  for  several  of  the  Mach  4 test  conditions.  Values  of  overall  sound  pressure 
level  (OASPL)  shown  therein  have  been  adjusted  from  raw  sound  pressure  level  data 
to  account  for  obviously  spurious  contributions  to  the  total  acoustic  signal  resulting 
from  vibration-induced  sensor  output  at  low  frequency  (i.e. , < 1000  Hz).  This  pro- 
cedure reduced  scatter  in  most  data,  but  high  frequency  contributions,  which  appeared 
to  affect  sensor  A- 9 (also  A-3  in  Run  10),  were  not  explicitly  eliminated  in  Figures  27 
to  32.  The  transition  zones  shown  in  these  figures  were  inferred  from  experimental 
data  as  well  as  analysis.  This  was  required  to  insure  that  measured  heat  flux  levels 
were  consistent  with  the  assumed  state  of  tne  boundary  layer  flow.  In  particular,  the 
complexity  of  flow  ever  the  control  surface  region  yielded  heat  flux  levels  that  were 
difficult  to  interpret  relative  to  transition  effects.  For  all  test  conditions  examined, 
the  flap  region  experienced  turbulent  flow  behavior.  Finally,  it  should  be  noted  that 
when  values  of  acoustic  data  from  run  to  run  at  identical  angle  of  attack  were  within 
experimental  scatter  (i.e. , approximately  ± IdE),  average  values  of  OASPL  are  shown, 
in  order  to  simplify  graphical  representation. 

Data  in  Figures  27  through  29  illustrate  the  effect  of  increased  flap  deflection. 

While  both  the  acoustic  and  thermal  levels  over  the  flap  increased  with  flap  deflection, 
there  does  not  appear  to  be  any  upstream  effects  along  the  cone  and  slice  regions.  This 
condition  prevailed  at  both  0*  and  7*  angles  of  attack  (at  Re^/ft.  « 3.6  x 106)  for  flap 
angles  ranging  from  0*  to  20*  in  spite  of  increases  in  flop  acoustic  levels  of  up  to  17  dB 
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for  the  higher  flap  angles.  At  Res  /ft.  =2.5x10  , data  along  the  conical  portion  of 

CD 

the  model  exhibited  similar  behavior;  however,  for  the  lower  Reynolds  number  a 
variation  in  acoustic  and  heat  transfer  level  over  the  slice  is  noted  with  increasing  flap 
angle  (Figure  29).  It  should  be  noted  that  this  region  is  in  a transitional  flow  state 
such  that  increased  flap  deflections  could  accentuate  the  flow  instability  in  the  area  of 
the  joint  between  the  slice  and  flap. 

Figures  27  and  28  exhibit  a coincidence  of  peak  acoustic  level  with  transition  zone 
along  the  conical  portion  of  the  model.  This  maximum  level  occurs  near  the  end  of 
transition  where  the  heat  transfer  is  also  maximum.  A clear  indication  of  flow  expan- 
sion from  the  conical  part  of  the  model  to  the  slice  is  noted  in  lower  levels  of  both 
acoustic  and  thermal  data.  This  is  a consequence  of  the  pressure  drop  as  the  flow 
turns  through  the  expansion  wave  at  the  intersection  of  the  cone-slice  region.  More- 
over, non-center  line  acoustic  gages  (located  on  the  cone)  show  a higher  level  due  to 
the  local  pressure  of  the  conical  region.  As  might  be  expected,  increased  levels  in 
the  acoustic  and  thermal  data  with  increasing  flap  deflections  are  a consequence  of  the 
higher  local  pressure  associated  with  these  deflections.  The  same  characteristic  can 
be  noted  on  the  conical  portion  of  the  model  when  comparing  the  0°  and  7°  angle  of 
attack  data. 

Figure  30  illustrates  the  Reynolds  number  effect  for  identical  model  conditions 
( a - Rj^j/Rg  = \a=  0 for  6F  = -15® ).  It  is  evident  that  the  transition  rone  over  the  cone 
moves  forward  for  increasing  values  of  Re  . As  an  example,  acoustic  sensor  A-19 

OD 

is  measuring  levels  near  the  end  of  transition  for  Run  69,  while  the  peak  value  has 
moved  forward  to  approximately  the  location  of  A-10  in  Run  27.  On  the  slice,  although 
there  Is  no  clear  trend  in  acoustic  data,  thermal  measurements  indicate  that  heat  trans- 
fer levels  associated  with  turbulent  flow  in  Run  27  are  higher  than  levels  of  early  transi  - 
tion exhibited  by  Run  69.  On  the  flap,  this  relationship  is  reversed.  Inasmuch  as  the 

-1/2 

heat  transfer  coefficient  is  proportional  to  Re^  , the  lower  Reynolds  number  case 
will  yield  higher  Stanton  numbers  along  the  oonlcal  portion  of  the  model.  On  the  other 
hand,  lower  values  of  Reynolds  number  delay  transition;  hence,  the  reason  for  lower 
heat  flux  values  along  the  slioe  (Figure  29  shows  the  slice  in  a transitional  mode). 

Since  the  flap  region  is  in  a turbulent  flow,  the  lower  Reynolds  number  case  again 

-1/5 

yields  higher  Stanton  numbers  (here,  the  heat  transfer  coefficient  is  proportional  to  Re^  ). 
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As  indicated  by  the  acoustic  and  thermal,  data  of  Figure  31,  the  effect  of  Increased 
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angle  of  attack  (R.  ./R_  = \ - 6_  = 0 for  Re  /ft.  * 3. 6 x 10  ) is  to  drive  maximum 
transitional  levels  forward  on  the  cone.  Increasing  peak  values  of  OASPL  as  angle 
of  attack  increases  (i.  o. , A- 10/9/3  for  Runs  6,9,10,  respectively)  primarily  result 
from  high  frequency  contributions  to  the  sound  pressure  level.  Turbulent  levels  of 
fluctuating  pressure  and  heat  transfer  on  the  slice  and  on  the  flap  also  increase  with 
increasing  angle  of  attack. 

Figure  32  demonstrates  the  effect  of  nose  bluntness  including  mass  transfer  at 
Mach  4 (i.e. , injection  of  mass  into  the  boundary  layer  at  the  nosetip).  Shown  for 
comparison  is  a fairing  of  the  sharp  nose  data  of  Run  42.  The  reason  for  injecting 
mass  at  the  nosetip  is  twofold.  First  it  destabilizes  the  boundary  layer  thereby  ren- 
dering an  effective  boundary  layer  control  mechanism.  Second,  it  simulates  an  actual 
ablation  process  experienced  on  a re-entry  vehicle.  In  Run  102,  it  is  evident  that  the 
relatively  constant  levels  on  the  cone  (^  124  dB)  are  indicative  of  laminar  flow  (tunnel 
nois^S  while  sensors  A-19,  20  and  ?1  are  monitoring  higi  levels  near  the  end  of 
transition.  The  slice  is  in  transitional  flow,  a fact  borne  out  by  the  monotonically  in- 
creasing acoustic  levels  from  sensors  A-13,  14,  15  and  17  as  well  as  from  theoretical 
predictions  of  the  heat  transfer  distributions  (when  compared  to  measured  data). 

The  laminar  sharp  cone  heat  transfer  data  fairing  (of  Run  42)  agrees  very  well 
with  th  at  of  the  blunt  nose  data  indicating  that  the  entropy  layer  has  recovered  to  the 
sharp  co’  i case  prior  to  the  first  recording  heat  flux  sensor.  Moreover,  the  blowing 
case  exhibits  heat  flux  levels  identical  with  the  sharp  cone  values  at  the  end  of  transi- 
tion (conical)  as  well  as  along  the  slice  and  flap.  The  acoustic  data  present  a different 
picture.  Sharp  cone  values  show  a significantly  higher  dB  level  than  the  blunt  cases 
along  the  mid-cone  region,  which  ir  a consequence  of  transition.  The  slice  characteristics 
are  similar  while  the  sharp  cone  values  appear  much  lower  along  the  flap  than  the  blunt 
values.  Of  particular  interest  along  the  slice  is  the  apparent  lower  acoustic  levels 
than  tunnel  noise  values  (as  measured  on  the  conical  frustum).  It  Is  believed  that  this 
phenomenon  is  a result  of  deflecting  the  noise  pressure  waves  from  the  tunnel  wall  as 
they  pass  through  the  bow  shock  and  expansion  wave  (at  the  intersection  of  cone-slioe). 
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The  latter  eauses  the  pressure  waves  to  accelerate  toward  the  freestream  direction  as 
a consequence  of  increasing  velocity  gradient  through  the  expansion  wave. 

2,2  Thermal  and  Acoustic  Distributions  - Mach  8 

Figure  33  presents  data  corresponding  to  the  basic  zero  angle  of  attack,  sharp 

cone  configuration  for  values  of  flap  deflection  between  0*  and  -20°  at  Re  = 3. 7 x 

6 w 
10  /ft.  It  is  evident  from  a comparison  with  Figure  27  that  the  transition  zone  has 

moved  forward  on  the  model  for  the  higher  Mach  number  condition.  Also,  laminar 

acoustic  levels  (A-2  in  Figure  33)  are  higher  in  Tunnel  B than  those  monitored  in 

Tunnel  A.  Both  these  trends  were  noted  in  the  1972  experiments  conducted  at  AEDC 

discussed  in  Reference  1. 

Fluctuating  pressure  levels  measured  in  turbulent  flow  on  the  slice  and  presented 
in  Figure  33  are  approximately  5dB  lower  at  Mach  8 than  corresponding  OASPL  data 
at  Mach  4.  This  relationship  in  the  acoustic  data  is  consistent  with  slice  values  of 
heat  flux,  respectively  shown  for  Mach  4 and  Mach  8 in  Figures  27  and  33.  Flap  de- 
flection has  virtually  no  effect  on  slice  acoustic  and  thermal  environments  for  turbulent 
flow  conditions,  a situation  also  noted  in  the  Mach  4 data  of  Figures  27  and  28.  However, 
fluctuating  pressure  loading  on  the  flap  itself  at  Mach  6 is  more  affected  by  flap  deflec- 
tion than  is  loading  at  Mach  4.  This  is  apparent  in  Figures  27  and  33,  where  flap  OASPL 
values  at  Mach  4 and  8 are  seen  to  increase  by  less  than  20  dB,  and  by  more  than  25  dB, 
respectively,  for  8 between  0*  and  -20* . 

Acoustic  and  thermal  measurements  at  Re  = 2. 5 x 10  /ft. , shown  in  Figure  34, 
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indicate  that  the  transition  zone  has  moved  aft  on  the  cone  from  its  location  at  the  higher 
Reynolds  number  condition.  The  same  trend  was  noted  in  Mach  4 data  of  Figure  30. 

Peak  acoustic  levels  in  transition  apparently  were  not  measured  in  these  runs,  with 
the  end  of  transition  occurring  between  frustum  sensors  A-10  and  A-19  through  A -22. 
Examination  of  heat  transfer  data  in  this  region  suggests  that  acoustic  gages  A-9  and 
A-10  are  registering  respectively  too  high  and  too  low.  Increasing  the  flap  deflection 
from  -15*  to  -20*  has  no  significant  effect  on  cor  leal  frustum  or  slice  aooustic/tbermal 
measurements,  for  example  at  6p  ■ 0* , heat  transfer  levels  are  virtually  coincident 
with  other  data  (6p  / 0*)  everywhere  except  on  the  flap  itself.  At  the  lower  Reynolds 
number  (2.5  x 10*/ft.),  values  of  flap  OASPL  at  « -20*  are  slightly  below  those 
measured  at  Rea  * 3. 7 x 106/ft. , however,  noting  the  flap  data  of  Figure  29  it  Is  evident 
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that  the  Maoh  8 levels  are  significantly  less  than  Mach  4 acoustic  data  at  He  =2.5x 
106/ft. 

Effects  of  model  angle  of  attack  at  Mach  8 and  x = 0°  are  indicated  in  Figure  35. 

r 

For  conical  frustum  measurements  in  Tunnel  A(i.  e. , Figure  31),  the  end  of  transition 
apparently  moves  forward  with  increasing  angle  of  attack.  However,  a less  clear 
trend  is  exhibited  by  the  Mach  8 data  of  Figure  35.  The  beginning  of  transition  remains 
approximately  constant  while  the  end  moves  aft  on  the  model  at  o - 14° . Contributing 
to  the  uncertain  situation  are  the  anomalous  outputs  of  sensors  A-9  (too  high)  and  A-10 
(too  low).  It  is  evident,  however,  that  peak  levels  are  5-10 dB  higher  for  the  w=  7° 
and  14°  cases  than  for  the  a = 0°  case.  Heat  transfer  data  on  the  cone  similarly  ex- 
hibit an  upward  trend  with  increasing  angle  of  attack  at  Mach  8.  This  presents  a 
different  picture  from  the  corresponding  Mach  4 runs  (Figure  31),  where  levels  of  both 
peak  acoustic  and  heat  transfer  data  show  less  dependence  on  model  angle  of  attack. 

Slice  and  flap  measurements  at  Mach  8 are  also  more  sensitive  to  cr  changes  than  are 
those  at  Mach  4 (as  a comparison  between  Figures  35  and  31  will  reveal). 

Data  measured  utilizing  the  blunt  nosetip  (R^Rg  = 0. 1)  are  shown  in  Figure  36, 
which  presents  OASPL  distributions  for  both  a blowing  and  a non-blowing  case.  At 
Mach  8 (X^  - . 018),  mass  transfer  trips  the  boundary  layer  resulting  in  a transition 
zone  on  the  conical  frustum  that  extends  approximately  from  X/L  = .2  to  .7.  In  addition, 
acoustic  levels  in  turbulent  flow  (sensors  A-19  through  A-22)  are  only  slightly  greater 

than  the  laminar  levels  indicated  by  Run  87.  However,  at  Mach  4 (X  = . 021)  virtually 
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the  entire  model  is  emersed  in  turbulent  flow,  acoustic  levels  of  which  are  approxi- 
mately 5dB  greater  than  laminar  measurements  of  Run  102  in  Tunnel  A. 

On  the  slice  at  Mach  8,  turbulent  fluctuating  pressure  levels  for  the  blunt  cone  are 
almost  identical  to  the  sharp  cone  acoustic  environment,  indicated  by  data  in  Figures 
36  and  33.  As  on  the  cone,  these  levels  are  only  slightly  higher  than  the  laminar  en- 
vironment of  Run  87,  although  thermal  data  do  indicate  increased  levels  of  heat  transfer. 
High  flap  levels  in  Run  100,  approximately  the  same  as  those  associated  with  the  sharp 
model  (Run  58),  evidently  arise  from  the  -20*  flap  deflection.  Similarity  of  charac- 
teristics between  the  blowing  case  and  the  sharp  cone  case  at  Maoh  8 is  evident  in  the 
heat  transfer  data  as  well,  wherein  levels  aft  of  the  transition  zone  and  on  the  slice 
and  flap  are  almost  identical  for  Runs  100  and  58. 


Figure  37  shows  the  results  of  a combination  of  geometric  effects.  For  the  blunt 

case,  Run  80  exhibited  the  same  characteristics  as  Run  100  along  the  conical  and  slice 

regions.  Inasmuch  as  the  flap  was  on  the  zero  angle  position,  both  the  heat  flux  and 

OASPL  there  indicated  levels  consistent  with  the  slice  values.  In  the  same  manner. 

Run  99  showed  conditions  quite  similar  to  Run  87  along  the  conical  and  slice  regions. 

It  is  interesting  to  note  the  change  in  heat  transfer  and  acoustic  characteristics  along 

the  flap  region.  Here  the  20*  flap  angle  for  the  non-blowing  nosetip  (Run  99)  shows 

levels  for  both  heat  and  OASPL  that  are  less  than  those  of  the  blowing  case  (Runs  100 

and  102).  This  phenomenon  is  somewhat  surprising  since  the  injectant  gas  at  the  nose- 

tip  should  have  negligible  effect  (compared  to  that  entrained  from  the  freestream)  at 

a downstream  position  of  x >>  R.T.  Finally,  Run  50  (<*  = 7°)  again  reflects  variations 

N 

in  the  transition  zone  when  compared  to  the  comparable  case  at  Mach  4 (Figure  28), 
relative  to  heat  flux  and  OASPL. 

(a)  2. 3 Sound  Pressure  Level  Data  - Mach  4 
Conical  Frustum  Measurements 

Figures  38  through  40  present  fluctuating  pressure  data  corresponding  to  fully 
turbulent  and  transitional  flow  at  zero  and  7*  angle  of  attack,  as  well  as  transitional/ 
turbulent  flow  at  14°  angle  of  attack.  Figure  38  demonstrates  the  essentially  similar 
spectral  content  of  the  fluctuating  pressures  in  turbulent  flow  over  a relatively  large 
distance  (X/L  - . 5 to  . 86).  In  Figure  39,  progression  of  transitional  flow  through 
onset,  peak  and  decay  into  turbulent  flow  is  presented  in  terms  of  sound  pressure  level 
(SPL)  curves  from  the  mid-cone  acoustic  sensor  array.  It  appears  that  the  maximum 
influence  of  transition  occurs  in  the  high  frequency  portion  of  the  fluctuating  pressure 
spectrum.  At  the  point  of  peak  transition,  levels  across  the  complete  spectrum  (20  Hz- 
20  KHz)  are  increased  4-15d£  when  compared  to  Tunnel  A tare  noise  and  3-10 dB  when 
compared  to  turbulent  flow.  Figure  40  illustrates  the  decay  of  transitional  flow  spectra 
at  14*  angle  of  attack.  In  this  oase,  there  appears  to  be  much  more  high  frequency 
scatter  in  the  SPL  data  for  sensors  A-3,  8,  and  9.  Levels  below  8 KHz,  however,  are 
approximately  coincident  for  these  three  sensors,  and  drop  off  several  dB  to  the  values 
indicated  by  A-10  and  A- 19. 
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Effects  of  Reynolds  number  and  model  angle  of  attack  variation  on  acoustic  levels 

forward  of  transition  (indicative  of  tunnel  tare)  are  illustrated  in  Figure  41.  These 

SPL  data  reveal  that  measured  levels  increase  with  both  Re  and  a. 

00 

(b)  Slice  Measurements 

Figures  42  through  44  present  acoustic  data  for  sensors  A-13,  14,  and  15  in  the 

slice  longitudinal  array  and  for  A-  '7  in  the  lateral  array  (A-16  and  A-18  provided 

erroneous  data  in  Tunnel  A).  Sound  pressure  level  curves  in  Figure  42  for  A-17  (Runs 

6,  23,  27,  and  42)  demonstrate  that  in  fully  turbulent  flow,  slice  acoustic  levels  are 

essentially  independent  of  flap  deflection,  for  0 s 1 6 | s 20" . In  these  zero  angle  of 
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attack  cases,  differences  evident  in  the  A-17  spectra  (for  frequencies  < 5 KHz,  there 
is  agreement  to  within  ± . 5dB)  are  a consequence  of  experimental  data  scatter,  since 
no  definitive  trends  with  flap  deflection  are  exhibited. 

Figure  43  illustrates  the  spatial  distribution  of  Mach  4 turbulent  spectra  along  the 
slice  at  a flap  deflection  of  -20°  for  the  sharp  nosetip.  Generally  decreasing  levels 
are  evident  from  slice  forward  to  slice  aft  locations.  A similar  trend  for  turbulent 
flow  was  noted  in  the  slice  spectra  for  the  blunt  nosetip  configuration  (with  blowing). 
This  situation  is  reversed,  however,  for  slice  transitional  cases.  The  wide  variation 
among  slice  SPL  data  indicated  by  Figure  44  is  typical  of  all  spectra  measured  in 
transitional  flow  over  the  slice  (Runs  69,  84,  89  and  102  - Mach  4). 

(c)  Flap  Measurements 

The  flap  was  exposed  to  turbulent  flow  conditions  for  all  Mach  4 cases  in  the  data 
sub-set  (Table  II);  given  detailed  examination  the  following  trends  were  noted  in 
flap  SPL  data: 

1.  ) At  a given  station  (e.  g. , sensor  A-2?)  and  for  identical  model 

and  freestream  conditions,  spectral  levels  increase  with  in- 
creasing flap  deflection. 

2. )  In  a given  run,  acoustic  tevele  as  measured  by  A-26,  27, 

and  28  wrre  monotonlcally  decreasing  across  the  flap. 

These  twc  trends  are  illustrated  respectively  by  Figures  44  and  45. 
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2.4  Sound  Pressure  Level  Data  - Mach  8 


2.4. 1 Comparison  of  Laminar  and  Turbulent/Transitional  Measurements 

Figure  47  presents  conical  frustum  data  in  laminar  flow  for  various  conditions  of 
freestream  Reynolds  number  and  model  angle  of  attack.  These  measurements  reveal 
the  same  trends  as  analogous  Mach  4 data  (Figure  41),  i.e. , sound  pressure  levels 

increase  with  larger  values  of  Re  and  a . Mach  8 laminar  data  also  suggest  that  a 
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higher  level  of  freestream  noise  exists  in  Tunnel  B than  in  Tunnel  A. 

A comparison  of  laminar  and  turbulent  fluctuating  pressure  data  along  the  model 
surface  at  zero  angle  of  attack  is  shown  in  Figures  48-50.  The  first  of  these  figures 
presents  Mach  8 sound  pressure  levels  measured  over  the  blunt  configuration  without 
mass  injection.  Two  characteristics  are  evident  in  these  laminar  data.  There  is 
essentially  similar  spectral  content  for  conical  frustum  measurements  (A-6  and  A-19). 
An  analogous  situation  exists  for  the  control  region  which  presents  an  effective  con- 
tinuous flat  surface  from  station  24. 25  aft  due  to  the  zero  flap  angle  for  Run  87.  A 
general  collapse  of  these  latter  data  (for  A-13,  16,  26,  27)  is  evident  between  fre- 
quencies of  250  Hz  and  12500  Hz  at  a level  approximately  5-10 dB  lower  than  the  conical 
frustum  laminar  sound  pressure  levels.  This  situation  suggests  that  the  expansion 
wave  which  exists  at  the  cone-slice  interface  may  be  responsible  for  attenuation  of  tunnel 
freestream  noise  impinging  on  the  model.  This  drop  in  level  occurs  both  longitudinally 
(compare  A-6  and  A--13)  as  well  as  circumferentially  (compare  A-19  and  A-16). 

Figures  49  and  50  (Runs  90,44,  respectively)  present  measurements  corresponding 
to  primarily  turbulent  flow  over  both  the  conical  frustum  and  the  slice/flap  region. 

(In  both  plots  sensor  A-6  is  registering  in  different  regimes  of  the  transition  zone). 
Comparing  data  for  the  two  blunt  configurations  (with  and  without  mass  injection),  it 
is  apparent  that  conical  frustum  levels  are  from  0-2  dB  greater  than  corresponding 
laminar  data  of  Figure  48.  Turbulent  levels  on  the  flap  (6p  ■ 0)  indicate  a 0-1  dB  in- 
crease over  laminar  spectra.  These  relationships  hold  except  at  the  high  frequencies 
where  turbulent  data  generally  exhibit  an  upswing.  Analogous  turbulent  measurements 
over  the  sharp  cone  model  shown  in  Figure  50  (Run  44)  are  quite  similar  to  those  of 
Run  90. 

Figures  51  and  52  respectively,  present  longitudinal  and  circumferential  variation 
of  sound  pressure  level  on  the  sharp  cooe  model  at  7*  angle  of  attaok.  The  data  of 
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Figure  51  indicate  that  the  maximum  influence  of  transition  occurs  in  the  high  fre- 
quency portion  of  the  spectrum,  an  effect  also  noted  in  the  Mach  4 sound  pressure 
level  data.  Measurements  from  A-19  in  this  figure  show  both  the  effect  of  lower 
levels  in  the  turbulent  boundary  layer  .region  and  the  effect  of  the  sensor  being  off  the 
windward  meridian.  Figure  52  indicates  the  variation  of  sound  pressure  level  in  the 
circumferential  direction  at  a = 7° . The  applicable  sensors  (i.e. , 19,  22,  21)  are 
located  at  angles  328°  , 90°  , 180° , respectively,  and  register  OASPL  values  of  130.2, 
126. 6 and  120. 5dB.  These  data  indicate  that  a minimum  of  10  dB  variation  can 
exist  between  windward  and  leeward  fluctuating  pressure  levels  at  model  angles  equal 
to  or  greater  than  7° . 

(a)  Slice  Measurements 

Figure  53  presents  slice  longitudinal  array  data  for  the  sharp  cone  at  a flap  deflec- 
tion of  -20° . (Sensors  A-17  and  A-18  in  the  lateral  array  provided  erratic  data  in 
Tunnel  B. ) These  measurements  reveal  that  a general  reduction  in  level  is  exhibited 
across  the  slice  going  aft.  However,  A-15  indicates  a much  greater  reduction  and 
appears  to  be  inconsistent  with  the  rest  of  the  data.  This  same  trend  occurs  in  Figure 
54  which  compares  sound  pressure  levels  from  A-15  and  A-16  for  various  values  of 

6 . Hence  explicit  levels  associated  with  sensor  A-15  should  be  treated  with  caution, 
r 

recognizing  that  the  measurements  are  perhaps  5-10  dB  low  at  Mach  8.  The  ?lice  data 
are  generally  insensitive  to  flap  deflection,  a fact  that  was  similarly  noted  in  the  Mach 
4 data. 

(b)  Flap  Measurements 

The  data  trend  with  flap  deflection  apparent  in  flap  acoustic  measurements  at  Maoh 
4 was  also  exhibited  by  sound  pressure  level  data  at  Mach  8.  Figure  55  reveals  that 
flap  acoustic  levels  increased  by  approximately  25  dB  throughout  the  frequency  range 
when  was  increased  from  0 to  -20* . At  Maoh  8,  the  general  trend  in  spectrum 
level  with  flap  location  was  to  increase  across  the  flap  going  aft.  However,  OASPL 
values  do  not  necessarily  follow  this  trend  due  to  the  difference  in  high  frequency  data 
for  sensors  A-27  and  A-28.  This  effect  is  Indicated  by  the  data  of  Figure  56. 
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2.5  Root-Mean-Square  Pressure 

Figure  57  illustrates  rms  fluctuating  pressure  data  along  the  conical  portion  of 
the  model  for  both  test  Mach  numbers.  The  attached  turbulent  boundary  layer  pre- 
dictions of  Lowson  and  the  present  method  (Equation  (61);  are  shown  for  comparison 
as  well  as  the  transitional  boundary  layer  prediction  of  Reference  1.  With  exception 
of  the  mass  injection  cases,  the  conical  region  acoustic  gages  were  experiencing 
transitional  flow  characteristics,  thereby  limiting  data  from  fully  developed  turbulent 
flow  conditions.  Also  shown  in  the  figure  are  typical  normalized  rms  pressure  levels 
in  laminar  flow  (indicative  of  tunnel  noise).  Higher  values  of  rms  pressure  are  noted 
for  the  blunt  configuration,  a consequeuce  of  the  lower  edge  dynamic  pressure  (entropy 
effect).  The  data  are  generally  in  g<>od  agreement  with  other  reported  results  for 
zero  angle  of  attack  conditions;  in  particular,  the  data  of  Reference  1 which  was  ac- 
quired in  the  same  facility.  Increasing  angle  of  attack  tends  to  introduce  significant 
data  scatter,  particularly  for  transitional  flow.  The  data  appear  to  indicate  decreasing 
values  of  normalized  rms  pressure  with  increasing  a.  It  is  difficult  to  discern  Rey- 
nolds number  effect  with  the  small  data  sample  available.  The  blunt  nose  data  repre- 
sent conditions  where  mass  injection  was  introduced  at  the  nosetip  for  the  purpose  of 
inducing  turbulent  flow  characteristics.  Here,  both  transitional  and  turbulent  data 
appear  to  be  consistent  with  sharp  measurements  acquired  through  natural  transition. 

Normalized  rms  pressures  for  the  slice  region  are  shown  es  a function  of  local 
Mach  number  in  Figure  58.  The  data  for  both  test  Maoh  number  conditions  appear  to 
be  in  general  agreement  with  other  reported  results  (Figure  1)  and  show  a continued 
decrease  in  Prmg/( for  values  ot  Me>6.  Also  shown  in  this  figure  are  the  flap  sensors 
A-26  and  A-27  for  zero  flap  deflection  (continuation  of  the  slice).  There  appears  to 
be  a decrease  in  rms  pressure  with  increasing  angle  of  attack  as  experienced  on  the 
conical  region.  However,  the  data  (at  least  tor  M = 8)  appear  to  bo  consistent  with 
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prediction  for  this  condition.  Slice  measurements  appear  to  reflect  a Refolds  number 

effect,  such  that  higher  levels  of  normalized  rms  pressure  occur  with  decreasing 

87 

freestresm  unit  Reynolds  number.  Donaldson  has  noted  higher  free  stream  noise 
content  at  this  unit  Reynolds  number  condition  than  for  the  other  reported  value 

Re  /ft.  - 3.6  x 10"6). 
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It  is  evident  that  measured  laminar  boundary  layer  levels  on  the  slice  (Run  87, 

M =8)  are  lower  than  the  conical  values,  a fact  similarly  noted  in  the  discussion 

CD 

of  SPL  in  Section  2. 4 above.  This  would  indicate  that  the  tunnel  radiated  noise  is 
attenuated  as  a result  of  the  expansion  wave  from  the  geometry  change  of  the  cone  to 
the  slice. 

Normalized  rms  pressure  fluctuations  for  the  flap  environment  are  shown  in 

Figure  59.  Increasing  the  flap  deflection  increases  the  rms  pressure  for  both  zero 

and  non-zero  angle  of  attack  conditions.  It  is  quite  interesting  to  note,  however, 

that  for  zero  angle  of  attack  and  0 s -f.  £ 15°  data  tend  to  follow  the  pattern  ex- 
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perienced  on  the  cone  and  slice.  The  most  significant  depature  appears  to  be  for  flap 
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deflections  >-15° , lower  freestream  unit  Reynolds  number  (i.  e. , 2. 5 x 10  ),  and 
blunt  body  tests.  Concerning  the  bluntness  effect,  the  non-injectant  data  tend  to  yield 
higher  values  of  rms  pressure  than  the  injectant  cases.  The  former  might  be  ex- 
periencing transitional  effects  in  the  flap  region.  It  is  also  interesting  to  note  that 
the  entropy  layer  significantly  reduces  the  local  pressure  field  for  the  Mach  8 condition 

(-6  > 10° ) compared  to  a negligible  effect  at  the  Mach  4 condition.  This  effect  ap- 
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pears  evident  in  the  data,  such  that  the  rms  pressure  at  Mach  8 tend  to  be  lower  than 
those  measured  at  Mach  4.  One  must  also  keep  in  mind  that  the  shock  emanating  from 
the  hinge  joint  together  with  the  expansion  wave  from  the  conlcal/sllce  intersection 
can  interact,  causing  high  rms  pressure  levels  at  the  extended  flap  positions.  Although 
not  an  obvious  trend,  increasing  model  angle  of  attack  tends  to  decrease  the  rms  pres- 
sure level.  However,  as  in  the  slice  oase,  flap  data  tend  to  follow  the  Mach  number 
dependent  predictions  indicated  in  Figure  59. 


2.6  Power  Spectral  Penalty  - Mach  4 


(a)  Conical  Frustum  Measurements 

Figure  60  presents  normalized  power  spectral  density  as  a function  of  Strouhal 
number  for  the  Mach  4 test  conditions.  Normalizing  parameters  can  be  obtained 
from  Tables  m and  V or  from  third  octave  PSD  plots  of  actual  data  presented  in 
Volume  n.  The  data  are  generally  separated  into  two  categories  consisting  of  transi- 
tional and  turbulent  flow.  For  transitional  flow,  the  data  reveal  an  angle  of  attack 
effect,  l.e. , a decrease  in  normalized  PSD  with  increasing  a.  The  same  trend  is 
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exhibited  for  the  turbulent  flow  behavior,  in  particular  for  f6*/u  > 10  . This  latter 
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observation  should  be  considered  with  caution  inasmuch  as  the  reporting  sensor 
(A-19)  is  not  in  line  with  the  windward  meridian  for  which  local  flow  properties  at 
angle  of  attack  were  evaluated.  The  turbulent  boundary  layer  data  also  indicate  a 
possible  bluntness  effect  which  is  believed  to  be  a consequence  of  the  entropy  layer 
yielding  lower  values  of  dynamic  pressure.  Insufficient  data  do  not  allow  for  specula- 
tion of  Reynolds  number  effects.  In  general,  these  normalized  spectral  measurements 
agree  with  data  obtained  in  a previous  study1  in  the  same  facility. 

(b)  Slice  Measurements 

Normalized  power  spectral  density  data  for  the  slice  portion  of  the  model  are  shown 
in  Figure  61.  These  measurements  reflect  turbulent  boundary  layer  conditions  for  the 
high  unit  free  Reynolds  number  cases.  Expressed  in  this  format,  the  data  do  not  ap- 
pear to  indicate  any  specific  trends  with  angle  of  attaok.  However,  bluntness  effects 
are  reflected  by  the  higher  levels  in  PSD.  In  general,  the  magnitude  of  normalized 
PSD  along  the  slice  tends  to  be  less  than  values  measured  on  the  conical  frustum.  As 
a consequence  of  flow  expansion  resulting  from  the  conical  to  slice  geometry  change, 
the  local  pressure  decreases  and  subsequently  could  attenuate  the  fluctuating  pressure 
field.  One  should  also  keep  in  mind  that  the  expansion  wave  has  the  effect  of  also 
attenuating  the  tunnel  radiated  noise  thereby  rendering  less  contaminated  noise  It  els 
on  the  slice  region. 

(c)  Flap  Measurements 

Figures  62  and  63  display  the  normalised,  power  spectra  for  the  flap  region  of  the 
model.  In  Figure  62,  a significant  Reynolds  number  effect  is  exhibited.  Two 
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phenomena  relative  to  this  condition  should  be  noted.  First,  for  Re  /ft.  + 2.  5 x 10  . 
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possible  shock  interaction  between  the  weak  expansion  wave  (cone/slice  intersection) 
and  strong  hinge  shock  may  be  contributing  to  measured  values.  Second,  the  lower 
freestream  Reynolds  number  runs  exhibited  transition  zones  over  the  control  surface 
region.  Figure  29  shows  the  OASPL  for  this  Reynolds  number  condition  to  be  greater 
than  other  reported  measurements  when  the  slice  was  experiencing  turbulent  flow 
characteristics. 

An  examination  of  the  zero  flap  deflection  condition  indicates  no  discernable  angle 
of  attack  effects.  However,  with  increasing  flap  angles,  a slight  effect  of  a is  noted, 
such  that,  within  a given  flap  angle  data  set  (6  * 0),  an  increase  in  angle  of  attack 
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tends  to  decrease  the  normalized  power  spectral  density.  A significant  increase  in 
normalized  PSD,  however,  is  experienced  with  increasing  flap  angle.  This  condition 
can  be  expected  as  a consequence  of  higher  local  pressure  and  subsequently  increased 
rms  fluctuating  pressure  levels.  It  should  be  noted  that  the  zero  flap  angle  condition 
(which  is  a continuation  of  the  slice)  yields  levels  that  are  consistent  with  those  ob- 
tained on  the  conical  frustum  as  well  as  the  slice.  However,  the  higher  flap  angle 
conditions  yield  levels  of  normalized  PSD  up  to  an  order  of  magnitude  higher.  This 
situation  mandates  a modification  to  the  prediction  algorithms  developed  in  the  litera- 
ture for  more  conventional  geometric  configurations. 

Figure  63  shows  the  normalized  PSD  distribution  along  the  flap  surface  for  a 
6 = -20°  condition  at  the  higher  unit  freestream  Reynolds  number  condition  (3.6  x 
10  ).  Data  from  both  sharp  and  blunt  nosetip  configurations  are  represented.  The 
latter  features  conditions  with  and  without  mass  Injection  at  the  nosetip.  While  the 
data  reflect  a significantly  higher  level  than  that  experienced  on  the  cone  or  slice  region, 
a still  further  departure  is  exhibited  for  the  non-blowing  blunt  case  (Run  102).  Here, 
the  normalized  PSD  level  appears  to  be  in  the  range  of  data  reported  for  lower  free- 
stream unit  Reynolds  number  of  Figure  62  which  featured  transitional  flow  along  the 
slice.  An  examination  of  Figure  32  indicates  that  the  non-blowing  blunt  run  is  also 
experiencing  transitional  flow  along  the  slice  and  has  OASPL's  exceeding  values  where 
turbulent  flow  prevailed  on  the  slice.  The  condition  when  the  slice  experiences  tur- 
bulent boundary  layer  flow,  as  in  Runs  42  (sharp)  and  103  (blunt  with  mass  injection). 
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normalized  flap  measurements  fall  generally  into  the  same  range  (Figure  63).  In  order 
to  ensure  that  the  effect  being  observed  is  a consequence  of  transitional  flow  along 
the  slice,  an  examination  was  made  relative  to  the  shock  angle  interaction.  The  bow 
shook  angle  15. 7s ) did  not  intersect  with  either  the  expansion  wave  (conical/slice) 
or  flap  shock  for  the  sub-set  conditions  analyzed  in  this  study  (Table  Q).  Hence,  the 
high  levels  of  normalized  PSD's  shown  in  Figures  62  and  63  appear  to  be  a consequence 
of  transition  occurring  on  the  slice  foward  of  the  flap. 

2.7  Power  Spectral  Density  - Mach  8 

(a)  Conioal  Frustum  Measurements 

Figure  64  presents  normalized  power  spectral  density  as  a function  of  Strouhal 
number.  As  in  the  Mach  4 tests,  data  generally  separate  into  two  categories  corres- 
pond^ to  transitional  and  turbulent  boundary  layer  flow  conditions.  Moreover,  the 
blunt  nose  run  (100)  features  higher  levels  of  nor  nallzed  PSD  for  turbulent  flow,  a 
condition  that  also  existed  at  Mach  4 (see  Figure  60).  However,  it  should  be  noted 
that  the  transitional  da«:a  reflect  an  angle  of  attack  trend  that  is  opposite  to  that  exhibited 
at  Mach  4;  i.  e.  , an  increase  in  level  with  increasing  a.  For  turbulent  flow  conditions, 
no  discemable  trend  in  normalized  PSD  with  angle  of  attack  is  evident.  The  data  are 

in  general  agreement  with  results  of  the  Mach  4 tests  as  well  as  data  with  other  transi- 
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ttonal  and  turbulent  reported  in  the  literature.  ' ' 

(b)  Slice  Measurements 

Normalized  power  spectra  for  the  slice  sector  of  the  model  are  shown  in  Figures 
65  and  66.  Unlike  the  Mach  4 data,  a more  significant  variation  in  the  zero  intercept 
values  (flat  portion  of  spectra)  is  evident.  The  level  of  sensor  A-15  appeared  to  be 
much  lower  than  measurements  for  comparable  gages  along  the  slice  (see  2. 3 above). 
Althot«h  a careful  examination  of  the  raw  data  was  made,  which  did  not  include  any 
malfunction,  the  absolute  data  for  this  sensor  are  considered  of  dubious  validity. 
However,  the  general  trend  of  the  measurements  relative  to  bluntaess,  angle  of  attack, 
and  Reynolds  number  can  be  used  to  determine  fluid  dynamic  effects.  The  blunt  nose- 
tip  runs  again  reflect  the  significant  effects  of  transition.  Run  102  represents  s con- 
trolled test  that  featured  a transition  zone  over  the  conical  frustum  of  the  model. 
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Figure  36  show-  the  transition  zones  for  both  Runs  100  and  102  based  on  thermal 
measurements.  Here  a longer  transition  zone  yielded  higher  values  of  OASPL  and 
ultimately  affected  flap  readings  as  well.  It  should  be  noted  that  the  predicted  end  of 
transition  for  both  runs  was  approximately  the  same  (x«  23  inches);  however,  the 
heat  transfer  data  of  Run  102  were  much  more  difficult  to  assess.  The  sharp  cone 
■•alues  of  normalized  PSD  over  the  slice  (Run  58)  indicate  levels  that  are  lower  than 
the  Mach  4 resultf . 

Figure  66  displays  normalized  power  spectral  density  at  several  conditions  for 
sensor  /. -15  (sharp  cone  model).  The  most  significant  effect  exhibited  by  these  slice 
data  appears  to  be  tne  angle  of  attack  variation  for  the  zero  intercept  value.  More- 
over, a comparison  of  the  a = 0“  data  indicate  higher  levels  of  nor  malized  PSD  for 
the  lower  freestream  unit  Reynolds  number,  a result  that  is  consistent  with  the  Mach 
4 data.  However  nc  dispai  lty  in  angle  of  attack  was  observed  at  the  Mach  4 condition 
as  noted  at  Mach  8. 

'c)  Flap  Measurements 

Figures  67  and  68  present  normalized  PSD's  for  the  fla*.  region  of  the  model.  As 
for  the  Mach  4 runs,  a significant  effect  of  Reynolds  number  is  notrd  (l.e. , Run  77 
at  Re^/ft.  - 2.  5 x 10  ).  However,  unlike  the  Mach  4 tests,  the  slico  at  Mach  8 was 
not  exposed  tc  '.ransitional  flow  characteristics.  Moreover,  the  zero  intercept  value 
appears  to  be  lower  than  the  comparable  Mach  4 case  (Run  84).  It  is  interesting  to 
note  the  disparity  of  angle  of  attack  effects  for  zero  flap  deflection.  There  is  a marked 
increase  in  normalized  spectral  levels  with  increasing  a which  did  not  occur  at  Mach 
4.  Data  for  a * 0*  reveal  zero  intercept  values  that  are  an  order  of  magnitude  less 
than  the  Mach  4 test  (Run  6).  The  data  also  show  the  same  trend  as  Mach  4 measure- 
ments with  increasing  flap  deflection.  For  the  20*  flap  deflection,  normalized  PSD 
levels  are  generally  of  the  same  magnitude.  Finally,  it  should  be  noted  that  normalized 
levels  on  the  flap  at  zero  deflection  (continuation  of  the  slice)  indicate  higher  values 
than  the  slioe  (Figure  60).  It  la  speculated  that  the  field  joint  (flap  hinge)  might  act 
as  a roughness  mechanism  which  in  turn  enhances  the  fluctuating  pressure  field. 
However,  this  situation  appeared  to  occur  only  for  the  Mach  8 teats. 

Figure  68  displays  normalized  power  spectral  density  along  the  flap  surface  for 
the  higher  unit  freestream  Reynolds  number  and  for  20“  flap  angle.  Levels  of  normalized 
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PSD  are  generally  the  same  as  those  experienced  at  Mach  4 (Figure  63).  One 
will  recall  that  Run  102  (Mach  4)  exhibited  high  values  of  PSD  as  a consequence  of 
transitional  flow  characteristics  along  the  slice,  thereby,  enhancing  the  flap  values. 

As  previously  indicated  Run  102  (Mach  8)  exhibited  a transition  zone  over  the  conical 
frustum  and  subsequently  influenced  the  slice  and  flap  region  OASPL's  (see  Figure 
36). 

2.8  Broad-Band  Cross-Correlation  Functions 
(a)  Conical  Frustum  Measurements 

Figures  69  and  70  present  longitudinal  broad-band  correlation  functions  for  laminar 
transitional  and  turbulent  flow.  The  laminar  data  at  Mach  4 and  8 provide  an  indication 
of  the  correlation  characteristics  of  the  tare  acoustic  environment  present  in  Tunnels 
A and  B respectively.  For  Mach  4 and  Mach  8 transitional  data  over  the  conical  frustum 
levels  of  correlation  are  generally  compatible  with  results  presented  in  Reference  1, 
even  to  the  extent  that  Mach  8 levels  are  higher  than  those  of  Mach  4.  Current  data, 
however,  exhibit  a much  more  marked  decrease  with  increasing  model  angle  of  attack. 
This  may  be  due  to  shifting  of  the  transition  zone  relative  to  the  conical  sensor  array 
with  changes  in  model  orientation.  Longitudinal  correlation  data  corresponding  to  tur- 
bulent flow  on  the  cone  (Rims  103  and  10  - Tunnel  A)  are  shown  in  Figure  69,  and  re- 
sults demonstrate  quite  good  agreement  with  turbulent  data  of  Reference  1.  Reduced 
levels  of  correlation  present  in  Run  10  are  evidently  related  to  the  fact  that  the  acoustic 
environment  is  in  an  intermediate  state  between  transitional  and  fully  turbulent  flow 
(see  Figure  28).  In  the  data  of  Run  103,  turbulence  was  induced  through  tripping  the 
boundary  layer  utilizing  mass  injection  over  the  nosetip,  whereas  transition  occurred 
naturally  in  the  1972  program1  and  turbulent  measurements  were  taken  at  model  loca- 
tions aft  of  the  transition  zone. 

Figure  71  presents  circumferential  cross -correlation  functions  relating  sensors 
A-5  and  A-6  only,  due  to  the  fact  that  A-7  did  not  provide  valid  data.  As  in  the  case 
of  longitudinal  functions,  transitional  data  exhibit  slightly  higher  correlation  than  tur- 
bulent, and  Mach  8 transitional  data  are  higher  than  analogous  Mach  4 circumferential 
functions.  At  Mach  8,  there  was  little  difference  between  circumferential  correlation 


79 


values  obtained  over  the  blunt  cone  with  mass  injection,  hence  only  examples  of  this 
former  are  shown  in  Figure  71. 

(b)  Slice  Measurements 

Figure  72  presents  Mach  4 longitudinal  correlation  functions  applicable  to  transi- 
tional and  turbulent  flow  on  the  slice.  For  the  latter  data  at  zero  angle  of  attack, 
there  is  only  slight  variation  in  the  slice  dross -correlation  characteristics  between 

6 =0°  and  6_  = -20® . It  is  apparent  that  increasing  the  model  angle  of  attack  tends 
r r 

to  increase  both  the  rate  of  decay  with  separation  distance  and  the  time  of  maximum 
correlation.  These  effects  were  also  noted  in  measurements  over  the  conical  frus- 
tum (Figure  69). 

Longitudinal  cross-correlation  functions  measured  at  Mach  8 are  shown  in  Figure 
73.  Laminar  flow  occurs  over  the  slice  for  the  blunt  nosetip  configuration  without 
mass  injection  at  the  nosetip.  A typical  set  of  broad-band  functions  corresponding  to 
laminar  flow  is  presented  and  values  of  correlation  are  seen  to  be  generally  higher 
than  turbulent  levels.  This  phenomenon  was  also  evident  in  the  conical  frustum  data 
of  Mach  4.  For  the  Mach  8 condition,  it  i.8  apparent  that  for  Run  100  (where  turbulent 
flow  is  induced  by  mass  injection)  the  rate  of  cross-correlation  decay  is  greater  than 
for  those  cases  wherein  natural  transition  occurs  forward  of  the  slice  at  zero  angle 
of  attack  (e.g. , Run  58).  As  is  the  case  of  Mach  4 turbulent  data,  the  rate  of  corre- 
lation decay  with  separation  distance  increases  for  o = 7°  and  14° . 

(o)  Flap  Measurements 

Longitudinal  cross -correlation  functions  measured  on  the  flap  are  presented  in 
Figures  74  and  75  for  Mach  4 and  Mach  8,  respectively.  The  initial  set  of  functions 
in  Figure  74  corresponds  to  those  Tunnel  A cases  wherein  the  slice  environment  was 
designated  as  transitional  (Runs  69,  84,  and  102).  Flap  overall  sound  pressure  level 
for  these  cases  (Figures  26  and  29)  are  the  highest  of  all  acoustic  measurements  in 
either  Tunnel  A or  B,  (f>  & 20°)  and  the  environments  are  further  characterized  by 
relatively  large  differences  in  level  among  sensors  A-26,  -27,  and  - 28  (AdBs»  12 
between  A-26  and  A-28).  Corresponding  cross -correlation  functions  shown  in  Figure 
74  depict  a rapid  decay  in  correlation  over  the  flap.  These  phenomena  suggest  a type 
of  Incipient  transitional/turbulent  flow  may  be  impinging  on  the  flap  for  these  oases. 
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analogous  to  the  "overshoot"  effect  on  acoustic  level  exhibited  by  trmsitlonal  flow  on 
a cone. 

Also  shown  in  Figure  74  are  data  indicating  the  effect  of  flap  deflection  on  longi- 
tudinal correlation  functions.  It  is  evident  that  the  rate  of  decay  with  separation 
distance  increases  for  increasing  flap  angle.  The  final  data  of  Figure  74  consist  of 
turbulent  cross-correlations  which  show  that  the  respective  effects  of  high  model 
angle  of  attack  and  mass  injection  induced  turbulence  are  apparently  to  increase  the 
rate  of  correlation  decay  with  separation  distance. 

Cross -correlation  data  for  the  flap  in  Run  87  (Mach  8 - R^/Rg  = 0. 1)  are  presented 
in  Figure  75.  Heat  transfer  data  as  shown  in  Figure  36  suggest  that  with  no  mass  in- 
jection over  the  nosetip,  the  boundary  layer  over  the  flap  is  laminar  or  perhaps  in  an 
early  transitional  state.  As  in  previous  conical  frustum  and  slice  data  corresr  ending 
to  laminar  or  transitional  boundary  layer  flow,  relatively  high  values  of  correlation 
are  indicated.  For  the  sharp  cone  at  Macb  8,  generally  lower  values  of  correlation 
over  the  flap  occurred  at  the  higher  flap  deflection  angles,  following  the  trend  for  Mach 
4 data.  The  final  data  of  Figure  75  consist  of  Mach  3 flap  cross-correlations  corre- 
sponding to  the  sharp  cone  model  at  14°  angle  of  attack  and  to  the  blunt  cone  model  mass 
injection  induced  turbulence. 

2. 9 Broad-Band  Convection  Velocity 

Values  of  broad -band  convection  velocity  were  computed  for  array  sensors  in  the 

conical,  slice,  and  flap  regions  of  the  model  following  procedures  outlined  in  Section 

V.  Based  on  conical  frustum  results  obtained  in  the  1972  wind  tunnel  program  at 

AEDC, 1 it  was  anticipated  that  transitional  convection  velocities  expressed  as 

“c  nr/u®  would  be  generally  loss  than  those  measured  in  turbulent  flow,  for  all  por- 
BB 

tions  of  the  model.  However,  data  obtained  in  the  current  program  do  not  follow  that 
relationship,  and  actually  exhibit  a wide  degree  of  scatter.  Indeed,  there  appear  to 
be  no  definitive  trends  with  model  geometry  or  with  boundary  layer  flow  characteristic 
and  normalized  convection  velocities  vary  between  . 5 and  1. 0. 
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SECTION  V 


DATA  PROCESSING  AND  ANALYSIS 
1.0  DATA  RECORDING  AND  PROCESSING 

A brief  description  of  the  recording  procedures  utilized  in  this  teat  program  is 
included  here.  For  additional  details  the  reader  is  referred  to  Volume  II  (Classified). 
The  transducer  channel  recording  network  is  shown  in  Figure  76.  Signals  from  the 
thirty  three  accustic  sensors  were  routed  through  multi-conductor  twisted  pair 
cables  from  the  model  to  a junction  box  outside  the  wind  tunnel  where  they  interfaced 
with  line  drivers  provided  by  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL) 
mobile  data  van.  The  microphone  acoustic  signals  were  further  amplified  within  the 
van  and  recorded  using  three  rack-mounted,  fourteen  track  (IRIG- format),  FM 
Record/Reproduce  Magnetic  Tape  Recorders.  The  data  were  recorded  at  a tape  speed 
of  60  inches  per  second  using  Medium  Band  FM  recording  to  achieve  frequency  response 
from  20  Hz  to  20  kHz  (limited  at  the  low  end  by  the  amplifiers  on-board  the  model). 
Magnetic  tape  format  was  set  up  to  assure  that  all  data  from  a given  microphone  array 
were  recorded  on  tracks  common  to  a particular  recording  head  stack  so  that  Bub- 
sequent  cross-correlation  could  be  done  with  as  little  attendant  phase  variation  as 
possible. 

When  the  wind  tunnel  model  in  the  derived  geometric  configuration  had  been  inserted 
into  the  flow  to  achieve  a specified  environment,  and  steady  state  aerodynamic  condi- 
tions were  established,  a sample  of  the  output  from  each  acoustic  sensor  was  recorded 
for  a minimum  of  5 seconds.  These  samples  of  fluctuating  pressure  versus  time  for 
the  data  subset  of  Table  n were  then  subjected  to  the  spectral  and  correlation  analyses 
shown  in  Table  VU  and  Figure  77,  using  the  digital  processing  system  at  AFFDL.  Both 
sound  pressure  level  and  power  spectral  density  were  obtained  for  each  desired  sensor 

output  using  a one  third  octave  band  analysis.  High  frequency  correlations  (a  5kHz) 
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to  these  data  were  included  on  the  baslB  of  the  Corcos  sensor  size  criterian.  The 
maximum  correction  so  applied  was  3 dB  at  20  kHz  for  Tunnel  A data.  Cross-spectral 
properties  (e.g. , normalized  co-  and  quad- spectrum,  phase  angle,  narrow-band  con- 
vection velocity,  and  coherence)  were  also  processed  as  functions  of  third-octave 
center  frequency  for  selected  cone,  slice  and  flap  sensor  arrays.  Broad-band  cross- 
correlation functions  were  computed  by  performing  the  inverse  Fourier  transform  of 


the  cross-power  spectral  density.  Values  of  broad-band  convection  velocity  were  cal- 
culated both  through  the  appropriate  relationship  expressed  in  Table  VQ,  and  by 
dividing  the  sensor  separation  distance  by  the  time  delay  corresponding  to  the  cross- 
correlation maximum  value.  Generally  close  agreement  between  the  two  methods 
was  achieved. 

2.0  DATA  ANALYSIS 

The  objective  of  this  investigation  is  to  develop  aeroacoustic  load  predictions 
associated  with  maneuvering  type  vehicles;  in  particular  along  control  surface  regions. 
An  experimental  and  subsequent  analytical  study  was  conducted  for  this  purpose.  The 
approach  adopted  considered  a careful  assessment  of  existing  prediction  capability 
relative  to  attached  turbulent  boundary  layer  flow  characteristics.  Moreover,  the  use 
of  incompressible  flow  data,  which  has  been  well  documented  and  examined  in  the 
scientific  community,  was  used  as  a baseline  to  test  the  various  prediction  schemes. 
Having  established  the  validity  of  the  prediction  technique  relative  to  incompressible, 
attached  turbulent  flow  behavior,  the  techniques  were  assessed  relative  to  compres- 
sible flow  conditions.  Inasmuch  as  the  techniques  appear  valid  for  various  flow  (in- 
compressible/compresjible)  and  geometric  conditions  (flat  plates,  cones,  wind  tunnel 
walls,  ogive  cylinder,  wings,  etc.),  it  appears  reasonable  to  expect  their  extension 
to  complex  flow  regions  featuring  control  surfaces.  One  must  keep  in  mind  that  attached 
turbulent  boundary  flow  is  the  primary  consideration. 

Section  III  of  this  report  presents  the  methodology  adopted  to  predict  environmental 
aeroacoustic  loads.  The  present  section  will  consider  these  techniques  together  with 
data  presented  in  Section  IV  as  the  basis  for  modification  of  present  capability  as  well 
as  extension  to  the  control  surface  area.  One  must  keep  in  mind  that  a priori  selec- 
tion of  the  data  sub-set  (Table  II)  was  influenced  by  realistic  angle  of  attack-flap  angle 
conditions.  Consequently,  limitations  in  data  restrict  interpretation  of  fluid  dynamic 

(a.  Re,  M,  T /T  ) and  geometric  conditions  (RM,  5 ).  This  section  will  also  con- 
w aw  w l? 

aider  use  of  normalizing  parameters  other  than  those  presented  in  Seetion  IV.  The 
various  boundary  layer  parameters  that  have  been  considered  in  the  literature  were 
discussed  in  Section  m.  For  this  purpose,  several  select  tests  will  be  used  where 
the  fully  developed  turbulent  boundary  layer  flow  prevailed  along  the  model  at  zero 
angle  of  attack  conditions. 
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2.1  Root-Mean-Square  Pressure 


Figures  57  through  59  present  the  rms  fluctuation  pressure,  normalized  with 

dynamic  pressure,  as  a function  of  Mach  number.  Inasmuch  as  the  dynamic  pressure 
. 2 

(q  = y/2  P M ) contains  the  Mach  number,  its  use  as  a normalizing  parameter 
6 6 6 

(when  plotting  as  a function  of  Mach  number),  although  apparently  valid,  is  question- 
able. When  the  local  static  pressure  was  used,  Equation  (64),  shown  in  Figure  29, 
it  also  appeared  valid  relative  to  the  data  required  prior  to  this  study.  Inasmuch  as 
the  prediction  technique  for  the  rms  pressure  was  developed  on  the  premise  that  the 
ratio  ranged  from  1.7  < prms/Tw  < 3 for  incompressible  flow  to  1.2  <Prms/Tw  <5.6  for 
compressible  flow,  ’ the  current  data  were  plotted  using  wall  shear  stress  and 
are  shown  in  Figure  78  as  a function  of  Mach  number.  The  data  show  a range  of  the 

ratio  0.  5 < p/r  <4  which  includes  turbulent  and  transitional  flow  behavior  for  both 
w 

the  conical  frustum  and  control  surface  region.  Also  shown  is  the  freestream  noise 

95 

level  measured  by  Laderman  in  the  AEDC  Tunnel  B facility. 

Some  interesting  results  are  observed  for  the  data  presentation  of  Figure  78.  It 

is  noted  that  the  conical  frustum  data  indicate  normalized  rms  pressure  levels  less 

than  unity.  The  same  phenomenon  Is  also  experienced  on  the  slice  and  flap  regions. 

These  measured  values  on  the  model  surface  are  less  than  the  freestream  noise  con- 

33 

tent  measured  by  Laderman.  Considering  the  works  of  Stainback  et  al  as  well  as 
98 

Beckwith,  who  compared  hot  wire  fluctuating  pressure  measurements  to  surface 
sensor  values  in  laminar  boundary  layers,  there  does  not  appear  to  be  a significant 
attenuation  effect  of  freestream  noise  through  the  shock  layer.  The  current  data  do 
not  exhibit  this  behavior;  in  particular,  data  representing  zero  angle  of  attack,  flap 
deflections  < 7°  and  sharp  cone  configurations.  For  p/t  > * » these  values  signify 
high  flap  deflections  > 7° , bluntness,  and  the  lower  freestream  unit  Reynolds  number 

g 

tests  (i.e. , 2.  5 x 10  ).  Inasmuch  as  the  wall  shear  values  are  consistent  with  values 
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obtained  in  previous  studies  ' at  AEDC,  and  considering  the  fact  that  measured 
levels  of  normalized  rms  pressure  in  laminar  flow  are  higher  than  turbulent/transi- 
tional levels,  one  must  speculate  on  possible  shock  attenuation  effects.  The  laminar 
flow  data  (indicative  of  freestream  noise)  shown  in  Figure  78  are  twice  as  high  as 
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reported  by  Laderman  (Tunnel  B).  Further  assessment  of  freestream  noise  influence 
on  reported  data  will  be  made  later  in  this  section. 

Figures  79  through  81  display  the  rms  fluctuating  pressure  normalized  with  the 
static  value  as  a function  of  local  Mach  number.  Figure  79  shows  the  data  on  the 

conical  frustum  for  both  Tunnel  A and  Tunnel  B test  conditions.  The  prediction  schemes 

55 

of  the  present  study  as  well  as  that  of  Lowson''  are  also  shown  for  comparison.  Al- 
though the  Lowson  concept  tends  to  fit  the  data  for  attached  turbulent  flow  along  the 
conical  frustum,  its  validity  is  deemed  inadequate  when  comparing  to  other  data  in 
the  literature  (Figure  20). 

When  the  data  are  presented  in  the  coordinates  of  Figure  79,  both  angle  of  attack 

and  Reynolds  number  effects  appear  distinquishable.  The  same  effect  was  noted  when 

using  the  dynamic  pressure  for  normalizing  the  rms  fluctuating  pressure.  Keeping 

in  mind  that  the  data  are  limited  and  allowing  for  speculative  license,  one  might 

account  for  the  angle  of  attack  effect  through  division  by  a parameter  such  as  ( 1 + 

a/9  ),  i.e. , in  Equations  (59)  or  (64).  Figure  79  shows  this  result  for  the  two  angle 
c 

of  attack  conditions.  In  each  case,  the  data  are  overpredicted;  however,  noting  the 
low  values  of  normalized  rms  pressure  indicates  the  possibility  of  the  approach. 

Figure  80  shows  normalized  rms  pressure  data  for  the  slice  region  of  the  model. 

The  high  Mach  number  data  appear  to  be  consistent  with  prediction.  Although  angle 
of  attack  effects  are  noted  (in  particular,  for  the  Mach  4 data),  the  data  tend  to  follow 
prediction  in  a more  orderly  fashion  than  those  measured  on  the  conical  frustum. 

One  must  keep  in  mind  again,  that  the  slice  measurements  indicated  lower  levels  of 
rms  pressure  for  laminar  boundary  layer  conditions  than  on  the  conical  frustum. 
Moreover,  the  larger  boundary  layer  thicknesses  experienced  on  the  slice  would  tend 
to  yield  smaller  errors  due  to  sensor  size  phenomena.  These  data  together  with  the 
data  of  Figure  20  lend  definite  support  to  the  present  prediction  method  (Equations  (59) 
or  (64)). 

Figure  81  presents  the  normalized  rms  fluctuating  pressure  for  the  flap  region. 

As  might  be  expected,  the  data  exhibit  the  same  variation  when  normalizing  with  dy- 
namic pressure  (Figure  59).  For  zero  angle  of  attack  and  flap  deflections  in  the  range 
0 < 5 < -15° , the  data  generally  exhibit  the  same  trend  as  slice  data.  If  we  again 
allow  for  possible  speculation,  the  higher  flap  deflections  might  be  accounted  by 
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multiplying  Equations  (59)  and  (64)  by  a function  suoh  as  (1+  £_/©  ) . Utilizing  such 

r C 

an  approach,  results  are  shown  in  Figure  81  for  three  flap  angles.  The  20“  flap 

deflection  happens  to  be  coincident  with  the  T «T  curve.  As  previously  noted, 

w aw 

the  highest  reported  rms  pressure  values  occurred  when  the  slice  experienced 
transition  flow  or  in  the  test  where  transition  was  developed  over  the  entire  conical 
frustum. 

When  assessing  the  prediction  techniques  for  rms  fluctuating  pressure.  Equation 
(71)  is  recommended;  namely 


where 
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In  Section  HI.  we  have  noted  that  the  velocity  power  law  exponent  (n)  does  not  influence 

the  ratio  K(n)  RF(n)/Re  2^+n\  Moreover,  attempting  to  define  the  value  of  A (1. 7 < 

s 

Prm/Tw=  A<3‘  2)from  incompressible  flow  was  not  possible.  Also,  it  was  determined 

that  the  incompressible  value  of  the  normalized  rms  pressure  p /q  ■ 0. 006  adequately 

rms  e 

represented  the  data  and  supported  the  various  prediction  schemes  for  power  magni- 
tude and  power  spectral  density.  Consequently,  one  can  eliminate  the  need  for  selecting 
the  parameter  A and  Equation  (71)  is  written  as 
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where  <Prma/cle){  " °* 006* 

Prediction  techniques  of  the  present  study  have  considered  the  viscous  power  law 

exponent  (m)  to  be  0. 8 and  the  velocity  power  law  exponent  (n)  to  be  7.  The  latter  value 

83 

assumes  fully  developed  turbulent  boundary  layer  flow.  However,  Laganelli  et  al 
have  shown  that  fully  developed  turbulent  conditions  (n  = 7)  could  not  be  attained  even 
on  models  five  feet  long.  More  realistically,  the  value  of  n,  as  determined  from  pro- 
file measurements,  was  approximately  9 and  this  parameter  reached  values  as  high  as 
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16  lust  after  transition.  If  a value  of  n = 9 and  m = 0. 8 were  used,  the  ratio  p Ar 

rms  e 

would  reduce  by  appro.. imately  15%.  If  higher  values  of  n were  chosen,  in  particular 
for  the  conical  frustum  region  where  gages  A-19/22/21/22  were  located  just  after 
transition, the  ratio  Prm^ could  reduce  by  almost  25%  and  approach  the  Lowson  pre- 
diction. 

To  include  both  angle  of  attack  and  flap  deflection  effects  on  the  basis  of  limited 
data  acquired  in  the  present  program,  Equation  (90)  might  be  expressed 
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Further  acquisition  and  evaluation  of  appropriate  data  are  required  to  confirm  the 
suggested  functional  form  of  the  first  term  on  the  right  hanc  side  of  Equation  (91). 

With  the  exception  of  data  reported  in  Reference  1 there  is  no  detailed  informa- 
tion available  for  developing  rms  pressure  characteristics  in  transitional  flow  regions. 
In  the  present  tests,  transitional  data  obtained  along  the  conical  frustum  are  shown 
in  Figures  57  and  79.  In  Figure  57,  the  measurements  of  Reference  1 are  shown  to- 
gether with  a recommended  data  fit.  The  algorithm  developed  empirically  from  the 
previous  study  is  expressed  as 

prms/qe  “ *0041A1  + .013  ] ; 3.7  sM^ss.l 
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As  previously  indicated,  Laganelli  et  al  examined  the  effects  of  the  power  law 

velocity  exponent  (n)  in  regions  of  laminar,  transitional  and  turbulent  flow.  It  was 

determined  that  the  power  law  exponent  was  in  the  range  2 <n<  6 during  transition 

flow.  These  results,  which  were  obtained  in  the  AEOC  Tunnel  B facility,  are  con- 
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slstent  with  the  works  of  other  investigators.  Considering  an  average  value  of  n=4, 
together  with  the  viscosity  exponent  m = 0. 8,  Equation  (90)  can  be  written  as 
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Inasmuch  as  rms  pressure  levels  are  not  available  for  incompressible  transitional 
flow,  the  turbulent  value  will  be  used.  Equation  (92)  is  shown  in  Figure  57  with 
<P  /q  ),  1 0. 006. 

The  present  prediction  method  is  seen  to  be  in  agreement  with  the  data  correla- 
tion developed  in  Reference  1.  Moreover,  the  measurements  of  Reference  1 together 
with  current  test  results  appear  to  be  adequately  represented  by  the  present  technique. 
It  -t'ould  be  noted  that  the  transitional  data  shown  in  Figure  57  are  considered  peak 
or  near  peak  transition  levels.  Consequently  the  velocity  power  law  exponent  may 
be  closer  to  6 than  the  average  value  of  4 selected.  In  this  situation,  the  predicted 
normalized  rms  pressure  would  be  reduced  by  20%  and  show  excellent  agreement  with 
the  present  data  as  well  as  with  measurements  of  Reference  1,  in  particular  for  a = 0° 
conditions.  For  this  condition  (n  = 6),  Equation  (90)  becomes 
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2.2  Power  Spectral  Density 

Figures  60  through  68  present  non-dimensional  power  spectral  density  measured 
on  the  cone  as  a function  of  the  far  field  Strouhal  number  (y  > 6^)  using  edge  dynamic 
pressure  as  a normalizing  parameter.  In  order  to  examine  othor  means  of  non- 
dimensionalizing  the  PSD,  Figures  82  through  84  were  constructed.  Figure  82  is  a 
replot  of  Figure  60,  at  Mach  4,  where  edge  dynamic  pressure  has  been  replaced  by 
the  local  static  value.  There  is  no  particular  improvement  in  format,  and  the  data 
do  not  separate  into  two  general  classes  of  transition  and  turbulent  conditions  as  was 
noted  with  Figure  60.  The  most  significant  departure  is  observed  for  the  blunt  transi- 
tion data  whose  levels  were  decreased  to  values  less  than  turbulent  boundary  layer 
values.  The  data  of  Figure  82  tend  to  have  the  same  angle  of  attack  characteristics 
as  displayed  in  Figure  60. 

Figure  83  shows  the  PSD  normalized  using  rms  pressure  as  a function  of  far  field 
Strouhal  number.  The  data  shown  are  a replot  of  Figures  60  and  64  for  zero  angle 

of  attack  conditions  at  Re  /ft.  "3.6x10  . As  in  the  previous  case,  the  data  do  not 

00 

generally  separate  into  turbulent/transitional  characteristics.  However,  a coalescence 
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of  the  measurements  is  observed  compared  to  the  cases  when  dynamic  and  static 
pressure  were  used.  A Mach  number  effect  is  also  apparent  for  the  zero  intercept 
value  in  that  Mach  4 levels  are  generally  higher  than  those  for  Mach  8.  Although  not 
shown,  the  data  spread  (one  decade)  also  includes  the  angle  of  attacK  measurements 
of  Figures  60  and  84. 

Figure  84  displays  the  zero  angle  of  attack  conical  frustum  data  from  the  preceding 

figure  replotted  as  a function  of  the  near  field  (y  < Strouhal  number.  The  absoissa 

represents  die  high  frequency  (low  scale)  contribution  of  the  fluctuating  presrure  field. 

Inasmuch  as  the  ordinate  (normalized  PSD)  remains  invariant,  a change  in  scale  along 

the  abscissa  is  the  only  notable  difference  to  Figure  33.  There  is  no  improvement  in 

collapse  of  the  data.  Moreover,  it  is  difficult  to  Interpret  any  consequences  of  the 

high  frequency  contributions  in  these  measurements. 

The  coordinates  shown  in  Figure  83  represent  the  normalized  functional  forms 
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suggested  for  power  spectral  density  by  several  authors.  ’ ’ ’ Where  measure- 

ments are  plotted  in  this  fashion,  normalizing  with  rms  pressure  tends  to  collapse 
the  data  more  uniformly  than  when  dynamic  pressure  is  used.  Moreover,  an  exami- 
nation of  measurements  of  Reference  1 indicated  the  same  effect.  Consequently,  data 
from  several  experiments,  compressible  and  incompressible,  were  plotted  using  the 
parameters  of  Figure  83  and  are  shown  in  Figure  85.  It  should  be  noted  that  the  circular 
frequency  (uj)  was  used  in  the  data  display. 

In  Figure  85,  it  is  apparent  that  measurements  for  any  one  tunnel  experiment 
tended  to  coalesce  the  data  using  the  rms  pressure.  However,  with  the  exception  ot 

the  high  frequency  range  > 0.  5),  when  many  data  sources  are  considered,  a 

© 

significant  scatter  exist  in  the  low  frequency  range  yielding  zero  intercept  values  with 

9 

no  Mach  number  consistency.  Except  for  measurements  of  Raman,  the  zero  intercept 

values  reflect  the  range  of  rms  pressure  ratio  p A , i.  e. , higher  zero  intercept 

rtDB  w 

values  are  consistent  with  the  lower  values  of  p A . 

rms  w 

Considering  Equation  (89)  together  with  Equation  (90),  there  results 

— — - - k *•  0. 6336,  for  li  ■ unity 
o*  6*  rr 
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An  inspection  of  Figure  85  indicates  significant  departures  from  the  above  value  for 
the  zero  intercept  levels.  Again,  one  recognizes  the  use  of  the  incompressible  format 
for  representing  the  spectrum  which  can  be  used  in  a compressible  state  with  a trans- 
formation function.  The  compressibility  parameter,  e^,,  was  found  to  be  the  trans- 
formation function  required  for  rms  pressure.  On  the  basis  of  the  above  observation, 

(i.  e. , wide  scatter  in  PSD  normalized  by  rms  pressure)  together  with  the  arguments 
presented  in  Section  III,  it  appears  evident  that  the  most  appropriate  way  to  represent 
non-dimensionalized  power  spectral  density  is  to  use  the  dynamic  pressure  as  a 
normalizing  parameter. 

W*  en  rms  pressure  is  plotted  as  a function  of  Mach  number  both  the  ordinate 

(p  /q  ) and  the  abscissa  (M  ) represent  dimensionless  groups  that  can  be  used 
rms  g g 

either  for  incompressible  or  compressible  conditions.  However,  no  such  provision 
is  available  when  displaying  PSD  as  a function  of  Strouhal  number.  Here,  the  indepen- 
dent variable  (Strouhal  number)  can  differ  by  over  an  order  of  magnitude  from  compres- 
sible to  incompressible  conditions.  This  is  primarily  a consequence  of  the  characteristic 
velocity  and  is  reflected  in  the  data  when  comparing  incompressible  and  compressible 

measurements,  i.e. , uj6*/u  ) < < u>  6*/u  )..  It  would  be  desirable  to  extend  the  logic 

e c e i 

employed  in  the  development  of  rms  pressure  with  the  compressibility  factor  to  power 
spectral  density. 

Consider  again  Equations  (89)  and  (90),  such  that 
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It  has  been  observed  that  with  Increasing  Mach  number,  the  zero  intercept  PSD  decreases 
(Figure  24  and  Figures  60-68).  Consequently,  the  right  hand  side  must  also  decrease 
with  increasing  Mach  number.  Inasmuch  as  the  term  (0/qe)0  behaves  In  this  fashion 
(Figures  1,  18,  57-59),  continuity  of  the  above  expression  Is  maintained.  The  com- 
pressible state  can  be  transformed  into  the  incompressible  plane  through  a transformation 
function,  a practice  commonly  employed  in  boundary  layer  theory.  For  example,  the 
skin  friction  coefficient  can  be  expressed  as 
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where  F„  is  the  transformation  function.  The  rms  pressure  has  been  transformed  by 
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such  a process  and  is  represented  by  Equation  (90)  i.  e. , (p  /q  ).  = (P  /q  ) /e 

rms  6 l rms  6 c i 

where  F_  = 1/e  . The  above  zero  intercept  level  can  then  be  written  as 
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For  incompressible  flow,  one  has 
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Hence,  comparing  the  above  expressions  allows  for  the  arbitrary  parameter  k to  be 
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equal  to  the  compressibility  parameter  e_  . The  Strouhal  number  can  also  be  trans- 


formed to 
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whore  F = l/e  = l/(it)  . Equation  (87)  becomes 
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where  it  is  noted  that  the  Strouhal  number  maintains  compressible  values.  The  above 
formulation  represents  the  power  spectral  density,  in  which  the  Houbolt  functional 
relation  has  been  generalized  for  application  to  attached  compressible  turbulent  boundary 
layer  flow.  Here,  one  notes  that  both  the  dependent  variable  (PSD)  and  independent 
variable  (Strouhal  number)  have  the  required  transformation  functions.  Equation  (94) 
is  shown  in  Figure  86  for  several  values  of  as  well  as  measured  data  on  the  conical 
frustum  and  slice  for  the  two  Mach  number  conditions  of  the  present  experiment.  The 
data  are  shown  with  projected  levels  of  the  zero  intercept  value.  When  compared  to 
the  present  theory,  the  data  do  not  precisely  match  the  zero  intercept  values;  however, 
the  general  agreement  is  considered  quite  good  recognizing  the  existence  of  experi- 
mental errors.  It  is  quite  interesting  to  note  that  the  decrease  of  normalized  PSD 
with  Strouhal  number  is  predicted  very  well. 
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Figure  87  shows  the  present  prediction  scheme  compared  to  the  data  of  several 
experiments  (see  Figure  24)  including  conical  measurements  of  Figure  86  and  an 
excellent  match  to  data  is  observed.  As  previously  indicated,  the  zero  intercept  value 
of  measured  normalized  PSD  decreases  with  increasing  Mach  number  and  the  theory  is 
shown  to  be  in  agreement  with  this  trend.  However,  the  decrease  of  the  spectrum 
with  increasing  Strouhal  number  displays  different  characteristics  among  the  various 
investiga*ors. 

For  predicting  levels  of  power  spectral  density,  normalization  by  edge  dynamic 
pressure  and  the  Houbolt  functional  form  are  recommended.  The  prediction  method 
is  shown  as  Equation  (94)  where  the  compressibility  parameter  has  been  employed  to 
transform  the  incompressible  theory  into  a compressible  plane.  Of  particular  impor- 
tance is  the  transformation  on  the  independent  variable  (Strouhal  number).  Concerning 
angle  of  attack  effects,  the  disparity  in  normalized  PSD  measurements  and  trends 
which  are  opposite  from  one  Mach  number  condition  to  another  (Figures  60  through 
68),  it  is  recommended  that  Equation  (94)  be  used  directly.  When  angle  of  attack  data 
for  either  the  cone  or  slice  are  represented  m the  coordinates  of  Figures  86  and  87, 
both  the  zero  intercept  levels  and  frequency  dependence  appear  to  be  predicted  within 
the  scatter  of  the  zero  angle  of  attack  data.  If  one  compares  the  slice  normalized  PSD 
values  of  Figures  61  and  62,  it  is  observed  that  angle  of  attack  effects  appear  more 
sensitive  for  the  Mach  8 tunnel  test.  However,  the  compressibility  parameter  accounts 
for  the  variation  in  angle  of  attack;  in  particular  for  the  higher  Mach  number  tests. 

For  the  flap  region,  it  is  recommended  that  Eq ration  (94)  be  used  with  modification 
to  account  for  flap  deflections.  Good  agreement  between  data  and  theory  are  accorded, 
at  zero  angle  of  attack  and  zero  flap  angle,  for  both  the  zero  intercept  value  and  spec- 
trum variation  with  frequency.  Again,  variation  of  measurements  with  angle  of  attack 
is  accounted  for  through  changes  in  the  compressibility  parameter.  Examination  of 
Figures  62  and  67  indicate  a significant  departure  of  angle  of  attack  effects  from  Mach 
4 to  Mach  8 test  conditions.  In  Tunnel  B the  higher  local  Mach  number  (at  6=0*) 
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yields  lower  levels  of  the  zero  intercept  value  which  subsequently  increase  with  in- 
creasing angle  of  attack.  However,  for  the  Mach  4 test  conditions  at  zero  flap  angle, 
the  variation  in  local  Mach  number  with  a was  not  sufficient  to  cause  significant  level 
changes  in  the  zero  intercept  value. 
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When  considering  flap  deflections,  it  is  estimated  that  the  same  type  of  function 

can  be  used  as  that  employed  for  rms  pressure;  namely,  f(l+6_/9  ).  In  particular, 

r c 

the  normalized  power  spectral  density  for  the  flap  region  can  be  approximated  by 
modifying  Equation  (94)  to  be 
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Equation  (95)  predicts  the  normalized  PSD  for  the  flap  region  over  the  range  of  measured 
flap  deflections  (0  £-6^,  £20°),  including  angle  of  attack  effects,  to  within  the  same  de- 
gree of  accuracy  as  determined  on  the  conical  frustum  and  slice  regions. 

Concerning  transition  prediction  capability  for  power  spectral  density,  direct  use 
of  fluid  dynamic  principles  was  not  possible  as  was  the  case  for  rms  pressure.  An 
inspection  of  the  conical  frustum  transition  data  of  Figures  GO  and  64  show  that  transi- 
tion measurements  yield  higher  zero  intercept  values  than  experienced  for  turbulent 
flow.  It  appeared  that  direct  substitution  of  either  Equation  (92)  or  (93)  into  Equation 
(94)  would  give  the  desired  result.  However,  use  of  these  equations  considered  the  in- 
compressible zero  intercept  value  as  the  maximum  level  (this  was  done  for  lack  of 
incompressible  transition  information).  Consequently,  increasing  compressible  effects 
tended  to  decrease  the  normalized  PSD  zero  intercept  value  regardless  of  choice  of 
velocity  power  law  exponent. 

An  examination  of  transition  data  of  Reference  I together  with  the  present  test 
results  indicated  that  the  zero  intercept  level  tends  to  increase  with  increasing  Mach 
number  which  includes  angle  of  attack  effects.  Hence,  to  accommodate  an  increase 
in  normalized  PSD  with  compressibility,  Equation  (88)  will  be  modified.  For  this 
situation,  k is  assumed  to  be  the  inverse  of  the  compressibility  function  such  that 
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which  represents  the  normalized  power  spectral  density  transitional  flow  zero  intercept 
level.  The  choice  of  the  above  transformation  has  no  direct  fluid  dynamic  foundation 
as  previously  employed  for  rms  pressure.  Moreover,  the  spectral  dependence  with 
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frequency  is  difficult  to  assess.  Inasmuch  as  the  data  of  Reference  1 aa  well  as  the 

current  measurements  tend  to  have  steep  rates  of  decrease  for  Strouhal  number 
-2 

f6*/u  > 10  , it  1b  suggested  that  the  spectral  distribution  with  frequency  remain 
© 

unchanged  in  Equation  <94),  thus  giving 
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Equation  (97)  is  shown  plotted  in  Figures  60  and  64  over  the  low  frequency  range 
only.  Also  represented  is  the  range  of  the  compressibility  parameter  for  the  current 
test  conditions.  It  is  noted  that  the  velocity  power  law  exponent  of  n = 4 was  assumed. 
The  data  are  in  good  agreement  with  theory  (Equation  (97))  and  more  importantly  are 
correctly  predicted  relative  to  compressible  effects. 

2.3  Cross-Power  Spectral  Density 
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Detailed  analysis  was  performed  for  those  data  considered  most  applicable  to  the 

flight  application  situation  to  be  discussed  in  Section  II  of  Volume  n.  It  was  desired 

to  have  as  much  of  the  model  as  possible  emersed  in  fully  turbulent  flow  (at  0=0), 

with  primary  emphasis  placed  on  those  runs  wherein  the  flap  deflection  was  consistent 

with  such  a zero  angle  of  attack  condition  (i.e. , -6  = 15°  or  20°).  For  evaluation  of 

r 

conical  frustum  data,  Run  103  at  Mach  4 proved  to  be  the  only  case  where  the  conical 
sensor  array  was  emersed  in  fully  turbulent  flow.  However,  for  measurements  on 
the  slice  and  flap,  Mach  8 Runs  58  (sharp  cone)  and  100  (blunt  nosetip  with  mass  in- 
jection) had  transitional  flow  sufficiently  far  forward  that  adequate  turbulent  data  were 
acquired  over  these  regions.  Results  from  wind  tunnel  runs  other  than  the  above 
three  will  be  introduced  as  ne:essary  to  evaluate  effects  of  flap  deflection  and  angle 
of  attack  variation,  as  well  as  transition.  Cross-spectral  data  corresponding  to  addi- 
tional runs  are  presented  in  Volume  n. 

The  initial  step  in  data  evaluation  consisted  of  compiling  individual  cross-power 
spectral  densities  into  a set  of  functions  corresponding  to  each  applicable  transducer 
array  of  the  wind  tunnel  run  under  investigation.  For  longitudinal  correlations,  peaks 
of  the  normalized  co-spectral  density  (C)  and  the  quad-spectral  density  (Q)  are  ex- 
pressed as  functions  of  non-dimensionalized  frequency,  ^f/uQ(0,  by  plotting  C 6 Q 
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versus  G/2tt  since  by  definition  u (f)  ~ 2rrf§/0  (f).  It  should  be  noted  that  for 

c 

§f/u  (f)  - 0 asymptotic  value  of  normalized  co-spectral  density  decreases  for  in- 
c 

creasing  values  of  transducer  separation  distance.  Hence,  the  effect  of  both  spatially 

dependent  and  frequency  depevident  cross-spectral  attenuation  is  evident. 

In  the  case  of  circumferential  correlations,  only  the  in-phase  (co)  cross -spectral 

data  are  relevant.  These  functions  are  expressed  in  terms  of  the  quantity  nf/uc(f), 

explicit  values  of  which  are  generated  as  follows:  the  longitudinal  convection  velocity, 

u (f),  is  defined  as  above  by  2nf!L/0,(f)  where  §,  is  a specific  longitudinal  separation 
c 111 

distance  and  0 (f)  is  the  corresponding  phase  angle.  Therefore,  7)  f/u  (f)  = 7]  f r 

1 1C  1 

2irf?1/8(f)-T)1/51xejL(l)/aTT.  As  an  example,  consider  the  circumferential  co- 

spectral  density  associated  with  transducers  5-6.  For  this  case 

^ 5-6  95-8  (f)  . 75  e5-8(f)  1.2e5-8 

T'5-6f/uC,Long  W ■ §5  ' 2n  “.625  2n  “ 3601 

In  effect,  then,  to  provide  the  appropriate  representation  as  functions  of  T)f/u  , cir- 

c 

cumferential  co-power  spectral  densities  are  plotted  against  scaled  values  of  longi- 
tudinal phase  angle. 

The  approach  taken  in  the  present  investigation  was  to  consider  each  of  the  three 
geometric  regions  of  the  model  in  turn  (cone,  slice,  flap)  and  relate  cross-spectral 
characteristics  to  correlation  properties  measured  in  the  study  of  Reference  7.  In  that 
program,  it  was  proposed  that  appropriate  representations  for  the  turbulent  cross- 
spectral  coefficients  and  over  a cone  are  Equations  (37)  and  (38),  respectively. 

2.3.1  Conical  Frustum  Measurements 

Figures  12,  14,  15,  and  88  present  applicable  cross-spectral  data  for  evaluation 
of  correlation  characteristics  over  the  cone.  Measurements  expressed  both  in  terms 
of  non-dimensional  frequency  and  asymptotic  values  (as  f - 0)  versus  normalized 
separation  distance  are  shown. 

(a)  Longitudinal  Correlation 

(1)  Frequency  Dependence 

Equation  (37)  indicated  that  even  if  f - 0,  such  that  f§/ue  - 0,  as  long  as  l ft  0 
(l.e. , for  any  actual  cross-power  spectral  density  measurement),  the  asymptotic  value 
for  fS/uc  - 0 will  be  less  than  1.  Therefore,  measured  cross -speotral  data  are  shown 
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compared  to  an  exponential  envelope  which  has  a factor  in  it  to  account  for  the  spatial 

attenuation.  For  this  reason,  in  Figure  88  a factor  of  . 7 has  been  Incorporated  to 

acoount  for  attenuation  corresponding  to  a separation  distance  of  approximately  1. 5". 

It  was  noted  in  Reference  7 that  the  effect  of  noise  contribution  (i.  e. , tunnel  tare)  on 

cross-PSD  measurements  is  apparently  to  inhibit  decay  of  the  cross- spectral  envelope. 

Hence,  in  this  figure,  the  relatively  constant  correlation  indicated  for  sensors  A-5 

and  A-10  (§/6*  = 70. 2)  in  the  range  f?/u  (f)  > 1.  5 is  probably  an  effect  of  tunnel  tare. 

c 

The  comparison  of  available  turbulent  data  shown  in  Figure  88  with  the  function 
e ' ^TOfS/Ug)  reveai8  this  representation  does  provide  a reasonable  envelope  of 
measured  values  when  factors  to  account  for  spatial  attenuation  are  implicitly  incor- 
porated. 

(2)  Spatial  Dependence 

Figure  12  presents  longitudinal  data  for  Run  103  (Mach  4)  in  which  asymptotic 
values  from  Figure  88  (as  f - 0)  are  plotted  versus  %/ 6*.  The  levels  of  present  data 
are  slightly  higher  than  both  the  primary  spread  of  measurements  from  Reference  7 
and  the  prediction  curve  derived  therein.  Redefinition  of  the  mathematical  represen- 
tation would  not  appear  to  be  justified,  solely  on  the  basis  of  the  present  (limited)  data. 
However,  when  the  new  representation  shown  in  Figure  12  (.  58  + .42  e ’ ^ ) is 

incorporated  into  the  R/V  response  routine,  it  provides  a better  match  to  measured 
re-entry  vibration  data  (Volume  II).  Until  this  ambiguity  is  resolved  through  further 
investigation,  it  is  recommended  that  the  new  equation  be  employed  in  design  analyses 
on  the  grounds  that  it  constitutes  a slightly  more  conservative  representation. 

(b)  Circumferential  Correlation 

(1)  Frequency  Dependence 

Figure  14  presents  a comparison  of  circumferential  co  cross-power  spectral 

7 

density  from  Run  103  (Mach  4)  with  data  of  previous  compressible  measurements. 

This  figure  includes  from  current  data  only  the  correlation  data  for  sensors  A-5  and 
A-6  (7)  - . 75"),  since  sensor  A-7  did  not  provide  valid  measurements  in  either  Tunnel 
A or  B.  Also  shown  in  Figure  14  is  the  exponential  recommended  in  Reference  7 to 
represent  circumferential  correlation  (e  ^(^TTjfn/u^  multiplied  by  .85  to  acoount 
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for  average  spatial  attenuation  of  the  data.  There  is  generally  good  agreement  among 
all  the  measurements  and  the  empirically  developed  function;  hence  it  is  reoommended 
this  exponential  be  maintained  as  an  appropriate  representation  for  decay  of  circum- 
ferential correlation  with  non-dimensional  frequency. 

(2)  Spatial  Dependence 

A comparison  of  asymptotic  measurements  for  circumferential  cross-power 
spectral  density  is  provided  in  Figure  15.  As  with  the  longitudinal  correlation,  the 
present  data  (one  point)  are  above  the  previous  prediction  curve,  and  provide  an  ex- 
ponential expression  which  ultimately  results  in  better  agreement  with  measured  flight 
vibration  data  when  incorporated  into  the  R/V  response  routine.  Therefore  it  is 
recommended  that  e”°’  019j  be  employed  in  design  analyses  until  further  data 
are  acquired. 

2.3.2  Slice  and  Flap  Measurements 

From  an  initial  evaluation  of  correlation  characteristics  of  both  slice  and  flap  data, 
it  became  evident  that  a general  similarity  exists  among  these  measurements.  This 
is  not  surprising  considering  the  planar  nature  of  these  two  geometries.  Therefore, 
data  from  the  slice  and  the  flap  will  be  considered  together.  It  should  be  noted  that 
certain  anomalous  characteristics  were  evident  in  the  output  of  sensor  A-28  (e.  g. , 
the  spectrum  shown  in  Figure  56).  Therefore,  when  available,  the  correlation  relating 
sensors  A-23  and  A-25  (equivalent  flap  data)  was  substituted  for  cross-spectrum  A-26 
to  A-28. 

(a)  Longitudinal  Correlation 
(1)  Frequency  Dependence 

It  is  evident  from  an  examination  of  Figures  89  through  93  that  slice  and  flap 
cross-spectral  measurements  in  turbulent  flow  exhibit  a more  significant  loss  of 
correlation  with  non-dimensional  frequency  than  do  conical  frustum  data  discussed 
above.  Hence,  a modification  to  the  previous  prediction  curve  is  required.  One 
possible  exponential  representation  which  matches  available  data  reasonably  well  is 
e • 2<2r')^  /uo  K tB  apparel  that  the  rate  of  correlation  deoay  (i.  e. , the  coefficient 
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of  fS/u  ) In  cross-spectral  measurements  for  the  slice  and  the  flap  is  generally  in- 
c 

sensitive  to  flap  deflection.  Therefore  the  above  exponential  is  applicable  for  zero 
angle  of  attack  conditions  and  for  flap  deflection  in  the  range  0 £ -6_  £ 20°. 

r 

(2)  Spatial  Dependence 


Asymptotic  values  (as  f -•  0)  of  longitudinal  cross-PSD  on  the  slice  and  on  the 
flap  are  presented  in  Figure  94.  It  is  evident  that  for  small  values  of  normalized 
separation  distance  (i.e. , 5/4*  £ 7)  data  are  generally  coincident  with  the  prediction 
curve  which  provided  a good  match  to  conical  frustum  data  from  Reference  7.  How- 
ever, at  larger  values  of  5/6*  the  degree  of  asymptotic  correlation  decreases  relative 

to  this  curve.  This  trend  of  both  slice  and  flap  measurements  (0  £ -6  £ 20°)  mandates 

r 

a revision  to  the  conical  frustum  prediction  curve.  A possible  modification  which  provides 
better  agreement  with  present  data  is  0. 35  + 0. 65  e * . This  expression  is 

applicable  to  the  zero  angle  of  attack  situation  when  turbulent  boundary  layer  flow 
emerses  the  slice  and  flap  region.  Variation  with  Mach  number  and  flap  deflection 
is  automatically  accounted  for  in  the  functional  dependency  of  6*  on  these  parameters. 

(b)  Lateral  Correlation 


Due  to  the  erratic  behavior  of  sensors  A-16,  17,  18,  no  valid  lateral  cross-spectral 
data  over  the  slice  were  acquired.  From  a review  of  lateral/circumferential  measure- 
ments both  in  the  present  investigation  (cone)  and  throughout  the  literature  (cones, 
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plates,  wind  tunnel  walls),  ’ ’ ’ it  1b  apparent  that  the  general  level  of  lateral 

cross-correlation  is  much  less  sensitive  to  the  geometric  configuration  than  longitudinal 
correlation  (see  Figures  13  and  14).  Therefore  the  functional  form  for  lateral  cross- 
spectral  coefficient  proposed  in  this  study  is  identical  to  the  expression  for  circum- 
ferential correlation  discussed  in  Section  2. 3. 1 above. 

In  summary,  the  complete  representations  for  cross-spectral  coefficients  de- 
scribing correlation  in  fully  turbulent  flow  over  the  slice  and  flap  (0  £ -6„  £ 20s)  are 

r 

as  follows: 


A§(S..d)  = e--2u,5/Uo(u,)[.35  + .65  e’*073  ^ 


] 


(98) 


A^Ol,®) 


e-1.4«BTl/u0(tt,)[  3+  7e-.055V6*3> 


(99) 
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2. 3. 3 Transition  and  Angle  of  Attack  Effects 


(a)  Conical  Frustum 

Figure  95  presents  Mach  4 transitional  cross-spectra  from  the  cone  array  com- 
pared to  the  frequency-dependent  prediction  curve  proposed  in  Reference  7.  Good 
agreement  with  the  proposed  exponential  is  evident.  An  examination  of  other  oonloal 
data  reveals  that  the  location  of  the  transition  zone  with  respect  to  the  sensor  array 
can  influence  measured  cross-correlation  characteristics.  However,  no  definitive  trends 
in  this  respect  were  exhibited  throughout  the  data.  It  is  therefore  recommended  that 

7 

the  exponential  developed  by  Howe  on  the  basis  of  many  variations  in  angle  of  attack 
and  freestream  Reynolds  number  be  maintained. 

A comparison  of  asymptotic  longitudinal  data  (as  f -*  0)  from  the  present  program 
with  cross-spectral  measurements  from  Reference  7 indicates  generally  good  agree- 
ment, as  shown  in  Figure  96.  Model  angles  of  attack  up  to  14°  are  represented.  The 
prediction  curve  apparently  constitutes  an  upper  bound  to  current  data,  particularly 
at  larger  values  of  normalized  separation  distance.  The  greatest  discrepancy  occurs 
for  correlations  between  A-5  and  A-10.  However,  an  examination  of  OASPL  distributions 
shown  in  Section  IV  reveals  that  A-10  generally  registers  lower  values  than  might  be 
expected  (e.g. , Figure  34).  These  possibly  erroneous  readings  may  be  contributing 
to  the  low  longitudinal  correlations  in  Figure  96  at  the  largest  values  of  £/6*. 

Typical  high  angle  of  attack  (a  = 7°)  circumferential  data  measured  in  transitional 
flow  over  the  cone  are  given  in  Figure  97.  Also  presented  are  both  data  and  the 
empirically  developed  curve  from  Reference  7.  Current  data  reveal  approximately 
the  same  trend  with  normalized  frequency  as  previous  transitional  measurements 
at  high  model  angle  of  attack. 

On  the  basis  of  the  above  comparisons,  it  is  evident  that  conical  frustum  measure- 
ments from  the  present  program  exhibit  the  same  overall  characteristics  as  transi- 
tional cross-PSD  data  of  Reference  7.  It  is  therefore  recommended  that  the  empirical 
representations  derived  therein  be  maintained  for  cross-spectral  coefficients  in  transi- 
tional flow  over  conical  portions  of  a re-entry  vehicle,  and  that  the  same  angle  of  attack 
criterion  be  applied  for  measurements  at  high  a. 
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These  representations  are  as  follows: 


A (l.a))-e~*06u,5/uo(uj)[.6  + .4  e'*024  §/6*] 


(100) 


A,  CH.«)  - e-*68^/uc(-)[<6+  >4  e-.024T)/6*L  ^ 


£.7  9 


(100) 


A (T),<ju) 


-.446uuT|/u0(uu) 


[. 6 + .4  e 


024  T)/ 6* 


]»  a >.7  0 


(102) 


(b)  Slice  and  Flap 

Transitional  flow  over  the  control  surface  region  of  the  model  constitutes  an  ex- 
tremely complex  phenomenon.  An  example  of  this  is  exhibited  in  Figure  29  which 
presents  OASPL  data  for  two  Mach  4 runs  (69,  84)  wherein  transitional  flow  exists 
over  the  aft  portion  of  the  model.  Although  heat  transfer  measurements  apparently 
Indicate  nominally  turbulent  flow  on  the  slice  and  flap,  acoustic  data  (in  particular 
sensors  A- 17  and  A-26)  reveal  much  higher  than  anticipated  fluctuating  pressure  environ- 
ments. Measured  levels  on  the  slice  and  flap  may  be  a consequence  of  a type  of  "over- 
shoot" phenomenon  for  the  control  surface  region,  analogous  to  the  characteristic  of 
acoustic  environments  near  the  end  of  transition  on  a cone.  The  result  of  these  ob- 
servations is  that  limited  cross -spectral  measurements  corresponding  to  "transitional 
flow"  on  the  control  surface  region  are  inconclusive,  since  aerodynamic  characteristics 
change  so  rapidly  in  that  particular  region  of  the  scale  model. 

It  was  noted  in  Section  2.3.2  above  that  in  fully  turbulent  flow,  flap  longitudinal 
cross-spectral  properties  tend  to  be  independent  of  flap  deflection.  This  same  ob- 
servation can  be  made  with  respect  to  model  angle  of  attack.  Therefore,  Equations 
(98)  and  (99)  would  also  be  appropriate  for  0°  * a ± 14° . However,  since  in  practical 
applications  maximum  flap  deflection  generally  occurs  for  minimum  angle  of  attack, 
this  limits  the  range  of  & , 6 ^ combinations  for  which  Equations  (98)  and  (99)  are 
applicable  (i.e. , a = 14® , -6p  = 20®  would  be  excluded  from  these  equations). 
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Figures  98  and  99  present  slice  cross-PSD’s  In  turbulent  flow  for  two  angles  of 
attack  (0° , 14*).  The  asymptotic  values  (as  f - 0)  are  virtually  Identical,  and  both 
sets  of  data  agree  reasonably  well  with  the  exponential  recommended  above  for  slice 
cross-spectral  coefficient  at  a = 0.  Although  the  higher  angle  of  attack  measurements 
appear  to  decay  slightly  more  rapidly  with  non-dimensional  frequency  than  do  the  zero 
angle  of  attack  cross -spectra,  the  trend  appears  too  subtle  to  warrant  re-definltion 
of  the  correlation  coefficient  on  the  basis  of  these  limited  data.  Hence  it  is  similarly 
recommended  that  Equations  (98)  and  (99)  apply  to  slloe  measurements,  Independent 
of  angle  of  attack,  in  the  range  0“  s a s 14* . 
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Figure  13.  Comparison  of  Narrow  Bsnd  Lateral  Space  Correlation 
Coefficients  from  the  literature 
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Figure  21.  First  Moment  of  Power  Spectral  Density  with  Far  Field  (y  > 6^  Strouhal  Number  (a  = 
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Figure  25.  Surface  Sensor  Location  - Identification 
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Figure  29.  Acoustic  and  Hiermal  Distribution  (Mach  4)  - dp  Effect  at  Re  /ft  - 2.5  x 10 
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Figure  30.  Acoustic  and  Thermal  Distribution  (Mach  4)  - Re-  Effect  at  or  = 0 
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Figure  47  Comparison  of  Conical  Frustum  Spectra  in  Laminar  Flow  for  Various 
Conditions  - Mach  8 
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Figure  50  Comparison  of  Turbulent  Boundary  Layer  Speotra  Along  Modal 
8urfaoe  - Sharp  Body  - Maoh  8 
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Figure  55  Comparlaon  of  Fi«p  Spectra  in  Turbulent  Flow  at  Various  Value*  of 
Flap  Deflection  - Mach  8 
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Figure  56  Comparison  of  Flap  Spectra  in  Fully  Turbulent  Flow  for  Sensors 
A-26,  -27,  and  -28  - Mach  8 
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Figure  5".  Normalized  RMS  Pressure  Measurements  Compared  to  Prediction  - Conical  Frustum 
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Figure  58.  Normalized  RMS  Pressure  Measurements  Compared  to  Prediction  - Slice  Region 
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Figure  59.  Normalized  RMS  Pressure  Measurements  Compared  to  Prediction  - Flap  Region 
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Figure  60.  Normalized  Power  Spectral  Density  for  Turbulent/Transition  Flow 

on  Conical  Frustum  - Mach  4 
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Figure  65.  Normalized  Power  Spectral  Density  for  Tuibulent  Flow 
on  Slice  Region  (Mach  8)  - RN  and  Spatial  Effects 


0 <f)  Ue/q2e  6 


<fl  iyq2e  6 


102 


* 


9- 


103 


io-4 


10‘5 


io-® 


10'7 


- Moo' 

■ FLA 

: v 

o 0 ° o 
° ° _ 

o • ° ° 

■8  R#  /FT 

P REGION  TBL  » “ 

*B“0  SYM  RUN/SEN  a dF  x 10® 

A 42/27  7 0 3.6 

□ 44/27  0 0 3.6 

7 46/27  14  0 3.6 

o 58/27  0 20  3.6 

O 77/27  0 20  2.5 

— o _ 

_ o o o ° 

o o o o o 

0 O 
o o 

► o 

rat 

7 7 7 7 V V ^ 

>7 

A 

A A A A ^ 

7 

it  A £ V 

^ A.  7 

= A A A 

Q 

D a o 

fa  □ 

° D □ a 

& y 

4 <t 

V 

° a 0 o 

□ 

□ 

0 

« 

a 

..  < 

1 LI  Mill. 

0 a D 

L - 1 1 1 i 1.  J 1 

1 1 1 . LI.1.J1  . 

io-3 


10 


•2 


f6*/U# 


10'1 


10° 


Figure  67.  Normalized  Power  Spectral  Density  for  Turbulent  Flow 
on  Flap  Region  (Mach  8)  - a and  6p  Effects 
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Figure  68.  Normalized  Power  Spectral  Density  for  Turbulent  Flow 
on  Flap  Region  (Mach  8)  - R Effects 
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Figure  70.  Broad-Band  Space-Time  Correlation  Functiona  (Longitudinal)  Measured  on 

the  Conical  Fmatum  - Mach  8 
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Figure  72.  Broad-Band  Space-Time  Cross -Correlation  Functions  (Longitudinal)  Measured 

on  Slice  - Mach  4 


188 


'o  '3)a  'o  '3)a  v 'o  '5)a 


(ANALYSES  TO  BE  PERFORMED  ON  SELECTED  SENSORS*  TESTSI 


Figure  77.  Data  Processing  Analysis  Diagram 
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Mach  4 and  8 Test  Measurements 
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Figurs  79.  Normalied  HMS  l re 3 sure  (with  Static  Pressure)  Variation  with 
Ma«*h  Number  - Conical  Frustum  Measurements 
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Figure  80.  Normalised  RMS  Pressure  (with  Static  Pressure)  Variation  with 
Mach  Number  - Slice  Region  Measurements 
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Figure  81.  Normalized  RMS  Pressure  (with  Static  Pressure)  Variation  with 
Mach  Number  - Flap  Region  Measurements 
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Figure  82.  Normalised  Power  Spectral  Density  (with  Static  Pressure)  - Mach  4 

Conical  Frustum  Measurements 
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Figure  83.  Normalized  Power  Spectral  Density  (with  RMS  Pressure)  Variation 
with  High  Scale  Strouhal  Number  - Conical  Frustum  Measurements 


Figure  85.  Comparison  of  Normalized  Power  Spectral  Density  fc.-  Incompressible  and  Compressible  Conditions 
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Figure  87.  Comparison  of  Present  Normalized  Spectral  Measurements  with  Other  Spectral  Data  and  Theory 
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Figure  88.  Longitudinal  Cross-Power  Spectral  Density  Measured  in  Fully  Developed  Turbulent  Flow 

Conical  Frustum  (M  ^ -n 


4.  GROUP  = 32,  TEST  = 103,  a = 0.0,  5 


Figure  89.  Longitudinal  Cross-Power  Spectral  Density  Measured  in  Fully  Turbulent  Flow  - 

Slice  (M  = 4) 
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Figure  90.  Longitudinal  Cross-Power  Spectral  Density  Measured  in  Fully  Turbulent  Flow  - 

Slice  (M  = 8) 


Figure  91.  Longitudinal  Cross-Power  Spectral  Density  Measured  in  Fully  Turbulent  Flow 
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Figure  94.  Asymptotic  Values  of  Longitudinal  Cross -PSD  at  /uc  (U  - 0 Measured  on  the 

Slice  and  Flap  (M  = 4,  8) 
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f?/Uc  <f> 

Figure  95.  Longitudinal  Cross-Power  Spectral  Density  Measured  in  Transitional  Flow 

Conical  Frustum  (M  - 4) 
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Figure  96.  Asymptotic  Values  of  Longitudinal  Cross-Power  Spectral  Density  at  f?  / Uc  (f)  = 

Transitional  Flow  (M  = 4,  8,  10) 
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Figure  97.  Circumferential  Cross-Power  Spectral  Density  - Transitional  Flow  pt 

High  Angle  of  Attack  (M  = 4,  8) 
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Figure  98.  Longitudinal  Cross-Power  Spectral  Density  Measured  in  Fully  Turb'ilent  Fiow 

Slice,  Sharp  Cone  Model  (M  = 4) 
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APPENDIX:  APPLICATION  OF  VIBRATION  PREDICTION  METHODS  TO  A 
BALLISTIC  RE-ENTRY  VEHICLE 

1.0  INTRODUCTION 

Computation  oi  R/V  response  to  boundary  layer  fluctuating  pressures  is  performed 
by  extending  the  principles  developed  in  a previous  GE/AFFDL  contract  (Reference 
99).  Response  of  both  a ballistic  and  a maneuvering  design  are  investigated  in  Volume 
n,  utilizing  the  overall  approach  indicated  in  Figure  A-l.  This  appendix  summarizes 
the  unclassified  techniques/results  of  the  ballistic  R/V  analysis  which  incorporates 
key  parameters  of  the  fluctuating  pressure  environment  derived  in  Volume  I for  the 
conical  portion  of  the  model.  These  parameters  are  the  overall  rms  pressure,  power 
spectral  density,  and  cross-power  spectral  density  (narrow-band  spatial  correlation 
function).  Ballistic  re-entry  flight  data  and  corresponding  vibration  predictions  at  one 
altitude  are  compared  as  a baseline  which  establishes  overall  validity  of  the  upgraded 
computational  technique.  , 

The  particular  ballistic  configuration  investigated  consists  of  a prototype  opera- 
tional design.  Responses  at  various  internal  locations  (see  Figure  A-2)  were  monitored 
for  several  flights  having  approximately  the  same  trajectory  and  freestream  aero- 
dynamic characteristics.  A typical  vibration  response  power  spectral  density  compila- 
tion for  the  axial  direction  is  presented  in  Figure  A-3.  For  the  ESD,  lateral  flight 
data  are  generally  higher  than  those  vibration  levels  which  would  be  applicable  in  com- 
ponent ground  testing,  due  to  the  fact  that  the  component  extremity,  rather  than  the 
CG  or  the  individual  mounting  location,  was  instrumented.  Maximum  vibration  levels 

corresponding  to  fully  developed  turbulent  boundary  layer  flow  were  measured  to  be 

2 2 

approximately  the  following:  . 007  g /Hz  axially  on  the  payload,  . 02  g /Hz  laterally 

2 

on  the  aft  bulkhead,  and  .3  g /Hz  laterally  at  the  tip  of  the  ESD. 

2.0  ANALYSIS 

In  the  ballistic  vehicle  dynamic  model,  the  R/V  structure  is  idealized  as  an  axi- 
symmetrlc  lumped-parameter  (spring-mass)  system.  As  a consequence  of  basic 
vehicle  symmetry,  axial  motions  and  bending  motions  (lateral,  rotational)  are  mutually 
uncoupled  and  are  therefore  treated  separately;  further,  bending  motions  are  simplified 
by  considering  planar  motion  only.  Dynamic  motions  are  described  by  totals  of  46  axial 
and  80  lateral- rotation  degrees -of- freedom  at  vehicle  mass  points  (model  nodes). 
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The  structural  dynamic  model  is  comprised  of  the  following  subsystems:  the 
nosetip-forward  frustum,  mid  frustum,  aft  frustum,  aft  cover,  payload,  bulkhead  with 
components,  and  the  forward  component.  A beam  finite -element  formulation  is  the 
basis  for  the  frustum  subsystems,  which  together  contain  a total  of  27  nodes  (26 
elements)  represeuting  the  vehicle  external  structure  (shell)  and  include  the  effects  of 
the  two  field  joints.  The  remaining  subsystems  are  basically  semi-empirlcal  in  that  they 
were  originally  defined  by  analysis  and  subsequently  updated  on  the  basis  of  correlations 
with  ground  test  results.  The  aft  cover  subsystem  comprises  seven  model  nodes,  and 
couples  to  the  aftmost  vehicle  shell  node;  in  the  bending  model,  this  subsystem  is 
simplified  to  a mass  effect  only  at  this  shell  location  because  of  the  high  frequency 
nature  of  the  aft  cover  lateral  characteristics.  Ten  model  nodes  describe  the  motions 
of  the  payload  and  its  support  structure,  and  this  subsystem,  defined  by  the  payload 
contractor,  is  supported  at  three  points  on  the  vehicle  shell.  The  component  bulkhead 
model  subsystem  consists  of  seven  nodes,  representing  the  shell  support  point,  bulk- 
head center  and  five  principal  components  individually.  Finally,  the  forward  component 
subsystem  contains  five  nodes,  four  representing  vehicle  shell  support  locations  and 
one  the  motion  at  the  assembly  center. 

Total  model  stiffness  and  mass  matrices  were  developed  by  assembling  the  para- 
meters of  the  model  subsystems.  Then  the  eigenvalue  problem  posed  by  these  matrices 
in  the  equations  of  motion  was  solved  by  computer  to  give  the  R/V  natural  resonant 
frequencies  and  mode  shapes.  These  dynamic  characteristics,  summarized  in  Table 
A-I,  were  then  recorded  on  computer  tape  for  subsequent  use  in  response  analyses. 

Analysis  of  ballistic  vehicle  response  concentrated  on  the  point  in  the  trajectory 

where  dynamic  pressure  is  a maximum.  Boundary  layer  characteristic  parameters 

corresponding  to  this  particular  attitude  were  calculated  along  the  vehicle  surface  using 

93 

GE-RESD  computer  program  3DFF.  Applicable  boundary  layer  properties  so  computed 
are  summarized  for  R/V  aft  cone  locations  in  Table  A-n. 

Equations  utilized  to  convert  aerodynamic  characteristics  to  the  form  required 
for  computation  of  R/V  structural  response  UBing  the  approach  of  Figure  A-l  are 
described  below. 
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2.1  PSD 


Equation  (94)  of  Volume  I yields  the  following  expression  for  power  spectral 
density: 

cp(f)  =2tt<p(u>)=  <g/(*e>i  eT/ue  (Al) 

1 + e"4  (2rrf6*/u  )2 
X 0 


where 


(n\\ 


= .006 


and 

«,  = fT*/T  m 

1 0 

with 

T*/Te  = [1/2  (1  + Tw/Te)  + . 22  r (^)  ] (A3) 

For  fully  turbulent  conditions  during  re-entry,  in  which  the  R/V  wall  temperature  is 
very  much  less  than  the  adiabatic  wall  temperature,  (A2)  reduces  to 

2 - 

c = [0.5+  0.02  M~1  (A4) 

X 6 

with  r = turbulent  recovery  factor  = . 896, 

V - 1.2  for  re-entry, 
n = 7, 

and  m = .6  (from  Reference  85). 

Therefore  Equations  (Al)  and  (A4)  define  the  PSD  over  a ballistic  re-entry  vehicle 
subjected  to  fully  turbulent  flow. 

2.2  Cross-PSP 

It  was  discovered  that  the  longitudinal  cross-spectral  coefficient  representation 
which  was  generated  on  the  basis  of  extremely  limited  TBL  data  over  the  cone  (Run  103, 
Tunnel  A)  provided  a better  match  with  ballistic  flight  data.  Therefore  this  cross-PSD 
representation  is  presented  here.  Equations  utilized  in  the  computation  of  ballistic 
re  -entry  vibration  response  are: 
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(A  5) 


A?(?,a')  =e"0,1U^/^[0.58  + 0.42  e'0, 0486 ?/6*] 

= e_1‘4u,T'/Uc  e-0*  019511/6* 

3.  0 RE-ENTRY  VIBRATION  PREDICTIONS 

Table  A-in  summarizes  computations  for  the  ballistic  vehicle  and  presents  maxi- 
mum response  levels  for  various  internal  and  external  (nose,  shell)  locations.  It  is 
evident  that  there  is  quite  good  correlation  with  flight  measurements  derived  from 
envelopes  of  response  spectra  in  fully  turbulent  flow,  except  for  bulkhead  axial  responses. 
Here,  one  notes  an  order  of  magnitude  discrepancy  between  maximum  flight  data  and 

the  predictions  (Figure  A-3).  However,  it  should  be  pointed  out  that  the  particular 

2 . 

response  PSD  which  contributes  the  . 022  g /Hz  peak  is  approximately  one  decade  higher 
than  the  majority  of  power  spectra  obtained  in  other  similar  flights  for  the  same  in- 
strumented location.  Therefore  this  bulkhead  axial  measurement  represents  a 3a  or 
higher  response,  while  predictions  based  on  nominal  conditions  provide  generally  good 
agreement  with  other  data  shown. 

For  bulkhead  lateral  responses,  flight  data  are  respectively  slightly  less  than  or 
greater  than  predicted  levels,  depending  on  location.  That  is,  for  locations  designated 
NG  and  FS,  flight  measurement  stations  are  oloser  to  the  aft  bulkhead  than  the  model 
node  points,  such  that  the  moment  arms  which  affect  lateral  response  are  shorter. 

This  situation  is  reversed  for  the  ESD  where  measurements  were  acquired  on  the  tip 
of  the  component;  hence  at  lateral  rocking  frequencies,  flight  responses  are  magnified 
over  levels  predicted  at  the  ESD  center  of  gravity. 

Maximum  vibration  levels  computed  for  the  R/V  nose  and  for  the  shell  are  also 
presented  in  Table  A-m.  Axial  and  lateral  responses  correspond  respectively  to  the 
0th  and  1st  harmonios  of  the  vehicle  structural  dynamic  model.  If  effects  of  higher 
harmonics  were  accounted  for,  significantly  higher  shell  responses  would  result,  but 
component  levels  would  be  virtually  unchanged.  This  is  due  to  the  faot  that  component 
packages  are  generally  mounted  on  support  structure  such  that  their  primary  behavior 
is  simple  beam-column  motion.  The  support  structure  does  not  transmit  vibration 
levels  associated  with  high  harmonic  shell  response. 
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TABLE  A-I 

BALLISTIC  R/V  DYNAMIC  MODEL; 
SUMMARY  OF  DYNAMIC  CHARACTERISTICS 
TO**  2000  HZ 


1 

AXIAL 

LATERAL 

ORDER 

NAT.  FREQ. 
(HZ) 

MODE  SHAPE 
DESC. 

NAT.  FREQ 
(HZ) 

MODE  SHAPE 
DESC.* 

1 

0 

RIGID  BODY 

0 

RIGID  BODY 

2 

110.4 

AFT  COVER 

0 

RIGID  BODY 

3 

159.8 

AFT  BULKHEAD 

184.6 

SHELL 

4 

216.2 

FORWARD  COMP 

213.5 

SHELL-P/L-AFT  B/H 

5 

334.9 

AFT  BULKHEAD 

236.4 

AFT  BULKHEAD 

6 

459.7 

SHELL- PAYLOAD 

272.2 

PAYLOAD 

7 

531.8 

AFT  BULKHEAD 

378.1 

SHELL-P/L 

8 

568.9 

AFT  COVER 

410.7 

AFT  BULKHEAD 

9 

841.2 

AFT  COVER 

423.0 

AFT  BULKHEAD 

10 

875.7 

AFT  BULKHEAD 

459.9 

AFT  B/H-SHELL 

11 

921.0 

R/V  AXIAL 

494.7 

SHELL-P/L 

12 

966.8 

AFT  BULKHEAD 

564.7 

SHELL-F/C-AFT  1/11 

13 

1089.2 

SHELL- PAYLOAD 

596.4 

AFT  B/H-P/L 

14 

1102.1 

AFT  COVER 

634.2 

AFT  B/H-P/L 

15 

1206.5 

AFT  COVER 

692.9 

AFT  B/H-P/L 

16 

1638.4 

SHELL 

738.4 

P/L-AFT  B/H-SHELL 

17 

1755.1 

AFT  BULKHEAD 

772.1 

PAYLOAD 

18 

1815.6 

PAYLOAD 

900.3 

SHELL-F/C-P/L 

19 

2107.1 

AFT  COVER 

1059.8 

P/L-SHELL 

20 

— 

1324.3 

SHELL 

21 

— 

1363.0 

SHELL-F/C 

22 

— 

1626.0 

SHELL-AFT  B/H 

23 

— 

1746.1 

AFT  BULKHEAD 

24 

2065.3 

AFT  BULKHEAD 

*P/L  - PAYLOAD 

A/C  - AFT  COVER 
B/H  * BULKHEAD 
F/C  - FORWARD  CCMP. 
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TABLE  A-H.  SUMMARY  OF  TRAJECTORY  AND  MAXIMUM  AERO- 
ACOUSTIC  CONDITIONS  FOR  BALLISTIC  VEHICLE 


VALUES  AT 
R/V  AFT  CONE 

AERODYNAMIC 

6*  (IN) 

.112 

PARAMETERS 

Ufi  (IN/SEC) 

1.7  x 105 

(q<B=  MAX) 

M 

e 

9.23 

qe  (PSI) 

2130 

AEROACOUSTIC 

ep(0)  (PSI2/HZ) 

1.4  x 10"4 

ENVIRONMENTS 

OASPL 

188 

(dB) 
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* Responses  listed  in  descending  order  by  magnitude. 
**  Flight  location  as  Identified  in  sketch 
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FIGURE  A-l.  SUMMARY  OF  VIBRATION  RESPONSE  PREDICTION  TECHNIQUE 
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